Amplified tree structure technology for fiber optic sensor arrays

ABSTRACT

A folded Sagnac fiber optic sensor array uses a common delay path to reduce distributed pickup in downlead fibers. The sensor array is used to detect perturbations (e.g., acoustic waves in water). By basing the folded Sagnac sensor array on operating principles similar to the Sagnac interferometer the sensor array has a stable bias point, has reduced phase noise, and allows a broadband signal source to be used rather than requiring a more expensive narrowline laser. A large number of sensors can be multiplexed into the architecture of the sensor array. In one particularly preferred embodiment, the sensor array is configured with couplers and amplifiers selected to direct substantially equal portions of an input light signal to each sensor in the array. In another particularly preferred embodiment, first and second delay paths are provided for light at first and second wavelengths to increase the dynamic range of the sensors.

RELATED APPLICATIONS

This application is a continuation-in-part of U.S. Patent applicationSer. No. 09/912,074, filed on Jul. 24, 2001 (now U.S. Pat. No.6,529,444, issued on Mar. 4, 2003, which is a continuation of U.S. Patenapplication Ser. No. 09/520,277, filed on Mar. 6, 2000 (now U.S. Pat.No. 6,278,657, issued on Aug. 21, 2001), which is a continuation-in-partof U.S. Patent application Ser. No. 09/253,203, filed on Feb. 19, 1999(now U.S. Patent No. 6,034,924, issued on Mar. 7, 2000), which claimsthe benefit of priority under 35 U.S.C. §119(e) of U.S. ProvisionalApplication No. 60/080,677, filed on Apr. 3, 1998. This application alsoclaims the benefit of priority under 35 U.S.C. §119(e) of U.S.Provisional Application No. 60/311,749, filed on Aug. 10, 2001, U.S.Provisional Application No. 60/311,867, filed on Aug. 13, 2001, and U.S.Provisional Application No. 60/355,297, filed on Feb. 7, 2002. Theforegoing applications are incorporated by reference herein.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention is in the field of fiber optic acoustic sensorarrays wherein light is propagated in the arrays and the effects ofacoustic signals on the light returning from the arrays are analyzed todetermine the characteristics of the acoustic signals.

2. Description of the Related Art

Fiber optic based acoustic sensors are promising alternatives toconventional electronic sensors. Included among their advantages are ahigh sensitivity, large dynamic range, light weight, and compact size.The ability to easily multiplex a large number of fiber optic sensorsonto common busses also makes fiber optic sensors attractive forlarge-scale arrays. The recent successful incorporation of multiplesmall-gain erbium doped fiber amplifiers (EDFAs) into a fiber opticsensor array to increase the number of sensors that can be supported bya single fiber pair has made large-scale fiber optic sensor arrays evenmore competitive.

For acoustic detection, the fiber optic sensor of choice has been theMach-Zehnder interferometric sensor. In any interferometric sensor,phase modulation is mapped into an intensity modulation through a raisedcosine function. Because of this nonlinear transfer function, asinusoidal phase modulation will generate higher order harmonics. Aninterferometer biased at quadrature (interfering beams π/2 out of phase)has a maximized response at the first order harmonic and a minimizedresponse at the second order harmonic. For this reason, quadrature isthe preferred bias point. As the bias point drifts away from quadrature(for example, due to external temperature changes), the response at thefirst order harmonic decreases and the response at the second orderharmonic increases. When the interferometer is biased at 0 or π out ofphase, the first order harmonic disappears completely. This decreasedresponse at the first order harmonic (resulting from the bias pointsaway from quadrature) is referred to as signal fading.

Because Mach-Zehnder interferometric sensors have an unstable biaspoint, they are especially susceptible to the signal fading problem justmentioned. In order to overcome signal fading, a demodulation of thereturned signal is required. The typical demodulation technique is thePhase-Generated Carrier (PGC) scheme, which requires a path-mismatchedMach-Zehnder interferometric sensor. (See, for example, AnthonyDandridge, et al., Multiplexing of Interferometric Sensors Using PhaseCarrier Techniques, Journal of Lightwave Technology, Vol. LT-5, No. 7,July 1987, pp. 947-952.) This path imbalance also causes the conversionof laser phase noise to intensity noise, which limits the performance ofthe Mach-Zehnder interferometric sensor arrays at low frequencies andplaces stringent requirements on the linewidth of the source. Thisnarrow linewidth requirement has slowed the development of amplifiedMach-Zehnder interferometric sensor arrays at 1.55 μm.

The Sagnac interferometer has found widespread use in the fiber opticgyroscopes. (See, for example, B. Culshaw, et al., Fibre opticgyroscopes, Journal of Physics E (Scientific Instruments), Vol. 16, No.1, 1983, pp. 5-15.) It has been proposed that the Sagnac interferometercould be used to detect acoustic waves. (See, for example, E. Udd,Fiber-optic acoustic sensor based on the Sagnac interferometer,Proceedings of the SPIE-The International Society for OpticalEngineering, Vol. 425, 1983, pp. 90-91; Kjell Kråkenes, et al., Sagnacinterferometer for underwater sound detection: noise properties, OPTICSLETTERS, Vol. 14, No. 20, Oct. 15, 1989, pp. 1152-1145; and SverreKnudsen, et al., An Ultrasonic Fiber-Optic Hydrophone Incorporating aPush-Pull Transducer in a Sagnac Interferometer, JOURNAL OF LIGHTWAVETECHNOLOGY, Vol. 12, No. 9, September 1994, pp. 1696-1700.) Because ofits common-path design, the Sagnac interferometer is reciprocal andtherefore has a stable bias point, which eliminates signal fading andprevents the conversion of source phase noise into intensity noise.Therefore, the Sagnac interferometer is immune to the phase noise whichlimits the Mach-Zehnder interferometric sensors at low frequencies.

SUMMARY OF THE INVENTION

One aspect of the present invention is an acoustic sensor that comprisesa source of light pulses, a first coupler, a polarization dependentsecond coupler, an optical delay path and at least one detector. Thefirst coupler couples the light pulses to a first optical path having afirst optical length and to an array of sensors. The array of sensorscomprises at least a first sensor. The first sensor is in a secondoptical path having a second optical length different from the firstoptical length. The polarization dependent second coupler couples lightpulses received from the first optical path in a first polarization tothe optical delay path and couples light pulses received from the arrayin a second polarization to the optical delay path. The light pulsescoupled to the optical delay path in the first polarization return fromthe optical delay path to the second coupler in the second polarization.The light pulses coupled to the optical delay path in the secondpolarization return from the optical delay path to the second coupler inthe first polarization. The second coupler couples the light pulsesreturning to the second coupler from the optical delay path in the firstpolarization to the first optical path to propagate therein to the firstcoupler. The second coupler couples light pulses returning to the secondcoupler from the optical delay path in the second polarization to thearray to propagate therein to the first coupler. The first couplercombines the light pulses from the first optical path and the lightpulses from the array to cause light pulses traveling equal distancesthrough the first optical path and the array to interfere and togenerate a detectable output signal. The detectable output signal variesin response to acoustic energy impinging on the first sensor. Thedetector detects the detectable output signals to generate a detectoroutput signal responsive to variations in the detectable output signalfrom the first coupler. Preferably, the array includes a second sensor.The second sensor is in a third optical path having a third opticallength different from the first optical length and the second opticallength. Also preferably, the polarization dependent second couplercomprises a polarization beam splitter. In preferred embodiments, theoptical delay path comprises a length of optical waveguide and apolarization rotating reflector. The reflector causes light incident onthe reflector in the first polarization to be reflected as light in thesecond polarization, and causes light incident on the reflector in thesecond polarization to be reflected as light in the first polarization,the reflector advantageously comprises a Faraday rotating mirror. Inparticularly preferred embodiments, the first optical path includes anon-reciprocal phase shifter which causes light propagating through thefirst optical path in a first direction and light propagating throughthe first optical path in a second direction to experience a relativephase shift such that light combined in the first coupler has a phasebias. Preferably, in such embodiments, a third optical path ispositioned in parallel with the first optical path. One of the firstoptical path and the third optical path includes an optical delay tocause the first optical path to have an optical path length differentfrom an optical path length of the third optical path, such that lightpropagating through the first optical path has a propagation timedifferent from a propagation time of light propagating through thesecond optical path to thereby time multiplex the light pulses.Preferably, the non-reciprocal phase shifter comprises a first Faradayrotator, a quarter-wave plate and a second Faraday rotator, the firstFaraday rotator. The quarter-wave plate and the second Faraday rotatorare positioned such that light propagating in the first direction passesthrough the first Faraday rotator, then through the quarter-wave plate,and then through the second Faraday rotator, and such that lightpropagating in the second direction passes through the second Faradayrotator, then through the quarter-wave plate, and then through the firstFaraday rotator. Alternatively, the non-reciprocal phase shiftercomprises a first quarter-wave plate, a Faraday rotator, and a secondquarter-wave plate. The first quarter-wave plate, the Faraday rotator,and the second quarter-wave plate are positioned such that lightpropagating in the first direction passes through the first quarter-waveplate, then through the Faraday rotator, and then through the secondquarter-wave plate, and such that light propagating in the seconddirection passes through the second quarter-wave plate, then through theFaraday rotator, and then through the first quarter-wave plate.

Another aspect of the present invention is an acoustic sensor thatcomprises a source of input light pulses, an array of optical sensors;an optical delay path, an optical detector system; and an input/outputsystem. The input/output system receives the input light pulses anddirects a first portion of each light pulse having a first polarizationthrough the array of optical sensors in a first direction, then throughthe optical delay path, and then to the optical detector system. Theinput/output system directs a second portion of each light pulse in asecond polarization orthogonal to the first polarization through theoptical delay path, then through the optical sensor array in a seconddirection, and then to the optical detector system. The optical detectorsystem receives the light pulses in the first and second polarizationsand detects changes in the light pulses caused by perturbations in theoptical sensors.

Another aspect of the present invention is a method of detectingacoustic signals. The method comprises generating an input light signaland coupling the input light signal to at least first and secondpropagation paths to propagate in respective first directions therein.The first and second propagation paths have respective first and secondoptical lengths. The first and second propagation paths outputrespective first and second output light portions. The first and secondoutput light portions are output from the first and second propagationpaths at differing times in accordance with differences in the first andsecond optical path lengths. The second output light portion ismodulated by an acoustic signal impinging on the second propagationpath. The first light portion is coupled to a delay path in a firstpolarization, and the second light portion is coupled to the delay pathin a second polarization. The delay path outputs a first delayed lightportion corresponding to the first output light portion. The firstdelayed light portion has the second polarization. The delay pathoutputs a second delayed light portion corresponding to the secondoutput light portion. The second delayed light portion has the firstpolarization. The first and second delayed light portions are coupled tothe first and second propagation paths to propagate therein inrespective second directions opposite the respective first directions.The first propagation path outputs a first set of return light portions.The first set of return light portions comprise a respective returnlight portion for each of the first and second delayed light portions.The second propagation path outputs a second set of return lightportions. The second set of return light portions comprise a respectivereturn light portion for each of the first and second delayed lightportions. The first and second sets of return light portions are coupledto at least one detector. The return light portions in the first andsecond sets of return light portions result from output light portionsand delayed light portions which travel identical optical path lengthsand interfere to generate detectable output signals. The methodselectively detects the detectable output signals to detect only outputsignals resulting from interference of light portions which propagatedin the first propagation path in either the first direction or thesecond direction. The detectable output signals vary in response to theacoustic signal impinging on the second propagation path.

Another aspect of the present invention is a sensor that comprises asource of light and a first coupler that couples light to a common pathand to a sensing array to propagate in respective first directionstherein. The sensing array comprises a plurality of sensing paths. Apolarization dependent second coupler couples light from the common pathand from the sensing array to a delay path. The second coupler couplesonly light in a first polarization from the common path to the delaypath. The second coupler couples only light in a second polarizationfrom the sensing array to the delay path. The delay path rotates lightin the first polarization to the second polarization and rotates lightin the second polarization to the first polarization. The second couplerfurther couples light from the delay path in the first polarization tothe common path and couples light from the delay path in the secondpolarization to the sensing array to propagate in respective seconddirections therein to the first coupler. The first coupler providesoutput light responsive to the light propagating in the respectivesecond directions. A detector receives the output light from the firstcoupler and generates an output signal responsive to interference oflight in the first coupler. Preferably, the delay path comprises alength of optical fiber and a polarization rotating reflector. Thelength of optical fiber is selected to provide an optical delay time.The light propagates through the optical fiber from the second couplerto the reflector. The reflector reflects light into the optical fiber topropagate through the optical fiber to the second coupler. The reflectorfurther rotates light incident in the first polarization to the secondpolarization and rotates light incident in the second polarization tothe first polarization. Preferably, the reflector comprises a Faradayrotating mirror. Also preferably, the polarization dependent secondcoupler comprises a polarization beam splitter positioned so that thedelay path receives the light from a port of the polarization beamsplitter and returns light to the port of the polarization beamsplitter.

Another aspect of the present invention is a sensor array that comprisesa source of light and a first coupler that receives the light from thesource. The first coupler couples a first portion of the light to afirst coupler port and couples a second portion of the light to a secondcoupler port. An interferometric loop has a first end coupled to thefirst coupler port to receive the first portion of the light and has asecond end coupled to the second coupler port to receive the secondportion of the light. The interferometric loop propagates the firstportion of the light in a first direction to the second coupler port andpropagates the second portion of the light in a second directionopposite the first direction to the first coupler port. Theinterferometric loop comprises a plurality of sensors coupled inparallel between the first end of the interferometric loop and thesecond end of the interferometric loop. The sensors perturb lightpassing through the sensors in response to a sensed parameter (e.g.,acoustic signals). A first plurality of couplers distribute the firstportion of the light approximately equally to each of the sensors,collect the second portion of the light from each of the sensors, andpropagate the collected light to the first end of the interferometricloop. A second plurality of couplers distribute the second portion ofthe light approximately equally to each of the sensors, collect thefirst portion of the light from each of the sensors, and propagate thecollected light to the second end of the interferometric loop. At leastone first amplifier is coupled between the first end of theinterferometric loop and the first plurality of couplers. At least onesecond amplifier is coupled between the second end of theinterferometric loop and the second plurality of couplers. A pluralityof delay portions are connected between the first and second ends of theinterferometric loops and the sensors. The delay portions have delaysselected so that the light passing through each sensor is delayed by adifferent amount than the light passing through the other sensors.Preferably, the first plurality of couplers further comprise a firstdistribution coupler that receives the first portion of the light fromthe first amplifier and distributes the first portion of the light to afirst plurality of internal amplifiers that are coupled between thefirst distribution coupler and a first plurality of internal couplers.The first distribution coupler collects the second portion of the lightfrom the first plurality of internal amplifiers and propagates thesecond portion of the light to the first amplifier. Also preferably, thesecond plurality of couplers further comprise a second distributioncoupler that receives the second portion of the light from the secondamplifier and distributes the second portion of the light to a secondplurality of internal amplifiers that are coupled between the seconddistribution coupler and a second plurality of internal couplers. Thesecond distribution coupler collects the first portion of the light fromthe second plurality of internal amplifiers and propagates the firstportion of the light to the second amplifier. The first plurality ofinternal couplers distribute the first portion of the light to theplurality of sensors and collect the second portion of the light fromthe plurality of sensors. The second plurality of internal couplersdistribute the second portion of the light to the plurality of sensorsand collect the first portion of the light from the plurality ofsensors. Advantageously, the source of light is a broadband source, suchas, for example, a superfluorescent fiber source. Also advantageously,the first and second amplifiers and the first and second plurality ofinternal amplifiers are erbium-doped fiber amplifiers, and the first andsecond distribution couplers and the first and second pluralities ofinternal couplers comprise 4×4 couplers.

Another aspect of the present invention is a method of sensing aparameter that comprises propagating light from a source of lightthrough an interferometric loop such that approximately equal portionsof the light counterpropagate in first and second directions in theloop. The light propagating in the first direction of theinterferometric loop is amplified and coupled into a plurality ofsensors such that approximately equal portions of the light propagatingin the first direction are passed through each of the sensors. The lightpropagating in the second direction of the interferometric loop isamplified and coupled into the plurality of sensors such thatapproximately equal portions of the light propagating in the seconddirection are passed through each of the sensors. The light propagatingin the first direction is caused to interfere with the light propagatingin the second direction to generate a plurality of output signalsresponsive to light passing through each sensor in the first and seconddirections. Each of the sensors perturbs light passing therethrough inresponse to a sensed parameter (e.g., an acoustic signal), and each ofthe sensors has a unique optical path length such that the lightpropagating in the first direction interferes with the light propagatingin the second direction at a unique time.

Another aspect of the present invention is a sensor system that sensesperturbations over first and second dynamic ranges. The sensor systemcomprises a source of input light pulses at a first wavelength and asource of input light pulses at a second wavelength. The system includesan array of sensors, a first optical delay path at the first wavelength,and a second optical delay path at the second wavelength. A firstdetection system is responsive to light at the first wavelength, and asecond detection system is responsive to light at the second wavelength.An input/output system receives the input light pulses at the firstwavelength and the second wavelength. The input/output system directs afirst portion of each light pulse at the first wavelength having a firstpolarization through the array of sensors in a first direction, thenthrough the first optical delay path, and then to the first detectionsystem. The input/output system directs a second portion of each lightpulse at the first wavelength in a second polarization orthogonal to thefirst polarization through the first optical delay path, then throughthe array of sensors in a second direction, and then to the firstdetection system. The first detection system detects variations inreceived light caused by perturbations varying over the first dynamicrange. The input/output system directs a first portion of each lightpulse at the second wavelength having a first polarization through thearray of sensors in a first direction, then through the second opticaldelay path, and then to the second detection system. The input/outputsystem directs a second portion of each light pulse at the secondwavelength in a second polarization orthogonal to the first polarizationthrough the second optical delay path, then through the array of sensorsin a second direction, and then to the second detection system. Thesecond detection system detects variations in received light caused byperturbations varying over the second dynamic range.

Another aspect of the present invention is an acoustic sensor systemthat senses acoustic signals over first and second dynamic ranges. Theacoustic sensor system comprises a source of input light pulses at afirst wavelength and a source of input light pulses at a secondwavelength. The acoustic sensor system further includes an array ofacoustic sensors, a first optical delay path at the first wavelength,and a second optical delay path at the second wavelength. A firstdetection system is responsive to light at the first wavelength. Asecond detection system is responsive to light at the second wavelength.An input/output system receives the input light pulses at the firstwavelength and the second wavelength. The input/output system directs afirst portion of each light pulse at the first wavelength having a firstpolarization through the array of acoustic sensors in a first direction,then through the first optical delay path, and then to the firstdetection system. The input/output system directs a second portion ofeach light pulse at the first wavelength in a second polarizationorthogonal to the first polarization through the first optical delaypath, then through the array of acoustic sensors in a second direction,and then to the first detection system. The first detection systemdetects variations in received light caused by acoustic signals varyingover the first dynamic range. The input/output system directs a firstportion of each light pulse at the second wavelength having a firstpolarization through the array of acoustic sensors in a first direction,then through the second optical delay path, and then to the seconddetection system. The input/output system directs a second portion ofeach light pulse at the second wavelength in a second polarizationorthogonal to the first polarization through the second optical delaypath, then through the array of acoustic sensors in a second direction,and then to the second detection system. The second detection systemdetects variations in received light caused by acoustic signals varyingover the second dynamic range.

Another aspect of the present invention is a method of sensingperturbations. The method comprises inputting light pulses of a firstwavelength into an array of sensors that includes a first optical delaypath at the first wavelength. Light pulses of a second wavelength arealso input into the array of sensors. The array of sensors includes asecond optical delay path at the second wavelength. The second opticaldelay path has a different optical length than the first optical delaypath. A first portion of each light pulse at the first wavelength havinga first polarization is directed through the array of sensors in a firstdirection, then through the first optical delay path. A second portionof each light pulse at the first wavelength in a second polarizationorthogonal to the first polarization is directed through the firstoptical delay path, then through the array of sensors in a seconddirection. Variations in the first and second portions of each lightpulse at the first wavelength caused by perturbations varying over afirst dynamic range are detected. A first portion of each light pulse atthe second wavelength having a first polarization is directed throughthe array of sensors in a first direction, then through the secondoptical delay path. A second portion of each light pulse at the secondwavelength in a second polarization orthogonal to the first polarizationis directed through the second optical delay path, then through thearray of sensors in a second direction. Variations in the first andsecond portions of each light pulse at the second wavelength caused byperturbations varying over a second dynamic range are detected. Inparticular embodiments of the method, the perturbations are acousticsignals.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will be described below in connection with theaccompanying drawing figures in which:

FIG. 1 illustrates an exemplary Sagnac interferometer having a singlesensing loop;

FIG. 2 illustrates a Sagnac sensor array in accordance with the presentinvention wherein each rung of a sensor array forms an additional Sagnacinterferometer;

FIG. 3 illustrates a Sagnac sensor array which includes erbium-dopedfiber amplifiers to regenerate signal power lost to coupling anddissipative losses;

FIG. 4 illustrates a graph of the frequency response of a Sagnacinterferometer in accordance with present invention compared with thethree dominant ocean floor noises;

FIG. 5 illustrates graphs of the maximum and minimum acoustic signaldetectable by a Mach-Zehnder interferometer and detectable by a Sagnacinterferometer in accordance with the present invention, showing therelatively constant dynamic range of a Sagnac interferometer over a widerange of frequencies;

FIG. 6 illustrates graphs of the minimum detectable acoustic signalversus frequency for three Sagnac interferometer configurations havingdifferent lengths of fiber in the hydrophone and the delay loop;

FIG. 7 illustrates a Sagnac interferometer in accordance with thepresent invention which includes an additional delay loop to increasethe dynamic range of the interferometer;

FIG. 8 illustrates a graph of the dynamic range provided by theinterferometer of FIG. 7;

FIG. 9A illustrates the positioning of the delay loop of theinterferometer in the dry end of a sensor array system;

FIG. 9B illustrates the positioning of the delay loop of theinterferometer in the wet end of a sensor array system;

FIG. 10 illustrates the Sagnac interferometer of FIG. 9B withannotations showing the lengths used in calculations of the effects ofphase modulation;

FIG. 11 illustrates a technique for winding the delay loop so as toreduce the effects of the acoustic wave upon the delay loop;

FIG. 12 illustrates a Sagnac interferometer in accordance with thepresent invention which includes empty rungs which detect distributedpick-up noise which can be subtracted from the signals generated by thesensors;

FIG. 13 illustrates a Sagnac interferometer in accordance with thepresent invention which includes a depolarizer to reduce the effects ofpolarization induced fading;

FIG. 14 illustrates a Sagnac interferometer which utilizes frequencydivisional multiplexing;

FIG. 15 illustrates a graph which shows the generation of the beatsignals between the delayed modulation signal and the returning sensorsignals in the interferometer of FIG. 14;

FIG. 16 illustrates a Sagnac interferometer which utilizes code divisionmultiplexing;

FIG. 17 illustrates the architecture of a folded Sagnac acoustic fibersensor array;

FIG. 18 illustrates a graph of the number of returned pulses per timeinterval, showing the separation in time of signal pulses and noisepulses;

FIG. 19 illustrates a folded Sagnac acoustic fiber sensor array having asecond delay loop to provide extended dynamic range;

FIG. 20 illustrates a folded Sagnac acoustic fiber sensor array having aphase modulator and nulling circuitry in place of the reflector in FIG.17;

FIG. 21 illustrates a further alternative embodiment of FIG. 19 in whichthe two delay loops are connected to different ports of the coupler;

FIG. 22 illustrates an alternative embodiment of a fiber optic acousticsensor array system using a Faraday rotating mirror;

FIGS. 23A, 23B and 23C illustrate further alternative embodiments of afiber optic acoustic sensor array which utilize an unpolarized lightsource in combination with a depolarizer, a polarization beam splitterand a Faraday rotating mirror;

FIG. 24 illustrates an alternative embodiment of a folded fiber opticacoustic sensor array which utilizes an unpolarized light source incombination with an optical circulator, a 2×2 coupler, and anon-reciprocal phase shifter;

FIG. 25 illustrates an alternative embodiment of a folded fiber opticacoustic sensor array similar to FIG. 24 in which the depolarizer islocated in the second array input/output fiber;

FIG. 26 illustrates a first preferred embodiment of the non-reciprocalπ/2 phase shifter in FIGS. 24 and 25, which illustrates the effect onthe polarization of the light propagating in a first direction throughthe phase shifter;

FIG. 27 illustrates the effect on the polarization of the lightpropagating in a second (opposite) direction through the phase shifterof FIG. 26;

FIG. 28 illustrates an alternative preferred embodiment of thenon-reciprocal π/2 phase shifter in FIGS. 24 and 25, which illustratesthe effect on the polarization of the light propagating in a firstdirection through the phase shifter;

FIG. 29 illustrates the effect on the polarization of the lightpropagating in a second (opposite) direction through the phase shifterof FIG. 28;

FIG. 30 illustrates a further alternative embodiment of a folded fiberoptic acoustic sensor array, which utilizes polarization-based biasingfor multiple detectors, wherein each detector has a bias point which canbe set independently of the bias points of the other detectors;

FIG. 31 illustrates an alternative embodiment of a folded fiber opticacoustic sensor array similar to FIG. 30 in which the depolarizer islocated in the second array input/output fiber;

FIG. 32 illustrates an alternative embodiment of a folded fiber opticacoustic sensor array similar to FIG. 30 in which an optical circulatorreplaces the 2×2 coupler;

FIG. 33 illustrates an alternative embodiment of a folded fiber opticacoustic sensor array similar to FIG. 32 in which the depolarizer islocated in the second array input/output fiber;

FIG. 34 illustrates a further alternative embodiment of a folded Sagnacsensor array, which includes a combined input/output subsystem;

FIG. 35 illustrates an alternative embodiment of a folded fiber opticacoustic sensor array similar to FIG. 34 in which the depolarizer islocated in the second array input/output fiber;

FIG. 36 illustrates a further alternative embodiment of a folded fiberoptic acoustic sensor array similar to FIGS. 34 and 35 in which thedetectors are coupled to the input/output subsystem by optical fibers topermit the detectors to be located remotely;

FIG. 37 illustrates an alternative embodiment of a 16-sensor array foruse in combination with the front end systems described in FIGS. 30-36;

FIG. 38, comprising FIGS. 38A and 38B, illustrates an alternative sensorarray that includes 16 sensors in an amplified tree-structure telemetryusing a polarization-based front end similar to the system describedabove in connection with FIGS. 30-36;

FIG. 39A illustrates the returned pulse train measured at one of thedetectors in FIG. 38A for a pulse width of 50 nanoseconds and arepetition frequency of 0.942 MHz (1.06-microsecond spacing;

FIG. 39B illustrate the returned pulse train measured at one of thedetectors in FIG. 38A for a pulse width of 50 nanoseconds and arepetition frequency of 1.042 MHz (0.96-microsecond spacing);

FIG. 40A illustrates an oscilloscope trace representing the measuredpower detected from one of the sensors in FIG. 38B when a peal-to-peakphase modulation of greater than π is induced in a PZT transducer inseries with the sensor;

FIG. 40B illustrates the results of measuring the visibility of thesensor in accordance with FIG. 40A for eight random settings of a commonpolarization controller positioned in one of the input/output fibers ofthe array in FIG. 38B;

FIG. 41 illustrates the results of measurement of the detected responsefrom a single detector configuration when a fiber wrapped PZT is used tosimulate an acoustic signal and an amplitude modulation at 3.6 kHz wasplaced on the optical signal from the source using a lithium niobatemodulator in the signal path from the source to the collimator in FIG.38A;

FIGS. 42A, 42B and 42C illustrate the effects of source power on themeasured noise at the detector (receiver) as a function of sensor countfor DC source currents of 1.4 μA, 4.05 μA and 9.89 μA, respectively;

FIG. 43A illustrates a plot of the minimum detectable acoustic signal ina 1-second integration time and also illustrates a plot of the maximumdetectable signal defined as that which produces a 1 rad/{square rootover (Hz)} phase signal for a delay time through the delay loop of 50microseconds;

FIG. 43B illustrates a plot of the minimum detectable acoustic signal ina 1-second integration time and also illustrates a plot of the maximumdetectable signal defined as that which produces a 1 rad/{square rootover (Hz)} phase signal for the same parameters as in FIG. 43A but for adelay time through the delay loop of 500 nanoseconds;

FIG. 44A illustrates the phase variation in a sensor caused an acousticsignal applied to one of the sensors of FIG. 38B;

FIGS. 44B and 44C are plots of the induced phase difference due to thephase signal of FIG. 44A for the signal which travels through the longerdelay loop (FIG. 44B) and the shorter delay loop (FIG. 44C);

FIGS. 45A illustrates a front end for an acoustic sensor array, which issimilar to the front end of FIG. 38A except that the front end of FIG.45A includes a first signal source that outputs an first sequence ofoptical signal pulses at a first wavelength λ₁ and includes a secondsignal source that outputs a second sequence of optical signal pulses ata second wavelength λ₂, and includes a first delay loop for the opticalsignal pulses at the first wavelength and a second delay loop for theoptical signal pulses at the second wavelength;

FIG. 45B illustrates an acoustic sensor array connected to the front endof FIG. 45A;

FIG. 46A is a plot of the signal from a given sensor that passes throughthe longer delay loop of FIG. 45A;

FIG. 46B illustrates the use of linear extrapolation to predict a fringecount of β(2τ) from the plot of FIG. 46A; and

FIG. 47 is a plot of the maximum value of A before fringe count errorsoccur as a function of frequency for T_(d) ^(λ) ^(₁) =100 microseconds,τ=50 microseconds, and for T_(d) ^(λ) ^(₂) =100 nanoseconds.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The present invention is described below in connection with an array ofacoustic sensors (e.g., hydrophones) in a Sagnac loop. Before describingthe preferred embodiments, a brief review of the operation of a singleloop Sagnac acoustic sensor is provided.

Single Loop Sagnac Acoustic Sensor

A simple Sagnac-based acoustic sensor 100 is shown in FIG. 1. The Sagnacloop is divided into two portions, a delay loop 102 and a hydrophone104. The delay loop 102 is simply a large length of fiber, typicallygreater than 1 km. The hydrophone 104 is a portion of fiber in which anacoustic wave is transformed into a phase modulation of an opticalsignal propagating through the fiber. A high responsivity to acousticwaves is typically accomplished by selecting optimized coatings for thesection of fiber in the hydrophone 104, and wrapping the fiber around amandrel of suitable composition. (See, for example, J. A. Bucaro, etal., Optical fibre sensor coatings, Optical Fiber Sensors, Proceedingsof the NATO Advanced Study Institute, 1986, pp. 321-338.) The length offiber wrapped around the hydrophone 104 is typically 10 meters to 100meters. Light from a source 110, such as, for example, asuperfluorescent fiber source (SFS), is split into clockwise (CW) andcounter-clockwise (CCW) beams by a 3×3 coupler 112. The operation of the3×3 coupler 112 is well-known and is described, for example, in Sang K.Sheem, Fiber-optic gyroscope with [3×3]directional coupler, AppliedPhysics Letters, Vol. 37, No. 10, Nov. 15, 1980, pp. 869-871.

Although described herein as using a 3×3 coupler 112, other couplers(e.g., a 2×2 coupler, a 4×4 coupler, etc.) can be used with alternativeembodiments of the present invention. For example, to use a 2×2 coupler,both ports of one side are used to create the Sagnac interferometer. Oneport of the other side is a detection port. The remaining port is usedto launch light into the array and can also be used as a detection portif a coupler or circulator is employed (in a similar manner as is donewith fiber optic gyroscopes). In general, any n×m coupler can beemployed by using two ports of one side of the coupler to create theSagnac interferometer and using the ports on the other side of thecoupler as detection ports, launching ports, or both.

After splitting, the CW beam travels through the delay loop 102 firstand then through the hydrophone 104, while the CCW beam travels throughthe hydrophone 104 first and then through the delay loop 102. During atime delay T_(delay) between a time when the CW beam travels through thehydrophone 104 and a time when the CCW beam travels through thehydrophone 104, the acoustic signal and likewise the acousticallyinduced phase modulation in the hydrophone 104 changes. This change inphase modulation is mapped into a phase difference between thecounter-propagating beams, which is converted into an intensitymodulation when the beams recombine at the 3×3 coupler 112. Thisintensity modulation is then detected by a first detector 120 and asecond detector 122 or by only one of the two detectors.

More explicitly, if an acoustic signal induces a phase modulation φ_(h)cos(Ωt) in the fiber of the hydrophone 104, the resulting phasemodulation between the interfering beams at the hydrophone 104,φ_(int)(t), is given by: $\begin{matrix}\begin{matrix}{{\varphi_{int}(t)} = \quad {{\varphi_{h}{\cos \left( {\Omega \quad t} \right)}} - {\varphi_{h}{\cos \left( {\Omega \left( {t + T_{delay}} \right)} \right)}}}} \\{= \quad {2\varphi_{h}{\sin \left( \frac{\Omega \cdot T_{delay}}{2} \right)}{\sin \left( {{\Omega \quad t} + \frac{\Omega \cdot T_{delay}}{2}} \right)}}}\end{matrix} & (1)\end{matrix}$

where T_(delay) is the travel time through the delay loop. Thus,φ_(int)(t) is a function of the hydrophone modulation oh and the productof the acoustic modulation frequency, Ω, with the loop delay, T_(delay).This differs from a Mach-Zehnder interferometric sensor in whichφ_(int)(t) is a function of only the hydrophone modulation φ_(h).Maximum sensitivity is achieved in the Sagnac loop acoustic sensor whenthe product of the acoustic frequency, Ω, and the time delay, T_(delay),is an odd multiple of π (maximum value of the first sine term inEquation 1). The acoustic frequency which makes this product π is calledthe proper frequency of the loop, which is the lowest frequency at whichmaximum sensitivity is achieved. Most underwater sensing applicationsare concerned with the detection of acoustic frequencies below 10 kHz.For the proper loop frequency to be less than 10 kHz, a delay time of atleast 50 microseconds and therefore a delay loop length of at least 10km is required. Thus, the Sagnac acoustic sensor 100 requires a largeamount of fiber for the detection of low acoustic frequencies (<10 kHz).

The common-path design inherent to the Sagnac interferometer has manyadvantages over a Mach-Zehnder interferometer in addition to the stablebias point and elimination of phase noise already mentioned. A Sagnacinterferometer allows the use of a short-coherence length, broadbandsource, such as a superfluorescent fiber source (SFS), an example of anamplified spontaneous emission (ASE) source. Such sources areinexpensive and can readily provide high powers. It has been shown thatthe use of the 3×3 coupler passively biases the Sagnac acoustic sensornear quadrature. (See, Sang K. Sheem, Fiber-optic gyroscope with [3×3]directional coupler, Applied Physics Letters, Vol. 37, No. 10, Nov. 15,1980, pp. 868-871; and H. Poisel, et al., Low-cost fibre-opticgyroscope, Electronics Letters, Vol. 26, No. 1, Jan. 4, 1990, pp.69-70.) By subtracting the signals from the two detection ports of the3×3 coupler, the source excess noise, which is the limiting noise sourceof SFS sources, can be subtracted while phase-modulation inducedintensity variations due to the hydrophone are added. This allows aSagnac interferometer to approach near shot-noise limited performance.(See, Kjell Kråkenes, et al., Sagnac interferometer for underwater sounddetection: noise properties, OPTICS LETTERS, Vol. 14, No. 20, Oct. 15,1989, pp. 1152-1145.)

Previous work on Sagnac-based acoustic sensors has been limited to asingle sensor configuration. Because of the inherent advantages of theSagnac interferometer, Applicants have determined that it is desirableto replace the Mach-Zehnder interferometric sensors in a large-scalearray with Sagnac based sensors. Each Sagnac sensor 100 discussed aboverequires many kilometers of fiber, making the insertion of numerous suchsensors into a large-scale array impractical. Research into usingrecirculating delay loops to reduce the fiber length requirement hasproduced sensors which use significantly less fiber but suffer from highnoise due to the incorporation of EDFAs within the recirculating loop.(See, for example, J. T. Kringlebotn, et al., Sagnac InterferometerIncluding A Recirculating Ring With An Erbium-doped Fibre Amplifier, OFS'92 Conference Proceedings, pp. 6-9.) A novel approach for decreasingthe fiber required is described below.

Novel Sensor Array Based on the Sagnac Interferometer

As set forth below, Applicants have discovered a novel system whichreduces the amount of fiber needed for a Sagnac-based large scale arrayby multiplexing multiple sensors onto the same delay loop, producing apractical Sagnac sensor array (SSA). As illustrated in FIG. 2, a Sagnacsensor array 200 in accordance with the present invention includes anarray 210 of hydrophones 212(i) in a ladder configuration which areattached to a single delay loop 214. For example, FIG. 2 shows a Sagnacsensor array 210 having N hydrophones 212(1), 212(2) . . . 212(N) inrespective rungs 216(1), 216(2) . . . 216(N). Each rung 216(i) in theSagnac sensor array 210 comprises a single fiber wrapped around arespective hydrophone 212(i). Every path from a 3×3 coupler 220 throughthe delay loop 214 and array 210 and back to the coupler 220 comprises aseparate Sagnac interferometer. Therefore, for an array of N sensors212, there are N separate Sagnac interferometers, each of which behaveslike the single loop Sagnac sensor 100 shown in FIG. 1. Each Sagnacinterferometer measures the acoustic signal at a separate point inspace, i.e., the location of the hydrophone 212(i). For example, theSagnac interferometer comprising the delay loop 214 and the rung 216(1)measures the acoustic signal at hydrophone 212(1). In addition, eachSagnac interferometer also picks up acoustic signals (e.g., noise)elsewhere in the loop, which noise is advantageously reduced, as will bediscussed below.

The Sagnac sensor array 200 is easiest understood in a time-divisionmultiplexed (TDM) configuration (non-TDM schemes are discussed later). Asource 222 (which may advantageously comprise a conventional pulsedsource or may comprise a cw source with an external modulator) generatesa light pulse which enters the Sagnac loop via a third port of thecoupler 220 and propagates in both the CW and CCW directions asindicated in FIG. 2. Upon reaching the array 210, the CCW pulse is splitinto a train of N separate pulses. At this point, the CW input pulse hasnot yet reached the array 210 and is still a single pulse. When the CWpulse reaches the array 210, it also is split into a train of N pulses.Each pulse in the CW train returns to the 3×3 coupler 220 aftertraveling through a respective rung 216(i) and interferes with the pulsein the CCW train which has traveled the same rung 216(i) in the oppositedirection. Thus, N pulses are detected by a first detector 230 and asecond detector 232, and each pulse comprises the CW and CCW pulses ofone of the N Sagnac loops (i.e., the two pulses which have traveled inopposite directions through the same respective rung 216(i)). Becausethe pulses which travel through different combinations of rungs do nottravel identical optical paths, such pulses are not coincident in timeat the coupler 220, and thus do not interfere with each other at thecoupler 220. The pulse widths should be smaller than the differentialdelay between adjacent sensors so that the pulses from adjacent sensorsdo not overlap.

As illustrated in FIG. 3, small-gain erbium doped fiber amplifiers(EDFAs) 240 are advantageously added to the array portion 210 just asEDFAs have been added to Mach-Zehnder interferometric sensor arrays.(See, for example, Craig W. Hodgson, et al., Optimization of Large-ScaleFiber Sensor Arrays Incorporating Multiple Optical Amplifiers-Part I:Signal-to-Noise Ratio, JOURNAL OF LIGHTWAVE TECHNOLOGY, Vol. 16, No. 2,February 1998, pp. 218-223; Craig W. Hodgson, et al., Optimization ofLarge-Scale Fiber Sensor Arrays Incorporating Multiple OpticalAmplifiers-Part II: Pump Power, JOURNAL OF LIGHTWAVE TECHNOLOGY, Vol.16, No. 2, February 1998, pp. 224-231; Jefferson L. Wagener; et al.,Novel Fiber Sensor Arrays Using Erbium-Doped Fiber Amplifiers, JOURNALOF LIGHTWAVE TECHNOLOGY, Vol. 15, No. 9, September 1997, pp. 1681-1688;and C. W. Hodgson, et al., Large-scale interferometric fiber sensorarrays with multiple optical amplifiers, OPTICS LETTERS, Vol. 22, No.21, Nov. 21, 1997, pp. 1651-1653.) The EDFAs 240 increase the number ofsensors which can be supported by a single array 210 by regenerating thesignal power which is lost to coupling and dissipative losses. The EDFAsare advantageously pumped by one or more pump laser sources 242 via asplitting coupler 244 and via a first wavelength division multiplexing(WDM) coupler 246 and a second WDM coupler 248.

Because it uses the Sagnac architecture, the Sagnac sensor array 200 hasall of the advantages of the single loop Sagnac based sensor 100discussed above. The common-path design eliminates the conversion ofsource phase noise into intensity noise at the interfering coupler 220.The source 222 can be a fiber ASE (amplified spontaneous emission)source (i.e., the SFS discussed above), which provides high powersinexpensively at 1.55 μm. Passive biasing near quadrature is achievablefor all sensors by using the 3×3 coupler 220. Also, the 3×3 coupler 220provides a convenient means to detect two interferometric outputs at thedetectors 230, 232, and to use the outputs of the two detectors tosubtract source excess noise. (See, for example, K. Krakenes, et. al.,Sagnac interferometer for underwater sound detection: noise properties,OPTICS LETTERS, Vol. 14, 1989, pp. 1152-1154, which shows the use of twodetectors in combination with a single Sagnac interferometer.)

The properties of this novel Sagnac sensor array 200 will be discussedmore specifically below followed by a more detailed discussion of thefrequency response and dynamic range which result from the use of aSagnac interferometer. Thereafter, a calculation of the magnitude of thedistributed pick-up from the non-hydrophone fiber loop segments will bedescribed, along with a technique for reducing this pick-up magnitude.Polarization will also be addressed below. New sources of noise whichare introduced by the Sagnac design are then discussed. Finally,multiplexing schemes other than TDM for the Sagnac sensor array arepresented.

Although the present invention is described above with respect to asingle sensor in each rung 216(i) of the array 210, it should beunderstood that each rung 216(i) may advantageously comprise a subarrayhaving multiple sensors, such as are described, for example, in allowedU.S. patent application Ser. No. 08/814,548, filed on Mar. 11, 1997,which is incorporated by reference herein. (See, also, C. W. Hodgson, etal., Large-scale interferometric fiber sensor arrays with multipleoptical amplifiers, Optics Letters, Vol. 22, 1997, pp. 1651-1653; J. L.Wagener, et al., Novel fiber sensor arrays using erbium-doped fiberamplifiers, Journal of Lightwave Technology, Vol. 15, 1997, pp.1681-1688; C. W. Hodgson, et al., Optimization of large-scale fibersensor arrays incorporating multiple optical amplifiers, Part I:signal-to-noise ratio, Journal of Lightwave Technology, Vol. 16, 1998,pp. 218-223; and C. W. Hodgson, et al., Optimization of large-scalefiber sensor arrays incorporating multiple optical amplifiers, Part II:pump power, Journal of Lightwave Technology, Vol. 16, 1998, pp.224-231.)

Frequency Response

As set forth above, the Sagnac sensor has a frequency dependent responsegiven by Equation 1. At frequencies well below the proper frequency ofthe loop, defined as 1/(2·T_(delay)), the minimum detectable acousticsignal scales with the inverse of acoustic frequency. This decreasedacoustic sensitivity at low frequencies has been a major concern for theSagnac acoustic sensor. However, it has been pointed out that thisdecreased sensitivity at low frequencies is fortunately matched by anincreasing ocean noise floor (See, for example, Sverre Knudsen, Ambientand Optical Noise in Fiber-Optic Interferometric Acoustic Sensors,Fiber-Optic Sensors Based on the Michelson and Sagnac Interferometers:Responsivity and Noise Properties, Thesis, Chapter 3, NorwegianUniversity of Science and Technology, 1996, pp. 37-40.) Ideally, itwould be desirable if the minimum detectable acoustic signal of an arrayat a given frequency were to be a constant amount below the ocean noisefloor at that frequency. Thus, the minimum detectable acoustic signalwould also increase at lower frequencies to match the increasing oceannoise floor. The frequency response of the Sagnac sensor array 200 ofthe present invention in fact does provide a good match between theocean noise floor and acoustic sensitivity. This is illustrated in FIG.4, where the minimum detectable acoustic signal for a Sagnac sensorarray is plotted as a curve 250 assuming an optical noise floor of 10μrad/{square root over (Hz)}, a hydrophone phase responsivity of3.2×10⁻⁷ rad/μPa and a delay loop length of 20 km. (The vertical axis isin dB relative to a baseline of 1 μrad/{square root over (Hz)}.) Alsoplotted in FIG. 4 are the ocean noise floors for the three dominantocean noise sources at these frequencies and a resulting sum of thenoise from the three sources. A curve 252 represents the noise fromocean turbulence, earthquakes, volcanic eruptions, and the like. A curve253 represents light shipping noise. A curve 254 represents DSSO(distant shipping and storms) noise. A curve 256 represents the sum ofthe noise floors from the three dominant sources (i.e., the sum of thecurves 252, 253 and 254). (See, for example, Robert J. Urick, The noisebackground of the sea: ambient noise level, Principles of UnderwaterSound, 3rd Ed., Chapter 7, McGraw-Hill, 1983, pp. 202-236.) The minimumdetectable acoustic signal of the Sagnac sensor array 200 increases insuch a way as to provide a nearly constant amount of detectable signalbelow the ocean noise floor at all frequencies below 10 kHz. Thus, thefrequency-dependent response of the Sagnac sensor array 200 does notprohibit low-frequency acoustic detection. The Mach-Zehnder array showsthe same trend as the Sagnac sensor array, namely a decreasingsensitivity towards lower frequencies, but in the Mach-Zehnder array,the decreasing sensitivity is smaller than in the Sagnac-based sensor.

Although both the Mach-Zehnder interferometer and Sagnac sensor array200 have similar frequency-dependent responses, the source of theirfrequency responses is fundamentally different. The increasing minimumdetectable signal in the Mach-Zehnder interferometer sensor array is dueto an increasing optical noise floor. The cause of this increasingoptical noise floor is the phase noise introduced by the path-imbalancedMach-Zehnder interferometer. Thus, although the noise floor is 10μrad/{square root over (Hz)} at 10 kHz, it increases towards lowerfrequencies. In the Sagnac sensor array 200, the increasing minimumdetectable acoustic signal is due to the sin(ΩT_(delay)/2) term inEquation 1, and not to an increasing optical noise floor. The opticalnoise floor remains a constant 10 μrad/{square root over (Hz)} over theentire frequency range.

The significance of this difference can be seen by examining the dynamicrange of the Mach-Zehnder interferometric sensor array and Sagnac sensorarray 200, illustrated in FIG. 5. The dynamic range of a sensor islimited by the minimum and maximum detectable phase shifts. Forinterferometric sensors, the maximum detectable phase shift is limitedby the nonlinear response of the interferometer and the minimumdetectable phase shift by the optical noise floor. Both the Mach-Zehnderinterferometric sensor array and the Sagnac sensor array have maximumdetectable phase shifts which are constant over the acoustic frequencyrange. However, the Sagnac sensor array 200 also has a flat minimumdetectable phase shift because it has a flat optical noise floor, whilethe Mach-Zehnder interferometric sensor array suffers an increasingminimum detectable phase shift due to an increasing optical noise floorcaused by the phase noise introduced by the path imbalancedinterferometer. The Sagnac sensor array 200 thus has a constant dynamicrange at all acoustic frequencies, while the Mach-Zehnderinterferometric sensor array has a decreased dynamic range at lowacoustic frequencies. This is illustrated in FIG. 5, wherein the minimumand maximum detectable acoustic signals (in dB arbitrary units) areplotted for the Sagnac sensor array 200 and a Mach-Zehnderinterferometric sensor array. As shown in FIG. 5, both arrays have anapproximately 100 dB dynamic range above 1 kHz, where phase noise doesnot limit the Mach-Zehnder interferometric sensor array. At 10 Hz, phasenoise dominates the Mach-Zehnder interferometric sensor array, and itsdynamic range is reduced to approximately 74 dB. Meanwhile, the dynamicrange of the Sagnac sensor array 200 remains at approximately 100 dB.

It is interesting to examine the frequency response of the Sagnac sensorarray 200 at frequencies well below the loop proper frequency as afunction of the delay loop length and hydrophone responsivity. At thesefrequencies, the sin(ΩT_(delay)/2) factor in Equation 1 can beapproximated as ΩT_(delay)/2, showing that the responsivity of theSagnac sensor array 200 is proportional to the product of Ω_(h) andT_(delay). Ω_(h) itself is proportional to the amount of fiber in eachhydrophone 212(i), and T_(delay) is proportional to the amount of fiberin the delay loop 214. Thus, the responsivity at frequencies well belowthe loop proper frequency is proportional to the product of thehydrophone fiber length and delay fiber length. FIG. 6 plots the minimumdetectable acoustic signal for several Sagnac sensor arrayconfigurations in which the product of the length of the fiber in eachhydrophone 212(i) and the length of the fiber in the delay loop 214 isconstant, but the relative distribution of fiber between the delay loop214 and each hydrophone 212(i) changes. For example, a curve 260represents the frequency response of a Sagnac sensor array 200 having 45km of fiber in its delay loop 214 and 100 meters of fiber in eachhydrophone 212(i); a curve 262 represents the frequency response of aSagnac sensor array 200 having 30 km of fiber in its delay loop 214 and150 meters of fiber in each hydrophone 212(i); and a curve 264represents the frequency response of a Sagnac sensor array 200 having 15km of fiber in its delay loop 214 and 300 meters of fiber in eachhydrophone 212(i). As illustrated, each Sagnac sensor array 200 has thesame sensitivity at low frequencies, but approaches a maximumsensitivity at different frequencies given by their respective loopproper frequencies. Thus, for a given minimum detectable acoustic signalat low frequencies, there is still some freedom in choosing the fiberlengths of the delay loop 214 and the hydrophones 212(i). This freedommay be used to help the Sagnac sensor array 200 satisfy other criteria,such as minimizing the total amount of fiber required or minimizing thedelay loop length.

Increasing the Dynamic Range of the Sagnac Sensor Array

As discussed above, the Sagnac sensor array 200 has a larger dynamicrange at low acoustic frequencies than the Mach-Zehnder interferometricsensor array because it is immune to phase noise. Ideally, an array 200provides enough dynamic range to detect the strongest and weakestacoustic signal which are likely to be encountered. This requirementoften translates into a required dynamic range of approximately 150 dB.In order to achieve such a large dynamic range in a Mach-Zehnderinterferometric sensor array, two separate sensors with different phaseresponsivities are required, with each detecting a fraction of the total150 dB dynamic range. The obvious disadvantage to this scheme is that itrequires two sensor arrays (i.e., twice as many hydrophones, rungs,sources and detectors). Effectively, an array which can support Nhydrophones can detect the acoustic signal at only N/2 points.

In the Sagnac sensor array 200, it is possible to achieve a largedynamic range without using additional hydrophones 212. Because thephase responsivity in the Sagnac sensor array is a function of thehydrophone responsivity and delay loop length, as shown in Equation 1,the phase responsivity of the entire array of hydrophones can be changedby modifying the delay loop length. By simultaneously using two separatedelay loops 214(1) and 214(2) of length L₁ and L₂, respectively, asshown in a modified sensor array 266 in FIG. 7, the detection range ofthe array 266 can be dramatically increased. The array 266 now has 2Nseparate Sagnac loops. Each hydrophone 212(i) returns a separate signalfor each of the two delay loop paths, and the length of each delay loop214(1), 214(2) determines the acoustic detection range of that signal.The total acoustic detection range of each hydrophone 212(i) is theunion of the detection ranges of each of the two Sagnac loop sensorswhich enclose the hydrophone 212(i). The lengths of L₁ and L₂ set theacoustic detection range. The length L₁+L₂ is chosen to allow the array266 to detect the smallest acoustic signal of interest. The length L₁ ofthe delay loop 214(1) is then chosen to place the detection range of thesignals which travel only this shorter delay loop on top of thedetection range of the signals which travel both delay loops 214(1),214(2). In a TDM system, as a result of the insertion of a second loop,the repetition frequency of the source pulses are halved in order toallow time for 2N pulses to return, and the lengths of the delay loops214(1), 214(2) are chosen such that there is no pulse overlap. Becausethe repetition frequency is halved, the dynamic range of each individualsignal decreases by 3 dB. This decrease is more than offset by theincrease in the total dynamic range achieved by piggybacking the dynamicrange of two separate signals. In FIG. 7, the second delay loop 214(2)is positioned such that all light passing through the second delay loop214(2) passes through the first delay loop 212(1). It should beunderstood that, alternatively, the two delay loops 214(1), 214(2) canbe optically in parallel such that the light which passes through thesecond delay loop 214(2) does not pass through the first delay loop214(1). In such case, the fiber length of the second delay loop 214(2)would have to be the sum of the first length and the second length(i.e., L₁+L₂). But, since L₁ is considerably shorter than L₂, thisadjustment is not essential. The embodiment of FIG. 7 reduces the totalfiber requirements by adding the length of the first delay loop to thesecond delay loop.

FIG. 8 illustrates the extended dynamic range made possible by using thetwo delay loops 214(1), 214(2) in the array 266 in which the dynamicrange of each signal is 100 dB and the ratio L1/L2 was set to be 5000.As shown, the array 266 is now able to detect over the entire dynamicrange of interest (approximately a 160-dB range) without increasing thehydrophone count.

Distributed Sensing

In the Sagnac sensor array 266, any phase modulation in theinterferometer can be transferred into an intensity modulation at theinterfering 3×3 coupler 220. This distributed sensing over the entireSagnac loop is disadvantageous for an acoustic sensor array. In order tobe practical, the acoustic sensor array should sample the acousticsignal at a number of discrete points in space (i.e., at thehydrophones) and return these signals independently. Mach-Zehnderinterferometric sensor arrays achieve this because the interferometer isconfined within a small space and thus only senses at that point. Inorder for the Sagnac sensor array 266 to be practical, the distributedsensing of the Sagnac loop must be decreased.

The bulk of the fiber in the interferometer constitutes the delay loop214, which can be located in two positions. The first is with the source222 and the detection electronics (i.e., the detector 230 and thedetector 232) in the dry end (i.e., out of the water), as shown in FIG.9A. Here the delay loop 214 can be environmentally shielded to minimizeany external modulation. However, downlead fibers 270, 272 which connectthe wet end to the array portion 210 are part of the interferometer. Thesecond possibility is to locate the delay loop 214 in the wet end (i.e.,in the water) with the array 210, as shown in FIG. 9B. As such, thedelay loop 214 cannot be isolated to the same extent as it could if itwere located in the dry end, but the downlead fibers 270, 272, 274 areoutside of the interferometer and thus are non-sensing. The relativemagnitude of the downlead and delay loop distributed pick-up dictateswhich configuration is best suited for a particular application. Itshould be noted that if the delay loop 214 is located in the dry end(FIG. 9A), the downlead fibers 270, 272 must remain stationary toprevent physical movements, such as bending and vibrations, of thesefibers, which can induce extremely large phase modulations. These arefiber motion induced phase modulations as opposed toacoustically-induced phase modulations. (Such physical movements areproblems in towed arrays, but may not be significant problems instationary arrays.) Thus, if the delay loop 214 is located in the dryend (FIG. 9A), the entire wet end of the Sagnac sensor array 210 must bestationary. However, with the delay loop 214 located in the wet end(FIG. 9B), only the portion to the right of the 3×3 coupler 220 in FIG.9B must remain stationary since the downlead fibers 270, 272, 274 arenot then part of the interferometer. When the delay loop 214 is locatedin the wet end (FIG. 9B), the delay loop fiber must be desensitized. Thedelay loop 214 can be made stationary by wrapping the delay loop fibersaround a desensitized cylinder (not shown), thereby eliminating fibermotion and making acoustic pick-up the dominant source of distributedpick-up signal. Because it is easier to desensitize fiber toacoustically-induced phase modulation than it is to desensitize fiber tomovement-induced phase modulation, the configuration which locates thedelay loop 214 in the wet end (FIG. 9B) is preferable for towed arrayapplications and will be described in more detail below.

Calculation of the Acoustic Pick-up Noise Induced in the Delay Loop

In this section, estimates are derived for the magnitude of theacoustically induced distributed pick-up noise as compared to theacoustically induced hydrophone phase modulation in the Sagnac sensorarray 210 of FIG. 9(b). The intensity modulation due to the distributedphase modulations resulting from the pick-up of acoustic signals in thedelay loop and bus fiber (the fiber connecting each hydrophone to thedelay loop and the 3×3 coupler) can be considered a source of noise. Forthe following discussion, consider one loop of the Sagnac sensor arrayas comprising only delay fiber of length L_(d), a bus fiber of lengthL_(b), a hydrophone fiber of length L_(h), and a total length L, asshown in FIG. 10. Also assume that L_(d) is much larger than L_(b) andL_(h). The phase responsivity of fiber to acoustic signals results froma pressure dependent propagation constant, β. In general, the pressuredependent component of the propagation constant at a position l and timet can be written as:

β(l,t)=β₀ R(l)P(l,t)  (2)

where β₀ is the zero-pressure propagation constant, R(l) is thenormalized phase responsivity of the fiber, and P(l,t) is the pressureas a function of space and time. If a sinusoidal acoustic signal offrequency Ω is assumed, Equation 2 can be rewritten as:

β(l,t)=β₀ R(l)[P ₀ +P _(m) sin(Ωt+θ(l))]  (3)

where P₀ is the steady-state pressure, P_(m) is the amplitude of thepressure modulation (assumed to be independent of l), and θ(l) containsthe spatial phase variation of the acoustic wave. In general, theinduced phase difference between interfering beams in a Sagnac loop dueto acoustically induced phase modulation from l=l₁ to l=l₂ is given bythe integral: $\begin{matrix}{{\varphi_{int}(t)} = {\int_{l_{1}}^{l_{2}}{\left\lbrack {{\beta \left( {l,{t + \frac{\left( {l - L} \right)}{v}}} \right)} - {\beta \left( {l,{t - \frac{l}{v}}} \right)}} \right\rbrack {l}}}} & (4)\end{matrix}$

where v is the speed of light in the fiber, and L is the loop length.Substituting Equation 3 into Equation 4 yields: $\begin{matrix}{{\varphi_{int}(t)} = {\beta_{0}P_{m}{\int_{l_{1}}^{l_{2}}{{{R(l)}\left\lbrack {{\sin \left( {{\Omega \left( {t + \frac{l - L}{v}} \right)} + {\theta (l)}} \right)} - {\sin \left( {{\Omega \left( {t - \frac{l}{v}} \right)} + {\theta (l)}} \right)}} \right\rbrack}{l}}}}} & (5)\end{matrix}$

Equation 5 can be used to determine the phase difference betweeninterfering beams due to acoustic modulation of the hydrophone, bus, anddelay fibers.

For the hydrophone fiber, Equation 5 is integrated from l₁=l_(d)+l_(b)/2to l₂=l_(d)+l_(b)/2+l_(h). It is assumed that θ(1) is constant over thisrange (i.e., that the acoustic wavelength is much larger than thedimension of the hydrophone). It is also assumed that the normalizedphase responsivity of the fiber, R(l), is constant and is equal to R_(h)in this range. Equation 5 then gives a phase difference amplitudebetween interfering beams due to hydrophone fiber modulation:$\begin{matrix}{{\varphi_{int}^{h} = {2\beta_{0}R_{h}P_{m}L_{h}{\sin \left( \frac{\Omega \cdot T_{delay}}{2} \right)}}},} & (6)\end{matrix}$

where it is assumed that ΩL_(h)/2v<<1. Note that Equation 2 agrees withthe expression given in Equation 1.

For the bus fiber, Equation 5 is integrated first from l₁=l_(d) tol₂=l_(d)+l_(b)/2, and then from l₁=L-l_(b)/2 to 1₂=L to include both theupper and lower bus lines. Again, it is assumed that R(l) is constantand equal to R_(b) for all bus fiber, such that θ(l) is constant in theintegral of Equation 5. The phase difference amplitude betweeninterfering beams due to fiber modulation becomes: $\begin{matrix}{{\varphi_{int}^{b} = {2\beta_{0}R_{b}P_{m}L_{b}{\sin \left( \frac{\Omega \cdot T_{delay}}{2} \right)}}},} & (7)\end{matrix}$

where it is assumed that ΩL_(h)/2v<<1. It should be emphasized that theassumptions on the constancy of θ(l) and the amplitude of ΩL_(h)/2v actto increase φ_(int) ^(b), thus giving a worst case scenario for the busfiber.

For the delay fiber, Equation 5 is integrated from l₁=0 to l₂=l^(d),and, as before, it is assumed that θ(l) is constant over this range(i.e., the delay loop coil is much smaller than the acousticwavelength), and that R(l) is constant and equal to R_(d) over theintegral. Equation 5 then yields a phase difference amplitude betweeninterfering beams due to delay fiber modulation given by:$\begin{matrix}{{\varphi_{int}^{d} = {{2\beta_{0}R_{d}{P_{m}\left( {L - L_{d}} \right)}{\sin \left( \frac{\Omega \quad T_{delay}}{2} \right)}} = {2\beta_{0}{R_{d}\left( {L_{b} + L_{h}} \right)}\sin \quad \left( \frac{\Omega \quad T_{delay}}{2} \right)}}},} & (8)\end{matrix}$

where it is assumed that Ω(L_(b)+L_(h))/2v<<1.

With Equations 6-8, the relative magnitude of these phase modulationsamplitudes can be computed. First, it is noted that a standard plasticcoated fiber has a normalized phase responsivity, R, of −328 dB re1/μPa, as described, for example, in J. A. Bucaro, et al., Optical fibresensor coatings, Optical Fiber Sensors, Proceedings of the NATO AdvancedStudy Institute, 1986, pp. 321-338. On the other hand, as described, forexample, in C. C. Wang, et al., Very high responsivity fiber optichydrophones for commercial applications, Proceedings of the SPIE-TheInternational Society for Optical Engineering, Vol. 2360, 1994, pp.360-363, a fiber wrapped around current hydrophones made from air-backedmandrels has a normalized phase sensitivity of −298 dB re 1/μPa, anincrease of 30 dB over standard fiber. If we assume that the delay loopand the bus fiber have the normalized phase responsivity of standardplastic coated fiber, and that the hydrophone fiber is wrapped around anair-backed mandrel, then the ratio of R_(h) to R_(b) or R_(d) isapproximately 30 dB. Therefore, under the simplifying assumption made toreach Equations 6-8, it can be found that: $\begin{matrix}{{\frac{\varphi_{int}^{h}}{\varphi_{int}^{d}} \approx \left( \frac{31}{1 + \left( {L_{b}/L_{h}} \right)} \right)}\text{and}} & (9) \\{\frac{\varphi_{int}^{h}}{\varphi_{int}^{b}} \approx {\left( \frac{31}{L_{b}/L_{h}} \right).}} & (10)\end{matrix}$

and

The ratio L_(b)/L_(h) is a function of the hydrophone position. For thefirst hydrophone, L_(b)/L_(h)≈0 making φ_(int) ^(h)/φ_(int) ^(d)=31 andφ_(int) ^(h)/φ_(int) ^(b) extremely large. For the last hydrophone,typical values of 100 meters and 1 km for L_(h) and L_(b), respectively,are used to arrive at φ_(int) ^(h)/φ_(int) ^(d)≈φ_(int) ^(h)/φ_(int)^(b)≈3. Thus, despite the fact that the hydrophone fiber constitutes arelatively small amount of the overall Sagnac loop, the magnitude of theacoustically induced phase modulations in the hydrophone fiber aregreater than the acoustically induced phase modulations in the delayloop fiber and in the bus fiber for even the furthest hydrophone. Thefollowing section describes a means for dealing with this level ofdistributed pick-up noise using empty rungs.

In order to evaluate the integral in Equation 5 for the delay loopfiber, it is assumed that R(l)=R_(d) for all l less than L_(d). It wasthis constancy of R(l) which eliminated any contribution to the integralof Equation 5 from l=(L−L_(d)) to L_(d) (because the integrand became anodd function about L/2). However, coiling a long length of fiber willresult in some dependence in R(l) on l (possibly because the inner layerof fiber has a different R than the outer layer). These variations inR(l) increase the delay loop pick-up from l=L−L_(d) to L_(d). In orderto reduce this pick-up, it is first noted that R(l) need only be an evenfunction around L/2 to make the integrand of Equation 5 an odd functionabout L/2. R(l) can be forced to be more symmetric about L/2 by wrappingthe delay loop in such a way as to position symmetric points of thefiber loop next to each other as shown in FIG. 11. Such a wrappingensures that symmetric points of the delay loop are positioned inproximity to each other so that any variations in R(l) due to theposition of the fiber on the coil are as symmetric about L/2 aspossible, thereby making the delay loop pick-up as close to theexpression of Equation 8 as possible. Note that, because each Sagnacloop in the Sagnac sensor array has a different L/2 point, only one loopcan be wrapped exactly as shown in FIG. 11, thereby introducing a smalldegree of oddness in R(l) to all but one of the Sagnac loops.

It should also be mentioned that in addition to enhancing the acousticsensitivity of fiber with a hydrophone, it is possible to desensitizefibers by applying a metallic coating of a particular diameter. (See,for example, J. A. Bucaro, Optical fibre sensor coatings, cited above.)Measured normalized phase responsivities as low as −366 dB re 1/ μPahave been reported. If such fibers are used in the delay or bus lines,the ratio of R_(h) to R_(b) or the ratio of R_(h) to R_(d) approaches 68dB (instead of 30 dB with plastic coated delay and bus fibers),increasing the hydrophone induced signal over the delay and bus inducedsignal by 38 dB.

Reducing the Distributed Pick-up Noise by Using Empty Rungs

In order to further eliminate distributed pick-up signal, thehydrophone-induced acoustic modulation can be isolated from thedistributed pick-up modulation by placing empty rungs 300 that do notcontain a hydrophone in the array 210, as shown in FIG. 12. Each rung216(i) which contains a hydrophone 212(i), called a sensing rung, isproceeded by one of the empty rungs 300(i). The fact that thenon-sensing fiber of each loop which encloses an empty rung 300(i) isnearly identical to the non-sensing fiber of the loop which encloses thecorresponding sensing rung 212(i) means the empty rung 300(i) and thecorresponding sensing rung 212(i) will have nearly the same distributedpick-up signal. By treating this empty rung 300(i) as another sensor inthe array 210 and properly timing the pulses (in the TDM scheme) fromthe empty rungs 300(i) and the sensing rungs 212(i) so that they do notoverlap, the distributed pick-up signal present on each sensing rung212(i) can be measured. After detection, this signal can be subtractedfrom the sensing rung signal, leaving only intensity variations producedby phase modulations in the hydrophone fiber. Implementing such a schemerequires 2N rungs for an N sensor array 210, thereby reducing the dutycycle of individual signals by one half.

If desensitizing the bus portion of the array 210 is not required, asingle empty rung 300 can be placed in the array 210 to measure thedistributed pick-up signal associated with the delay loop 214, therebyrequiring only N+1 rungs (N sensing rungs 212(i) and one empty rung 300)for N sensors. If one empty rung 300 does not adequately measure thedistributed pick-up signal for each sensing rung 212(i), more emptyrungs 300 can be added at periodic intervals along the array, until thedistributed pick-up signal present on each sensing rung 212(i) can beadequately measured by the nearest of these empty rungs 300. Using fewerempty rungs results in a higher duty cycle for individual signals. FIG.12 depicts the extreme in which an empty rung was added for everysensing rung.

Polarization

For maximum contrast in any interferometric sensor, the state ofpolarization (SOP) of the interfering beams must be identical when theyrecombine. If they are orthogonal, there is no interference and thus noamplitude-modulated signal. This is referred to as polarization-inducedsignal fading. Because each sensor in the Sagnac sensor array is aSagnac loop, the research carried out so far on polarization-inducedsignal fading in the Sagnac fiber gyroscope applies to the Sagnac sensorarray as well. One promising solution is to place a depolarizer withinthe Sagnac loop. (See, for example, K. Böhm, et al., LOW-DRIFT FIBREGYRO USING A SUPERLUMINESCENT DIODE, ELECTRONICS LETTERS, Vol. 17, No.10, May 14, 1981, pp. 352-353.) The depolarizer ensures that at leasthalf of the optical power is returning to the 3×3 coupler in the correctSOP at all times. This general approach produces a constant visibilityregardless of the loop birefringence. (See, for example, William K.Burns, et al., Fiber-Optic Gyroscopes with Depolarized Light, JOURNAL OFLIGHTWAVE TECHNOLOGY, Vol. 10, No. 7, July 1992, pp. 992-999). Thesimplest configuration uses an unpolarized source such as a fibersuperfluorescence source and a depolarizer in the loop. As illustratedin FIG. 13, in the Sagnac sensor array 200, one depolarizer 310 isplaced at a point which is common to all the Sagnac loops. Thedepolarizer 310 ensures that each sensor 212(i) has this constantvisibility independent of birefringence as long as the loopbirefringence remains constant. This represents a great simplificationin the handling of polarization-induced signal fading over those methodsused in Mach-Zehnder interferometric sensor arrays.

Although slow changes in the birefringence will be sufficiently canceledby the reciprocal nature of the Sagnac interferometer, birefringencemodulations at frequencies in the acoustic range of interest willproduce polarization noise. Most birefringence modulation at thesefrequencies occurs as a result of physical fiber movement. Thus, theSagnac loop should remain stationary in order to reduce the polarizationnoise (as well as the distributed pick-up signal).

Noise Sources Introduced by the use of the Sagnac Interferometer

Thermal Phase Noise

Because the index of refraction of the fiber changes with temperature,thermal fluctuations in a fiber will produce phase fluctuations in thelight traveling through it. These index variations are uncorrelated overthe length of fiber, and thus the resulting phase fluctuations scale asthe square root of length. Because Mach-Zehnder interferometerstypically use less than 100 meters of fiber in each arm, the magnitudeof this thermal phase noise is negligible. The Sagnac interferometer hasa great deal more fiber in the interferometer and as a result, thermalphase noise can become a limiting noise source. The magnitude of thisthermal phase noise in a Sagnac interferometer has been describedtheoretically and confirmed by experiment. (See, for example, SverreKnudsen, et al., Measurements of Fundamental Thermal Induced PhaseFluctuations in the Fiber of a Sagnac Interferometer, IEEE PhotonicsTechnology Letters, Vol. 7, No. 1, 1995, pp. 90-93; and Kjell Kråkenes,et al., Comparison of Fiber-Optic Sagnac and Mach-ZehnderInterferometers with Respect to Thermal Processes in Fiber, JOURNAL OFLIGHTWAVE TECHNOLOGY, Vol. 13, No. 4, April 1995, pp. 682-686.). Forloops greater than 2 km, the thermal phase noise can exceed 1μrad/{square root over (Hz)} in the frequency range of interest, whichis on the order of the required array sensitivity.

The thermal phase noise can be considered as a source of distributedpick-up noise, akin to an external modulation to the delay loop, and assuch can be reduced by using empty rungs, as described above. Thermalphase noise can also be reduced by shortening the loop length. Asdiscussed above, the loop length can be shortened without changing thelow frequency sensitivity by increasing the hydrophone fiber length bythe same factor as that by which the delay loop was decreased. Forexample a 40-km delay loop with 50 meters of hydrophone fiber has thesame low-frequency response as a 20-km delay loop with 100 meters offiber. The latter combination however will suffer less thermal phasenoise because the total delay loop length is shorter by almost a factorof two.

Kerr Effect Induced Phase Noise

Kerr-induced phase shifts which can be generated in a Sagnacinterferometer have received a great deal of attention for the fiberoptic gyroscope. (See, for example, R. A. Bergh, et al., Sourcestatistics and the Kerr effect in fiber-optic gyroscopes, OPTICSLETTERS, Vol. 7, No. 11, November 1982, pp. 563-565; R. A. Bergh, etal., Compensation of the optical Kerr effect in fiber-optic gyroscopes,OPTICS LETTERS, Vol. 7, No. 6, June 1982, pp. 282-284; and N. J. Frigo,et al., Optical Kerr effect in fiber gyroscopes: effects ofnonmonochromatic sources, OPTICS LETTERS, Vol. 8, No. 2, February 1983,pp. 119-121.) The demands of the gyroscope and the acoustic sensor,however, are different because the gyroscope measures DC levels. SmallDC offsets created by Kerr-induced phase shifts which would limit afiber gyroscope are non-issues with an acoustic sensor. The Kerr-inducedDC phase shift is not a problem as long as it does not move the biaspoint too far away from quadrature. The intensity noise on the lightsource can produce a Kerr induced phase noise on the output. However,the magnitude of this Kerr-induced AC phase noise is small as long asthe Kerr-induced DC phase shift remains small. The origin ofKerr-induced phase shifts in the Sagnac sensor array is different thanin the fiber gyroscope. The asymmetry of the Sagnac sensor array invitessuch a Kerr phase shift much more readily than the nominally symmetricgyroscope does. That asymmetry results from the array portion as well asany placement of EDFAs which are asymmetric, in that one beam sees gainbefore propagating through the delay loop, then sees loss, while thecounter-propagating beam sees loss, then sees gain. It is possible tobalance these asymmetries and null the Kerr-induced phase shift bychoosing the proper location for EDFAs in the delay loop. The specificsdepend on the exact array configuration and which multiplexing scheme isused.

Non-linear Phase Modulation Resulting from the EDFAs

The population inversions created in the EDFAs induce a phase shift onthe signal light that passes through it. (See, for example, M. J. F.Digonnet, et al., Resonantly Enhanced Nonlinearity in Doped Fibers forLow-Power All-Optical Switching: A Review, OPTICAL FIBER TECHNOLOGY,Vol. 3, No. 1, January 1997, pp. 44-64.) This phenomenon has been usedto produce all-optical interferometric switches. In a Sagnac sensorarray, the EDFAs within the interferometer create a nonlinear phaseshift via the same mechanism. Variations in the population inversion dueto pump or signal power fluctuations will produce phase modulationswhich will be converted to an intensity noise.

In order to estimate the magnitude of this noise source, a determinationmust be first made as to how the inverted population responds to pumpand signal power fluctuations. This is relatively straightforward to doby invoking the rate equations for an erbium system:

N ₁ +N ₂ =N ₀,  (11)

$\begin{matrix}{{{\frac{}{t}N_{2}} = {{\frac{I_{p}\sigma_{p}^{a}}{h\quad v_{p}A_{eff}}N_{1}} + {\frac{I_{s}\sigma_{s}^{a}}{{hv}_{s}A_{eff}}N_{1}} - {\frac{I_{p}\sigma_{p}^{e}}{{hv}_{p}A_{eff}}N_{2}} - {\frac{I_{s}\sigma_{s}^{e}}{{hv}_{s}A_{eff}}N_{2}} - \frac{N_{2}}{\tau_{2}}}},} & (12)\end{matrix}$

where N₁ and N₂ are the population densities of the lower and excitedstates respectively, N₀ is the total population density, I is theintensity, σ is the cross section, A_(eff) is the effective mode area inthe fiber, and τ₂ is the lifetime of level two. The subscripts p and sdenote pump and signal, respectively, and the superscripts a and edenote absorption and emission, respectively.

By splitting N₁, N₂, I_(p), and I_(s) into their steady-state andtime-varying components, then substituting this into Equation 12 andcombining Equation 12 with Equation 11, the result is: $\begin{matrix}{{\frac{}{t}{N_{2}(t)}} = {{\left\lbrack \frac{{N_{o}\sigma_{p}^{a}} + {N_{2}^{ss}\left( {\sigma_{p}^{e} + \sigma_{p}^{a}} \right)}}{{hv}_{p}} \right\rbrack {I_{p}(t)}} + {\left\lbrack \frac{{N_{o}\sigma_{s}^{a}} + {N_{2}^{ss}\left( {\sigma_{s}^{e} + \sigma_{s}^{a}} \right)}}{{hv}_{s}} \right\rbrack {I_{s}(t)}} + {\left\lbrack {\frac{I_{p}^{ss}\left( {\sigma_{p}^{e} + \sigma_{p}^{a}} \right)}{{hv}_{p}} + \frac{I_{s}^{ss}\left( {\sigma_{s}^{e} + \sigma_{s}^{a}} \right)}{{hv}_{s}} + \frac{1}{\tau_{2}}} \right\rbrack {N_{2}(t)}} - {\left\lbrack {{{\left. \frac{\left( {\sigma_{p}^{e} + \sigma_{p}^{a}} \right)}{h\quad v_{p}} \right\rbrack{I_{p}(t)}{N_{2}(t)}} - {\left\lbrack \frac{\left( {\sigma_{s}^{e} + \sigma_{s}^{a}} \right)}{{hv}_{s}} \right\rbrack {I_{s}(t)}{N_{2}(t)}}},} \right.}}} & (13)\end{matrix}$

where the superscript ss denotes steady-state values, and thetime-varying components are now written as explicit functions of time(N₂=N₂ ^(ss)+N₂(t)). If it is assumed that N₂(t) is much smaller than N₂^(ss), then the last two terms in Equation 13 can be neglected. Bywriting I_(p)(t)=I_(p) ^(m) sin(f_(p)t) and I_(s)(t)=I_(s) ^(m)sin(f_(s)t) (where I_(p) ^(m) and I_(s) ^(m) denote the modulationamplitudes of I_(p)(t) and I_(s)(t), respectively, and f_(p) and f_(s)respectively denote the pump and signal modulation frequencies) andsolving the resulting differential equations, it can be found that:$\begin{matrix}{\frac{{N_{2}\left( f_{p} \right)}}{N_{2}^{ss}} \approx {\left( {\frac{{\sigma_{s}^{a}\left( {\sigma_{p}^{e} + \sigma_{p}^{a}} \right)} - {\sigma_{p}^{a}\left( {\sigma_{s}^{e} + \sigma_{s}^{a}} \right)}}{\left( {\sigma_{p}^{e} + \sigma_{p}^{a}} \right)^{2}} \cdot \frac{v_{p}}{v_{s}}} \right){\left( \frac{1}{\sqrt{1 + {f_{p}^{2}/f_{o}^{2}}}} \right) \cdot \frac{I_{s}^{ss}I_{p}^{m}}{I_{p}^{{ss}^{2}}}}}} & (14)\end{matrix}$

$\begin{matrix}{\frac{{N_{2}\left( f_{s} \right)}}{N_{2}^{ss}} \approx {\left( {\frac{\sigma_{s}^{a}}{\sigma_{p}^{a}} - \frac{\sigma_{s}^{e} + \sigma_{s}^{a}}{\sigma_{p}^{e} + \sigma_{p}^{a}}} \right) \cdot \frac{v_{p}}{v_{s}} \cdot \left( \frac{1}{\sqrt{1 + {f_{p}^{2}/f_{o}^{2}}}} \right) \cdot \frac{I_{s}^{m}}{I_{p}^{ss}}}} & (15)\end{matrix}$

where: $\begin{matrix}{{f_{o} = {{{\frac{\sigma_{p}^{e} + \sigma_{p}^{a}}{{hv}_{p}}I_{p}^{ss}} + {\frac{\sigma_{s}^{e} + \sigma_{s}^{a}}{{hv}_{s}}I_{s}^{ss}} + \frac{1}{\tau_{2}}} \approx {\frac{\sigma_{p}^{e} + \sigma_{p}^{a}}{{hv}_{p}}I_{p}^{ss}}}}\quad {{{when}\quad I_{p}^{ss}}\operatorname{>>}{I_{s}^{ss}.}}} & (16)\end{matrix}$

If it is assumed that λ_(p)=1480 nm, λ_(s)=1550 nm, and I_(p) ^(ss)=1 W,and if typical erbium-silica cross sections are assumed, then Equations14 and 15 simplify to: $\begin{matrix}{\frac{{N_{2}\left( f_{p} \right)}}{N_{2}^{ss}} \approx {\left( \frac{0.9}{\sqrt{1 + {{f_{p}^{2}/4.3}\quad {kHz}}}} \right) \cdot \frac{I_{s}^{ss}I_{p}^{m}}{I_{p}^{{ss}^{2}}}}} & (17) \\{\frac{{N_{2}\left( f_{s} \right)}}{N_{2}^{ss}} \approx {\left( \frac{1.2}{\sqrt{1 + {{f_{s}^{2}/4.3}\quad {kHz}}}} \right) \cdot {\frac{I_{s}^{m}}{I_{p}^{ss}}.}}} & (18)\end{matrix}$

The pump-induced population inversion fluctuations (Equation 17) will beanalyzed first. If I_(s) ^(ss)=1 mW, I_(p) ^(ss)=1 W, and it is assumedthat I_(p) ^(m)/I_(p) ^(ss)=10⁻⁶/{square root over (Hz)} (120 dB/{squareroot over (Hz)} electronic SNR), then |N₂(f_(p))|/N₂ ^(ss)=9×10⁻¹⁰{square root over (Hz)}⁻¹ at frequencies well below 4.3 kHz. In order toconvert this figure to a phase modulation, the fact that 10 mW of pumppower absorbed in an erbium-doped fiber induces approximately 7 radiansof phase shift at 1550 nm can be used. (See, for example, M. J. F.Digonnet, et al., Resonantly Enhanced Nonlinearity in Doped Fibers forLow-Power All-Optical Switching: A Review, OPTICAL FIBER TECHNOLOGY,Vol. 3, No. 1, January 1997, pp. 44-64.) Using simulations, 10 mW ofabsorbed pump power in a typical erbium-doped fiber providesapproximately 6 dB of small signal gain at 1550 nm, which is close tothe gain required by each amplifier in an array with distributed EDFAs.(See, for example, Craig W. Hodgson, et al., Optimization of Large-ScaleFiber Sensor Arrays Incorporating Multiple Optical Amplifiers-Part I:Signal-to-Noise Ratio; Craig W. Hodgson, et al., Optimization ofLarge-Scale Fiber Sensor Arrays Incorporating Multiple OpticalAmplifiers-Part II: Pump Power; Jefferson L. Wagener; et al., NovelFiber Sensor Arrays Using Erbium-Doped Fiber Amplifiers; and C. W.Hodgson, et al., Large-scale interferometric fiber sensor arrays withmultiple optical amplifiers, cited above.) Therefore, each amplifierprovides approximately 7 radians of DC phase shift. Since the nonlinearphase shift is proportional to the upper state population, N₂, it can bewritten that ΔN₂/N₂ ^(ss)=Δφ/φ^(ss). Using this relation and Equation 17again for I_(s) ^(ss)=1 mW, I_(p) ^(ss)=1 W, I_(p) ^(m)/I_(p)^(ss)=10⁻⁶/{square root over (Hz)} and f_(s)<<4.3 kHz, the low-frequencyphase noise induced by each EDFA is (7 radians)×(9×10⁻¹⁰){square rootover (Hz)}⁻¹=6.3×10⁻⁹ rad/{square root over (Hz)}. If it is assumed thatthere are a total of 500 such amplifiers and that the phase modulationsfrom all 500 amplifiers add coherently, the total pump noise inducedphase shift can be estimated to be 3.2 μrad/{square root over (Hz)}. Thetarget phase noise floor is typically set to 1 μrad/{square root over(Hz)}, indicating that the nonlinear phase-noise induced by the EDFAsdue to pump power fluctuations is close to but not significantly largerthan the required phase noise floor. In practice, the amplifiers' phasemodulations will not add coherently, which will reduce the 3.2μrad/{square root over (Hz)} figure.

Calculations of the induced phase shift due to signal power fluctuationsare more complicated because the signal power not only has intensitynoise but is also modulated by the multiplexing scheme. Againconsidering the TDM case, in general, while a given pulse is travelingthrough a particular EDFA, there may or may not be a counter-propagatingpulse traveling through that EDFA at the same time. Taking the worstcase in which there is always a counter-propagating pulse, I_(s) ^(m) istwice the intensity noise of each individual pulse. For the amplifiers,I_(s) ^(m) is typically 1.5 to 2 times the intensity noise of eachindividual pulse. Assuming the signal light has an electronic SNR of 120dB/{square root over (Hz)} at acoustic frequencies (i.e., I_(s)^(m)/I_(s) ^(ss)=10⁻⁶/{square root over (Hz)}), and inserting thisfigure into Equation 18 along with I_(p) ^(ss)=1 W and I_(s) ^(m)=2 mW,it can be calculated that |N₂(f_(s))|/N₂ ^(ss) is approximately 2.4×10⁻⁹{square root over (Hz)}⁻¹ at frequencies much lower than 4.3 kHz andthat the phase noise induced by signal intensity noise in each EDFA isthus 1.68×10⁻⁸ rad/{square root over (Hz)}. Again assuming 500amplifiers and coherent addition of all EDFA-induced phase modulation,the total EDFA induced phase noise on each pulse is 8.4 μrad/{squareroot over (Hz)}, a level which could again limit the performance of theSagnac sensor array. However, a more detailed study taking into accountthe multiplexing scheme and exact timing of the array is needed for amore accurate calculation.

Multiplexing Schemes in a Sagnac Array

Time-Division Multiplexing

It has been assumed thus far that the Sagnac sensor array is operated ina TDM configuration. It should be noted that, in the Sagnac sensorarray, the source requirements for such a TDM system are not asdemanding as those of a Mach-Zehnder interferometric sensor array in aTDM configuration. The reason for this is the use of the broadbandsource in the Sagnac sensor array. In the Mach-Zehnder interferometricsensor array, the light from adjacent rungs is coherent due to thenarrow linewidth source, and thus extremely high extinction ratios onthe input pulse are required to prevent multi-path coherentinterference. These high extinction ratio requirements are achieved byplacing multiple modulators in series, which results in a complicated,high loss, and expensive source. In the Sagnac sensor array, therequired extinction ratio need not be as high because the broadbandsource eliminates any possibility of multi-path coherent interference.In addition, the narrow linewidths required by the Mach-Zehnderinterferometric sensor array prevent the use of a pulsed laser source inplace of a continuous wave (cw) laser source which is externallymodulated with Lithium Niobate intensity modulators. In the Sagnacsensor array, either a continuous-wave ASE source which is externallymodulated, a pulsed ASE source, or some combination thereof could beused to construct the source. Again, the reason for this is that theSagnac sensor array does not require a narrow linewidth source. Althoughthe present invention does not require a narrow linewidth source, itshould be understood that the Sagnac sensor array of the presentinvention can be used with a narrow linewidth source, such as, forexample, a laser.

Frequency Division Multiplexing

The use of the broadband source also allows the Sagnac sensor array tooperate in non-TDM configurations without changing the design orrequiring additional sources. Frequency division multiplexing (FDM) iscommonly used with Mach-Zehnder interferometric sensor arrays using thePhase-Generated Carrier (PGC) scheme but is also compatible with theSagnac sensor array. FIG. 14 shows a basic Sagnac sensor array 400 usinga FDM scheme. A fiber superfluorescent source (SFS) 402 (or otherbroadband source, such as, for example, an LED) generates input light. Achirped intensity modulation is applied to the input light via anintensity modulator 404 which is controlled by a chirped frequencygenerator 406. The modulated light enters a sensor array 410 via a 3×3coupler 412. The light passes through a delay loop 414 and pluralsensing rungs 416(i) having respective sensors 418(i). Empty rungs (notshown) can also be included if desired. After passing through the delayloop 414 and the rungs 416(i), the light exits from the sensor array 410through the coupler 412 and is detected by a detector 420 whichgenerates an electrical output signal responsive to the detected light.The electrical output signal from the detector 420 is mixed in a mixer422 with the same chirped frequency which has been time delayed by adelay 424 which delays the chirped frequency by a time Δt. In the setupillustrated in FIG. 14, the output of the mixer 422 is applied to aspectrum analyzer 426. In an operational embodiment, the output of themixer 422 is applied to a signal processing subsystem (not shown) whichanalyzes the output of the mixer 422 to reproduce the acoustic signalsimpinging on the array 410.

The signals returning from the sensors 418(i) in the various rungs416(i) are further delayed with respect to the delayed chirp frequency.This is illustrated by the graphs in FIG. 15 by the original chirpedfrequency 450, the delayed chirped frequency 452 from the delay 424, thechirped return signal 460 from the first rung, the chirped return signal462 from the second rung and the chirped return signal 464 from thethird rung. In the mixer 422, separate beat frequencies f_(b) ₁ 470,f_(b) ₂ 472, f_(b) ₃ 474, respectively (shown in FIG. 14), are formedbetween the mixing chirped frequency 452 and each of the signalsreturning from the various rungs in the Sagnac sensor array 410. (See,for example, S. F. Collins, et al., A Multiplexing Scheme For OpticalFibre Interferometric Sensors Using An FMCW Generated Carrier, OFS '92Conference Proceedings, pp. 209-211.) Although only three chirped returnsignals 460, 462, 464 are illustrated in FIG. 15, it is contemplatedthat up to N return signals can be provided, where N is the number ofrungs in the array 410. The chirped return signals from the Nth rungcauses a beat frequency f_(bN) in the mixer 422.

As illustrated by a pictorial representation of a spectral output inFIG. 14, acoustic modulation of the signals will appear as uppersidebands 480, 481, 482 and lower sidebands 484, 485, 486 to the beatfrequencies. An advantage of this FDM scheme is that the demands on thearray timing are greatly relaxed over those required in a TDM system. ATDM system requires a specific delay between adjacent rungs in order toprevent pulses from overlapping, and this can present a demandingengineering problem. In FDM, variations in fiber lengths shift beatfrequencies but do not induce overlap between signals as long as thesebeat frequencies are separated by twice the acoustic detection range.The latter is accomplished by selecting the proper chirp rate. Unlike ina TDM system, all paths return light at all times, which can result inphase noise between the different incoherent signals. The broadband ASElight source minimizes the magnitude of this phase noise. (See, forexample, Moslehi, Analysis of Optical Phase Noise in Fiber-Optic SystemsEmploying a Laser Source with Arbitrary Coherence Time, Journal ofLightwave Technology, Vol. LT-4, No. 9, September 1986, pp. 1334-1351.)

Code Division Multiplexing

Code division multiplexing (CDM) has received increased attention latelyfor its use in sensor arrays. (See, for example, A. D. Kersey, et al.,Code-division Multiplexed Interferometric Array With Phase NoiseReduction And Low Crosstalk, OFS '92 Conference Proceedings, pp.266-269; and H. S. Al-Raweshidy, et al., Spread spectrum technique forpassive multiplexing of interferometric optical fibre sensors, SPIE,Vol. 1314 Fibre Optics '90, pp. 342-347.) As illustrated for a Sagnacsensor array 600 in FIG. 16, in CDM, the input light from a fibersuperfluorescent source 602 (or other broadband source, such as, forexample, an LED) is modulated in an intensity modulator 604 according toa pseudo-random code generated by a code generator 606. The modulatedlight is applied to an interferometric loop 608 via a 3×3 coupler 610and propagates through a delay loop 614 and a plurality of rungs 616(i)in an array 612. In the illustrated embodiment, each rung 616(i)includes a respective sensor 618(i). Empty rungs (not shown) can also beincluded if desired. The light returns from the loop via the 3×3 coupler610 and is detected by a detector 620. The electrical output of thedetector 620 is applied to a correlator 622 along with the output of thecode generator 606, which output is delayed for a duration τ_(cor) by adelay 624. The bit duration of the pseudo-random code is shorter thanthe propagation delay between adjacent rungs in the array 612. Whenτ_(cor) is equal to one of the loop travel times τ_(i), through arespective rung 616(i), then the signal returning from this sensor inthe rung 616(i) is correlated to the delayed pseudo-random code. Theother signals, which have delays τ_(j) where |τ_(j)−τ_(i)|>τ_(bit)correlate to zero. The correlation process involves, for example,multiplying the detected signal by 1 or −1 (or gating the signal in anelectronic gate 630 to the non-inverting and inverting inputs of adifferential amplifier 632) depending on whether the correlating code ison or off. The output of the differential amplifier on a line 634 is thecorrelated output. The signal is then time averaged over a periodt_(avg) equal to the duration of the code. The uncorrelated signals timeaverage to zero, thereby isolating the signal from sensor 618(i).τ_(cor) is scanned to retrieve sequentially the signals from allsensors.

An advantage of CDM over TDM is that the delay between sensors does nothave to be controlled accurately. Any loop delays τ_(j) in which|τ_(j)−τ_(j)±11|>τ_(bit) is acceptable (where τ_(bit) is the duration ofa pulse in the code). Correlating requires a knowledge of the τ_(j)'s,which are easily measured. As with FDM, the use of a broadband sourcebenefits reducing the phase noise which results from the addition of allthe signals together.

The foregoing described a novel design for an acoustic sensor arraybased on the Sagnac interferometer. The major advantages of this designare the use of common-path interferometers. This eliminates theconversion of source phase noise into intensity noise, which isprevalent in Mach-Zehnder interferometric sensors, and allows the use ofa cheap, high-power ASE source or other broadband source. The responseof the Sagnac sensor array as a function of acoustic frequency is shownto match the ocean noise floor. The design also allows the dynamic rangeto be dramatically increased without adding hydrophones by using oneadditional, very short delay loop. A technique for eliminatingpolarization-induced signal fading was discussed above. The Sagnacsensor array also allows the use of several multiplexing schemes in asimpler form than is achievable with a standard Mach-Zehnder array.Because of these features, the Sagnac sensor array design provides avery promising alternative to Mach-Zehnder-interferometer-based sensorarrays.

Folded Sagnac Sensor Array

FIGS. 17-20 illustrate alternative embodiments of a distributed acousticsensor array based upon the Sagnac effect which has an architecturemodified to reduce the distributed pick-up from the downlead fibers. Inparticular, FIG. 17 illustrates a basic folded Sagnac acoustic fibersensor array 700 which comprises a source 702, a first detector 704 anda second detector 706. Preferably, the source 702, the first detector704 and the second detector 706 are located in the dry end of the sensorarray 700 (e.g., on shore or on board a ship).

The source 702 generates light pulses which are coupled to a 3×3 coupler710 via a downlead fiber 708. As illustrated, the 3×3 coupler is locatedin the wet end (e.g., proximate to the ocean floor). The 3×3 coupler 710has a first output port coupled to one end of a common fiber rung (rung0) 712, has a second output port coupled to a first array input/outputfiber 714 of an array 716, and has a third output port which isnon-reflectively terminated. Approximately 33 percent of the light fromthe source 702 is coupled to each of the first and second ports of the3×3 coupler and thus approximately 33 percent of the light propagates tothe common fiber rung 712 and approximately 33 percent of the lightpropagates to the array 716. As discussed above, although describedherein as a 3×3 coupler 710, other n×m couplers (e.g., a 2×2 coupler, a4×4 coupler, etc.) can be used with the embodiment of FIG. 17 and thealternative embodiments of the present invention described below.

The array 716 comprises a plurality of rungs 718(i) (i.e., 718(1),718(2) . . . 718(N)) coupled between the first array input/output fiber714 and a second array input/output fiber 720. Each rung 718(i) includesa respective acoustic sensor (i.e., hydrophone) 722(i). The array 716advantageously includes distributed erbium doped fiber amplifiers(EDFAs) 724, such as described above in connection with FIG. 3. (Thepump source for the EDFAs 724 is not shown in FIG. 17.) Althoughdescribed herein with respect to the array 716, other arrayconfigurations can also advantageously be used in the present invention.

The second array input/output fiber 720 couples the array 716 to a firstport of a 2×2 coupler 730. A second end of the common rung (rung 0) 712is coupled to a second port of the 2×2 coupler 730. Although describedherein as an array 716 comprising plural sensors 722(i), it should beunderstood that the present invention has applications for a sensorsystem having only a single sensor 722.

A third port of the 2×2 coupler 730 is nonreflectively terminated at aterminal 732. A fourth port of the 2×2 coupler 730 is coupled to a delayloop downlead fiber 740. The delay loop downlead fiber 740 couples thefourth port of the 2×2 coupler to a first end of a delay loop 750. Thedelay loop 750 may be located either in the dry end as shown or in thewet end. A second end of the delay loop 750 is coupled to a reflector752 such that light exiting the second end of the delay loop 750 isreflected back into the delay loop 750, propagates through the delayloop 750 and propagates through the delay loop downlead fiber 740 backto the fourth port of the 2×2 coupler 730. The light returned from theloop downlead fiber 740 is divided by the 2×2 coupler 730 withsubstantially equal portions propagating in the common rung 712 and inthe array 716 with both portions propagating toward the 3×3 coupler 710.The two portions are combined in the 3×3 coupler 710 where light pulseswhich have traveled the same distance through the array 716 and throughthe common rung 712 interfere and light pulses which have traveleddifferent distances do not interfere. The signals resulting from theinterference are output from the 3×3 coupler 710 as first and secondoutput signals which respectively propagate to the first detector 704via a first detector downlead fiber 770 and propagate to the seconddetector 706 via a second detector downlead fiber 772. The detectors704, 706 generate electrical output signals which are analyzed byelectronics (not shown) in a conventional manner to reproduce theacoustic signals impinging on the sensors 722(i). As discussed below,the signals which interfere within the 3×3 coupler 710 return from eachsensor 722(i) at different times, and can therefore be separated by timedivision multiplexing, frequency multiplexing, code divisionmultiplexing, or the like, as discussed above. The non-interferingsignals do not generate detectable output signals and are ignored.

The embodiment of FIG. 17 can be further modified by inserting adepolarizer (not shown) in one of the fiber segments 712, 714 or 720 inconjunction with an unpolarized source, as described above in connectionwith the Sagnac interferometer. Such embodiments will be described belowin connection with FIGS. 23A, 23B and 23C.

The light in a single pulse from the source 702 will now be tracedthrough the sensor array 700. A source pulse from the source 702 islaunched and travels down the source downlead 708 and through the 3×3coupler 710 to the common rung 712 and to the array 716. Together, thecommon rung 712 and the N rungs 718(i) in the array 716 provide N+1separate paths for the source pulses to travel to the 2×2 coupler 730.Because there are N+1 separate paths for the source pulse to travel, thesource pulse is split into N+1 separate pulses which pass through the2×2 coupler 730 and travel down the delay loop downlead 740 to the delayloop 750. After passing through the delay loop 750, the N+1 pulses arereflected by the reflector 752 and then propagate back through the delayloop 750, down the delay loop downlead 740 to the 2×2 coupler 730 in thewet end, still as N+1 separate pulses. Each of the N+1 pulses is againsplit into N+1 pulses in the common rung 712 and the N rungs 718(i).After passing back through the common rung 712 and the rungs 718(i), the(N+1)² pulses are combined in the 3×3 coupler 710 and then return downthe detector downleads 770, 772 back to the dry end where the pulses aredetected by the first and second detectors 704, 706 and analyzed.

Because there are (N+1)² possible separate combinations of paths fromthe source 702 to the reflector 752 and back to the detectors 704, 706,there are (N+1)² returned pulses. The only pulses that will interfere ina useable manner are pairs of pulses which travel the same exact pathlength but in opposite order. For the purposes of the followingdiscussion, a pulse will be identified by two numbers where the firstnumber identifies the path taken by the pulse from the source 702 to thereflector 752, and the second number identifies the path taken by thepulse from the reflector 752 back to the detectors 704, 706. Forexample, the pulse 0, 1 travels through the common rung (rung 0) 712,then through the delay loop 750, to the reflector 752, back through thedelay loop 750, and then through rung 718(1). The pulse 1,0 travelsfirst through the rung 718(1), then through the delay loop 750, to thereflector 752, back through the delay loop 750, and then through thecommon rung (rung 0) 712. Because the distance traveled by the pulse 0,1is identical with the distance traveled by the pulse 1,0, the pulse 0,1and the pulse 1,0 interfere when combined at the 3×3 coupler 710 andtherefore define a common-path interferometer (i.e., a folded Sagnacinterferometer) in the same manner as the Sagnac interferometersdescribed above. Acoustic sensing results from the hydrophone 722(1)which is placed in rung 1 which responds to acoustic modulation. Theinterfering pulses 0,1 and 1,0 see the hydrophone 722(1) at differenttimes and thus pick-up a phase difference due to the time varyingacoustic modulation of the hydrophone 722(1). At the 3×3 coupler 710,this phase difference is converted into an intensity modulation which istransmitted down the detector downleads 770, 772 to the detectors 704,706. The same effect occurs for the pulses 0,2 and 2,0, for the pulses0,3 and 3,0, etc.

Because the folded Sagnac interferometer is common-path, the source 702can have a short coherence length, which means that interference willonly occur between pulses which have traveled nearly identical paths.Therefore, pulse i,j will interfere with pulse j,i only. As statedabove, there are N interferometers of interest (pulse 0,i interferingwith pulse i,0 for i=1 to N). There are also the many otherinterferometers which do not include the common rung (rung 0) 712 (e.g.,pulse 1,2 interfering with pulse 2,1, pulse 1,3 interfering with pulse3,1, etc.). Such interfering pulses contribute noise to the usefulpulses, and shall be referred to herein as noise pulses. These noisepulses carry two types of noise. As with all pulses, they carryadditional shot noise, ASE-signal beat noise (in an amplified array),phase noise, etc., which increase the detected noise. The noise pulseswhich form an unwanted interferometer (pulse 1,2 interfering with pulse2,1, etc.) also carry intensity modulation due to interferometricsensing of acoustic waves. This intensity modulation is an unwantedsignal and can be viewed as a source of noise. It is important to notethat these unwanted interferometers have as their interfering pointcouplers 780(1) through 780(N) where the rungs 718(1) through 718(N)couple to the first input/output fiber 714 of the array 716, whereas thesignal pulses interfere at the 3×3 coupler 710. Because the noise pulsesinterfere before they reach the 3×3 coupler 710, the intensitymodulation of the noise pulses is provided symmetrically to bothdetectors 704 and 706. The signal pulses which interfere at the 3×3coupler 710 however produce an asymmetric intensity modulation.Therefore, by differentially amplifying the currents from the detectors704, 706, the intensity modulation of the signal pulses adds and theintensity modulation of the noise pulses subtracts, thus reducing thenoise contribution of the unwanted interferometers.

To completely eliminate all the noise added by these noise pulses, thepulses of interest can be separated from the noise pulses by using atime division multiplexing scheme and properly choosing delay lengths.In particular, the optical path length from the 3×3 coupler 710 throughthe common rung 712 to the 2×2 coupler 730 is selected to correspond toa propagation time τ. The optical path length of a fiber portion fromthe 3×3 coupler to the coupler 780(1), through the first rung 718(1), toa corresponding coupler 790(1) and to the 2×2 coupler 730 is selected tobe (N+1)τ. A portion of the optical path length is a common path fromthe 3×3 coupler 710 to the coupler 780(1) and from the coupler 790(1) tothe 2×2 coupler 730, and a portion of the optical path length is throughthe rung 718(1). The optical path lengths through each of the rungs718(i) are preferably selected to be approximately equal. The totallength of the optical path from the coupler 780(1) to the coupler 780(2)and the optical path from a coupler 790(2) to the coupler 790(1) isselected to be τ such the total optical path length from the 3×3 coupler710 to the 2×2 coupler 730 through the second 718(2) is τ longer thanthe total optical path length from the 3×3 coupler 710 to the 2×2coupler 730 through the first rung 718(1) (i.e., the total optic pathlength between the two couplers 710, 730 through the second rung 718( )is (N+2)τ). The total additional optical path length for each successiverung is elected to be τ. Thus, the travel time of light from the 3×3coupler 710 through a rung 718(i) to the 2×2 coupler 730 is defined asthe delay time T_(i) of the rung 718(i).

In accordance with the foregoing description, T_(i) is determined by theoptical path lengths through the rungs as follows:

T_(i)=τ i=0 (for the common rung 712)

T_(i)=(N+i)τ 1≦i≦N (for each of the sensing rungs 718(1), 718(2), etc.

From the foregoing, it can be seen that the optical path length throughthe farthest rung N is (N+N), or 2Nτ.

The duration of each pulse is selected to be no more than τ. Thus, asillustrated in FIG. 18, the first pulse 800 returned to the 3×3 coupler710 will be the pulse which traveled through the common rung 712 (i.e.,rung 0) from the source 702 to the reflector 752 and back to thedetectors 704, 706. This pulse has a total propagation time of 2τ. (Incomparing propagation times, the propagation time of each pulse to thereflector 752 through the delay loop 750 and back is ignored because thepropagation time is common to all pulses and simply operates as anoffset (not shown) to the timing diagram in FIG. 18.) The next set 810of pulses returned to the detectors 702, 706 are the pulses which travelthrough the common rung 712 in one direction and travel through asensing rung 718(i) in the opposite direction (i.e., the pulses 0,1 and1,0; 0,2 and 2,0; 0,3 and 3,0, through 0,N and N,0). These pulses haverespective propagation times of 2τ+Nτ, 3τ+Nτ, 4τ+Nτ, through (N+1)τ+Nτ.Thus, all the useful pulses are received between a time (N+2)τ and atime (2N+2)τ (including the duration τ of the last pulse received). Incontrast, the interfering pulses which travel through a sensing rung718(i) in both directions (i.e., the pulses, 1,1, 1,2 and 2,1, 1,3 and3,1 . . . 2,2, 2,3 and 3,2, . . . etc.) are received as a set of pulses820 between a time 2(N+2)τ and a time (4N+1)τ. Thus, the signal pulsesare separated from the noise pulses.

For example, in FIG. 18, the number of returned pulses as a function oftime is plotted for N=50. As illustrated, a single pulse is received ata time 2τ. Thereafter, no pulses are received during the interval 3τthrough 52τ. Then, from 52τ through 102τ, two pulses are received duringeach time interval. The noise pulses then return from a time 102τ to atime 201τ. In this way, the signal pulses are separated in time from thenoise pulses, thus preventing the noise pulses from adding noise to thesignal pulses. The electronics (not shown) are readily synchronized toonly look at the pulses received between the time 52τ and the time 102τ.

It should be noted that the source 702 can be activated to send out thenext pulse at the at a time interval of 150τ relative to the previouspulse because the 0τ to 50τ interval in response to the next pulse canoverlap the 150τ to 200τ interval of noise pulses returning in responseto the previous source pulse. Thus, a next set 830 of useful pulses canbegin arriving at a time 201. Therefore, the embodiment of FIGS. 17 and18 has an overall duty cycle of roughly ⅓ for useable signalinformation.

The advantage of the folded Sagnac acoustic fiber sensor 700 over theSagnac loop illustrated in the previous figures is that the delay fiber750 is insensitive to modulation. Because the downleads are often quitelong and are subjected to large movements and vibrations, distributeddownlead pickup is a potentially serious limitation to a Sagnac acousticfiber sensor. In the folded Sagnac acoustic fiber sensor 700, the source708 and detector downleads 770, 772 are insensitive because they occuroutside the interferometer. The delay loop downlead 740 is insensitivebecause all the interfering pulses travel this same fiber separated bysmall time delays (approximately 1 microsecond) and thus see the sameperturbations. Any low frequency (much less than approximately 1 MHz)modulation to the delay loop downlead and delay loop itself is seensubstantially equally by both interfering pulses and thus does notcontribute to a phase difference. The array portion 716 and the commonrung 712 comprise the only sensitive fibers in the interferometer 700.

As shown in FIG. 17, the remotely pumped distributed erbium doped fiberamplifiers (EDFAs) 724 can be located throughout the array 716 toregenerate power, as discussed above.

The 3×3 coupler 710 is used to passively bias each sensor 722(i) nearquadrature and to allow source noise subtraction. Noise subtractionresults from the fact that each detector 704, 706 is biased on anopposite slope (because of the way the signals coming out of the 3×3coupler 710 are phased with respect to each other), causing phasemodulation to asymmetrically affect the intensity at each detector,while source excess noise symmetrically affects the intensity at eachdetector. Therefore, by differentially amplifying the detector outputs,the phase modulation induced intensity variations are added and thesource's intensity noise is subtracted in the same manner that thesignals from the unwanted interferometers would be subtracted.

It should be understood with respect to FIGS. 17 and 18 that a similartime divisional multiplexing effect can be accomplished by providing alonger optical path length through the common rung 712 and shorteroptical path lengths through the sensing rungs 718(i). For example, thecommon rung 712 can advantageously be selected to have an optical pathlength of 2Nτ (i.e., T₀=2N), and the optical paths through the rungs canadvantageously be selected to be τ, 2τ, 3τ, . . . Nτ. The foregoing canbe summarized as:

T_(i)=2Nτ i=0 (for the common rung 712)

T_(i)=iτ 1≦i≦N (for each of the sensing rungs 718(1), 718(2), etc.

Thus, the first signal to return will have an optical propagation time(again subtracting out the propagation time through the delay loop 750which is common to all signals) of 2τ which is the time required to passthrough the first rung 718(1) in both directions. The longest delay ofany signal which passes through one of the sensing rungs 718(i) in bothdirections is 2N for a signal pulse which travels both directionsthrough the farthest sensing rung 718(N). The first useable signal toreturn is a signal which results from the interference of a signal whichtravels in to the reflector 752 through the common rung 712 and returnsthrough the first sensing rung 718(1) with a signal which travels to thereflector 752 through the first sensing rung 718(1) and returns throughthe common rung 712. The interference signal will arrive at a time(2N+1)τ which is later than the last unwanted signal. The last useablesignal will arrive at a time (2N+N)τ (i.e., 3Nτ). Finally, a signalproduced by a pulse which traveled to and from the reflector 752 in thecommon rung 712 arrives at a time 4Nτ, which is well separated from theuseable interference signals.

It is desirable for acoustic sensors to have as large a dynamic range(range of detectable acoustic modulation amplitudes) as possible.Without using demodulation techniques such as the phase-generatedcarrier scheme, the minimum detectable phase modulation is set by thenoise performance of the array, and the maximum detectable phasemodulation (approximately 1 rad) is set by the nonlinear responsefunction of an interferometer. In a Mach-Zehnder sensor, the mapping ofacoustic modulation to phase modulation is a function of only thehydrophone's responsivity. Thus, these limits on the detectable phasemodulation along with this mapping of acoustic modulation into phasemodulation give the range of acoustic modulation the sensor can detect.

In a folded Sagnac acoustic fiber sensor array, the mapping of acousticmodulation into phase modulation is a function of both the responsivityof each of the hydrophones (sensors) 722(i) and the length of the delayloop 750. Thus by changing the length of the delay loop 750, the dynamicrange of the sensors 722(i) can be adjusted without modifying thehydrophones 722(i) themselves. In addition, if two reflectors 752(1) and752(2) are used, each sensor 718(i) can have two different delay loops750(1) and 750(2), as shown in a sensor 850 in FIG. 19. This allows eachsensor 722(i) to return two signals which have different dynamicsranges, as discussed above with respect to FIGS. 7 and 8, therebygreatly increasing the total dynamic range of each sensor 722(i). Thepenalty is a reduction in duty cycle for each individual signal by afactor of 1/(number of delay loops).

FIG. 20 illustrates a sensor 900 which implements a phase-nullingtechnique similar to techniques which have been used in fibergyroscopes. The delay loop reflector 752 of FIG. 17 is not used in thesensor 900 of FIG. 20. Rather, the pulses are instead returned via areturn downlead 910 into the previously unused port of the 2×2 coupler730. An optical isolator 912 is inserted in the return downlead 910 toprevent light from traveling the delay loop 750 in both directions. Thesensor 900 of FIG. 20 behaves identically to the sensor 700 of FIG. 17with the reflector 752. However, the sensor 900 allows the addition of aphase modulator 920 to be inserted into the return downlead 910. Thephase modulator 920 is activated to add a phase shift to each pulseindividually. By feeding the detected phase shift into the phasemodulator 920 via a differential amplifier 922, phase changes are nulledout, and the required applied phase shift in the phase modulator 920becomes the signal. In this phase nulling method, the dynamic range ofthe array 900 is limited only by the maximum phase shift that the phasemodulator 920 can provide.

FIG. 21 illustrates a further alternative embodiment of FIG. 19 in whichthe two delay loops 750(1) and 750(2) are not connected to the samedelay loop downlead. Rather, the first end of the first delay loop750(1) is connected to a first delay loop downlead 740(1) which isconnected to the fourth port of the 2×2 coupler 730 as in FIG. 19. Thesecond end of the first delay loop 750(1) is coupled to the firstreflector 752(1) as before. The first end of the second delay loop750(2) is coupled to the third port of the 2×2 coupler 730 via a seconddelay loop downlead 740(2), and the second end of the second delay loop750(2) is coupled to the second reflector 752(2). Approximately half thelight from the 2×2 coupler 730 is coupled to each of the downleads740(1), 740(2). The light in each downlead 740(1), 740(2) is delayed inthe respective delay loop 750(1), 750(2) and is reflected back to the2×2 coupler 730 as before. The reflected light is coupled to the commonrung 712 and to the array 716. The delays of the delay loops 750(1),750(2) are selected so none of the N+1 pulses which propagate from thefourth port of the 2×2 coupler 730 through the first delay loop 750(1)overlap in time with any of the N+1 pulses which propagate from thethird port of the 2×2 coupler 730 through the second delay loop 750(2).Thus, the embodiment of FIG. 21 provides similar functionality to theembodiment of FIG. 19; however, the embodiment of FIG. 21 utilizes thelight which was coupled out of the third port of the 2×2 coupler 730 inFIG. 19 and discarded.

FIG. 22 illustrates an alternative embodiment of a fiber optic acousticsensor system 1000 using a folded Sagnac sensor array. In the system1000, a source 1004 is coupled to a first port of a 2×2 polarizationmaintaining coupler 1006 by an X-polarizer 1008. A detector 1002 isconnected to a second port of the 2×2 coupler 1006 via a X-polarizer1010. A second detector (not shown) may advantageously be included inthe embodiment of FIG. 22 by coupling light from the fiber leading tothe source 1004. The X-polarizer 1008 only passes light from the source1004 having a first polarization (e.g., an X-polarization). Thus, thepolarization maintaining coupler 1006 receives light having anX-polarization from the source 1004 and couples the light to a commonrung 1020 via a third port and to a sensor array 1022 via a fourth port.The sensor array 1022 has a similar structure to the sensor array 716 ofFIG. 17, and like elements have been numbered accordingly.

Note that the two X-polarizers 1008, 1010 can be replaced by one or moreX-polarizers in alternative locations in the system 1000.

The common rung 1020 is coupled via an X-polarizer 1030 to a first portof a second polarization maintaining 2×2 coupler 1032. The lightpropagating to the array 1022 first passes through a depolarizer 1034and then to the first input/output fiber 714. The depolarizer 1034couples substantially equal amounts of the X polarized light to Xpolarized light and to Y polarized light. Thus, approximately 50 percentof the light propagates in the array 1022 as X-polarized light, andapproximately 50 percent propagates in the array 1022 as Y-polarizedlight.

After passing through the rungs of the array 1022, the light propagatesvia the second input/output fiber 720 and a Y-polarizer 1040 to a secondport of the second coupler 1032. The Y-polarizer 1040 allows onlyY-polarized light to enter the second coupler 1032. The coupler 1032combines the light from the array 1022 and from the common rung 1020Approximately half the light entering the coupler 1032 is coupled via athird port of the coupler 1032 to a light absorbing termination 1042,and approximately half of the light is coupled to a downlead fiber 1050which propagates the light to a first end of a delay loop 1052.

Light passes through the delay loop 1052 to a Faraday rotating mirror(FRM) 1054. The operation of the Faraday rotating mirror 1054 is wellknown and will not be described in detail. Basically, when light isincident onto the Faraday rotating mirror 1054 in one polarization, itis reflected in the orthogonal polarization. Thus, the X-polarized lightwhich passed through the common rung 1020 is reflected as Y-polarizedlight, and the Y-polarized light which passed through the array isreflected as X-polarized light.

The reflected light passes back through the delay 1052 and enters thefourth port of the coupler 1032. The light is coupled to the common rung1020 and to the array 1022. The X-polarizer 1030 in the common rungpasses only the light in the X-polarization which originally propagatedthrough the array 1022. Similarly, the Y-polarizer 1040 in the array1022 passes only Y-polarized light which originally propagated throughthe common rung 1020.

After propagating through the array 1022, the returning Y-polarizedlight is depolarized in the depolarizer 1034 to produce both X-polarizedlight and Y-polarized light. The light from the common rung 1020 entersthe third port of the coupler 1006, and light from the depolarizer 1034enters the fourth port of the coupler 1006. The light combines in thecoupler, and the X-polarized light from the two ports which has traveledthe same optical distance interferes and is coupled to the first andsecond ports. The portion coupled to the second port propagates throughthe X-polarizer 1010 to the detector 1002 where the interfering signalsare detected.

It should be understood that only the light which originally traveleddifferent paths to and from the Faraday rotating mirror 1054 interferesat the coupler 1006. The only light allowed to propagate through thecommon rung 1020 in the reflected direction is X-polarized light whichoriginally propagated in the array 1022 as Y-polarized light. Similarly,the only light allowed to propagate through any of the rungs of thearray 1022 in the reflected direction is Y-polarized light whichoriginally propagated in the common rung 1020 as X-polarized light.Potentially interfering light cannot travel through the rungs in bothdirections to produce the noise signals described above in connectionwith the above-described embodiments. Thus, each of the pulses generatedin the array 1022 from the reflected pulse that originally traveled inthe common rung 1020 can interfere with only a single one of the pulseswhich was originally generated in the array 1022 and which propagated inthe common rung 1020 after it was reflected. Thus, it is not necessaryin the embodiment of FIG. 22 to include additional delays to separatethe useable signal pulses from noise pulses.

FIGS. 23A, 23B and 23C illustrate further alternative embodiments of thepresent invention. A sensor array 1100 in the embodiments of FIGS. 23A,23B and 23C is similar to the sensor array 700 in the embodiment of FIG.17, and like elements have been numbered accordingly. The embodiments ofFIGS. 23A, 23B and 23C include an unpolarized source 1102. The 2×2coupler 730 of FIG. 17 is replaced with a polarization beam splitter(PBS) 1104 in FIGS. 23A, 23B and 23C. The use of the polarization beamsplitter 1104 saves approximately 6 dB of power compared to the coupler730 in FIG. 17 and the coupler 1130 in FIG. 22. The reflector 752 inFIG. 17 is replaced with a Faraday rotating mirror (FRM) 1106, which issimilar to the Faraday rotating mirror 1054 of FIG. 22. The 3×3 coupler710 in FIGS. 23A, 23B and 23C does not have to be a polarizationmaintaining coupler.

Each of FIGS. 23A, 23B and 23C includes a depolarizer 1110. In FIG. 23A,the depolarizer 1110 is located on the first array input/output fiber714. In FIG. 23B, the depolarizer 1110 is located on the common rung712. In FIG. 23C, the depolarizer 1110 is located on the second arrayinput/output fiber 720.

In the embodiment of FIG. 23A, light from the unpolarized source 1102enters the 3×3 coupler 710 and is coupled in approximately equalportions to the common rung 712 and to the first array input/outputfiber 714. As discussed above in connection with FIGS. 3 and 17, the useof the 3×3 coupler provides passive biasing near quadrature. The lightpropagating in the first array input/output fiber 714 passes through thedepolarizer 1110, which has the effect of causing substantially half ofthe light entering the array in one polarization (e.g., theX-polarization) to be coupled into the orthogonal polarization (e.g.,the Y-polarization), and likewise half of the light entering the arrayin the Y-polarization to be coupled to the X-polarization. Thus, afterthe depolarizer 1110, half of the light in the X-polarization originatedin the X-polarization and the other half of the light in theX-polarization originated in the Y-polarization. Likewise, after thedepolarizer 1110, half of the light in the Y-polarization originated inthe Y-polarization and the other half of the light in the Y-polarizationoriginated in the X-polarization. Effectively, the depolarizer 1110scrambles the unpolarized light.

The light passes through the array 716 in the manner described above inconnection with the other embodiments. The light exiting the array 716propagates through the second array input/output fiber 720 to a firstport 1121 of the polarization beam splitter 1104. The polarization beamsplitter 1104 splits the incident light into the two orthogonalpolarizations (i.e., the X-polarization and the Y-polarization). For thepurpose of this discussion, it is assumed that the polarization beamsplitter 1104 operates like a polarization-dependent mirror oriented at45°, wherein light entering the first port 1121 in one polarization(e.g., the X-polarization) is reflected to a second port 1122 and lightentering the first port 1121 in the other polarization (e.g., theY-polarization) is transmitted to a third port 1123. In the embodimentshown, the light exiting the second port 1122 is nonreflectivelyabsorbed by the terminator 732. The Y-polarized light exiting the thirdterminal 1123 propagates through the delay loop downlead fiber 740,through the delay loop 750 to the Faraday rotating mirror 1106. Notethat this Y-polarized light from the array portion 716 traveled throughthe depolarizer 1110 and half of it was originally X-polarized light andhalf of it was originally Y-polarized light. As discussed above, theFaraday rotating mirror 1106 causes the incident light to be coupled tothe orthogonal polarization. Thus, the Y-polarized light is coupled tothe X-polarization.

The X-polarized light reflected by the Faraday rotating mirror 1106passes through the delay loop 750 and the delay loop downlead fiber 740back to the third port 1123 of the polarization beam splitter. Becausethe light is now in the X-polarization, the light is reflected to afourth port 1124 rather than being transmitted to the first port 1121.Thus, the Y-polarized light which was originally incident on thepolarization beam splitter from the array 716 is coupled to the commonrung 712 to propagate back to the 3×3 coupler 710 in the X-polarization.

Unpolarized light which propagates from the 3×3 coupler 710 to thepolarization beam splitter 1104 via the common rung 712 enters thepolarization beam splitter 1104 via the fourth port 1124. The componentsof the light in the Y-polarization are transmitted to the second port1122 and are nonreflectively terminated by the terminator 732. Thecomponents of the light in the X-polarization are reflected to the thirdport 1123 and propagate to the Faraday rotating mirror 1106 via thedelay loop downlead fiber 740 and the delay loop 750. (The reason forincluding the depolarizer 1110 can now be understood. Because only theX-polarized light from the common rung 712 is coupled to the delay loopdownlead fiber 740, the depolarizer 1110 ensures that the light coupledfrom the array 716 to the delay loop downlead fiber 740 also includessome light which was originally X-polarized.) The Faraday rotatingmirror 1106 reflects the light as Y-polarized light, and the Y-polarizedlight propagates through the delay loop and the downlead fiber to thethird port 1123 of the polarization beam splitter 1104.

The Y-polarized light incident on the third port 1123 of thepolarization beam splitter 1104 is transmitted to the first port 1121and thus to the second array input/output fiber 720. The Y-polarizedlight propagates through the array 716 to the first array input/outputfiber 714 and then passes through the depolarizer 1110 to the 3×3coupler 710. The depolarizer 1110 operates to convert approximately 50percent of the Y-polarized light to X-polarized light. The X-polarizedlight from the depolarizer 1110 interferes with the X-polarized lightfrom the common rung 712. The resulting combined light is detected bythe detector 704 or the detector 706 in accordance with the phaserelationship between the interfering light signals in the 3×3 coupler710.

Note that the X-polarized light incident on the 3×3 coupler 710 from thedepolarizer 1110 and the X-polarized light from the common rung 712travel identical path lengths. For example, light which propagatesthrough the common rung 712 first, propagates in the X-polarizationthrough the common rung 712 and then propagates through the array 716 inthe Y-polarization. On the other hand, the light which propagatesthrough the array 716 first propagates in the Y-polarization through thearray 716 and then propagates in the X-polarization through the commonrung. Because the two “counterpropagating” light signals are in the samepolarizations when propagating through the corresponding portions of theinterferometric path, the propagation lengths are identical except forthe effect of incident noise sensed by the array 716.

It should be understood that the terminator 732 coupled to the secondport 1122 of the polarization beam splitter 1104 can be replaced with asecond delay loop (not shown) and a second Faraday rotating mirror (notshown) to provide a second interferometric path for light whichinterferes in the Y polarization. By adjusting the delay provided by thesecond delay loop, the return signals from the second interferometricpath can be precluded from overlapping with the return signals from thefirst interferometric path.

The embodiment of FIG. 23B is similar to the embodiment of FIG. 23Aexcept that the depolarizer 1110 is positioned in the common rung 712.The effect of the depolarizer 1110 in FIG. 23B is (1) to cause a portionof the light in the common rung 712 returning from the polarization beamsplitter 1104 in a single polarization (e.g., the X-polarization) to becoupled to the orthogonal polarization and (2) to scramble theunpolarized light which travels from the 3×3 coupler 710 through thecommon rung 712 towards the polarization beam splitter 1104. Thisensures that the light interferes when it recombines at the 3×3 coupler710 (the same reason the depolarizer 1110 was added to the fiber 714 ofFIG. 23A).

The embodiment of FIG. 23C is also similar to the embodiment of FIG. 23Aexcept that the depolarizer 1110 is positioned in the second arrayinput/output fiber 720. The embodiment of FIG. 23C is functionallyequivalent to the embodiment of FIG. 23A because it does not matterwhether the light passes through the array portion 716 and then passesthrough the depolarizer 1110 or passes through the depolarizer 1110 andthen passes through the array portion 716. Thus, the function of theembodiment of FIG. 23C is substantially the same as the function of theembodiment of FIG. 23A, as described above.

FIG. 24 illustrates a further alternative embodiment of the presentinvention in which a folded Sagnac sensor array 1200 includes thepolarization beam splitter (PBS) 1104, the Faraday rotating mirror (FRM)1106, and the depolarizer 1110 connected as shown in the array 1100 inFIG. 23A. Other components from FIG. 23A are also numbered as before.Unlike the array 1100 in FIG. 23A which has the 3×3 coupler 710, thefolded Sagnac sensor array 1200 has a polarization maintaining (PM) 2×2coupler 1220 connected in like manner as the 2×2 coupler 1006 in FIG.22. One port of the 2×2 coupler 1220 is connected to a first port of anoptical circulator 1222 via a first polarizer 1224. A second port of theoptical circulator 1222 is connected to a first detector 1226. A thirdport of the optical circulator 1222 is connected to an unpolarizedsource 1228 (e.g., an intensity modulated fiber superfluorescentsource). A second port of the 2×2 coupler 1220 is connected to a seconddetector 1230 via a second polarizer 1232. The detectors 1226 and 1230and the unpolarized source 1228 are connected to the circulator 1222 bystandard (not polarization maintaining) fibers. The polarizers 1224 and1232 are coupled to the polarization maintaining coupler 1220 via thepolarization maintaining fiber such that the polarizers 1224, 1232 arealigned with a same axis of the polarization maintaining 2×2 coupler1220. Alternatively, if a polarized source is used in place of theunpolarized source 1228, the polarized source (not shown) is connectedto a polarization maintaining circulator (not shown) by polarizationmaintaining fiber, and the polarization maintaining circulator isconnected to the polarizer 1224 by polarization maintaining fiber. Thepolarization maintaining components are connected such that thepolarized light from the source passes through the polarizer 1224. Theconnections from the polarization maintaining circulator to thedetectors 1226 and 1230 are provided by standard (not polarizationmaintaining) fibers.

The folded Sagnac sensor array 1200 further includes a non-reciprocalphase shifter 1250. The phase shifter 1250 is coupled to the common rung712 via a first optical fiber 1252 having a first end 1254 and a secondend 1256 and via a second optical fiber 1258 having a first end 1260 anda second end 1262. The first end 1254 of the first optical fiber 1252 iscoupled to the common rung 712 proximate to the 2×2 coupler 1220 via afirst coupler 1264. The first end 1260 of the second optical fiber 1258is coupled to the common rung 712 proximate to the polarization beamsplitter 1104 via a second coupler 1266. The respective second ends1256, 1262 of the first and second optical fibers 1252, 1258 are coupledto the phase shifter 1250, as discussed below in connection with FIGS.25 and 26.

Preferably, the common rung 712, the first fiber 1252 and the secondfiber 1258 are polarization maintaining (PM) fibers, and the firstcoupler 1264, the second coupler 1266 and the 2×2 coupler 1220 arepolarization maintaining (PM) couplers. Also, preferably, the firstcoupler 1264 and the second coupler 1266 are 50/50 couplers which coupleapproximately 50 percent of the light entering the common rung 712 ineither direction to the phase shifter 1250 while approximately 50percent of the light remains in the common rung. Thus, thenon-reciprocal phase shifter 1250 and the associated fibers form asecond rung 1268 in parallel with the common rung 712.

Preferably, one of the rungs 712, 1268 (e.g., the common rung 712)includes a delay element (e.g., a delay loop 1269) that introduces atime delay in one rung sufficient to prevent the pulses propagatingthrough the rungs from overlapping. Thus, the light returning to 2×2coupler 1220 from the sensor array 716 comprises two pulses for eachsensor that are spaced apart in time from each other. One pulsecomprises the combined light that passes through the common rung 712 ineach direction. The other pulse comprises the combined light that passesthrough the non-reciprocal phase shifter 1250 in each direction. Itshould be understood that the light pulse that passes through the phaseshifter 1250 in one direction and the light pulse that passes throughthe common rung 712 in the other direction have substantially differentpropagation times and will not overlap in the coupler 1220. Thus, theywill not interfere.

The light that passes though the common rung 712 in one direction doesnot undergo any phase shift within the common rung 712 relative to thelight that passes through the common rung in the other direction. Thus,the combined light that passes through the common rung 712 in bothdirections has a relative phase bias of zero. However, as discussedbelow, the non-reciprocal phase shifter 1250 does introduce a shift ofthe light in one direction with respect to the light in the otherdirection. In particular, in a preferred embodiment, the phase shifter1250 introduces a relative π/2 phase shift between the light in the twodirections. Thus, the light entering the coupler 1220 that haspropagated through the phase shifter 1250 in both directions willcombine in the coupler 1220 with a π/2 phase bias.

One skilled in the art will appreciate that the 50 percent coupler 1220in the interferometric configuration shown in FIG. 24 couples returninglight to the output port corresponding to the original input port whenthe returning light at the two input ports interferes in the coupler andhas a relative phase difference of 0, 2π, 4π, etc., and couplesreturning light to the other output port when the light has a relativephase difference of π, 3π, 5π, etc. When the returning light has arelative phase difference that is not a multiple of π, a portion of thereturning light is output from both ports. For example, when therelative phase difference is an odd multiple of π/2 (e.g., π/2, 3π/2,etc.), approximately 50 percent of the returning light is coupled toeach output port. By providing two independent propagation paths, eachdetector 1226, 1230 receives two signals that are spaced apart in timeand can therefore be separately detected. One signal has a 0 phase bias,and one signal has a π/2 phase bias so that when one signal is leastsensitive to perturbation, the other signal is most sensitive toperturbation, and vice versa. It should be understood that additionalrungs in parallel with the common rung 712 and with differing amounts ofrelative phase shift can be included to provide pulses with differentphase biasing.

FIG. 25 illustrates an alternative configuration of a folded Sagnacsensor array 1200′, which is substantially similar to the folded Sagnacsensor array 1200 of FIG. 24. In the folded Sagnac sensor array 1200′ ofFIG. 25, the depolarizer 1110 is located in the second arrayinput/output fiber 720 rather than in the first array input/output fiber714. Because of the reciprocal structure of the sensor array 716, therelocation of the depolarizer 1110 to the fiber 720 does not change theoverall operation of the folded Sagnac sensor array 1200′ with respectto the operation of the folded Sagnac sensor array 1200. Thus, theoperation of the folded Sagnac sensor array 1200′ will not be describedin detail herein.

The embodiments of FIGS. 24 and 25 include the sensor array 716, whichwas described in detail above. It should be understood that otherconfigurations of amplified sensor arrays can also be used in place ofthe sensor array 716 in the embodiments of FIGS. 24 and 25.

FIG. 26 illustrates a first preferred embodiment of the non-reciprocalπ/2 phase shifter 1250 of FIGS. 24 and 25. As illustrated in FIG. 26,the phase shifter 1250 comprises a first collimating lens 1270, a first45° Faraday rotator 1272, a quarter-wave plate 1274, a second 45°Faraday rotator 1276, and a second collimating lens 1278. In theillustrated embodiment, the first Faraday rotator 1272, the secondFaraday rotator 1276 and the quarter-wave plate 1274 comprise bulk opticdevices that are commercially available, but may advantageously comprisefiber optic or other waveguide devices. The collimating lenses 1270,1278 are positioned proximate to second ends 1256, 1262 of the PM fibers1252, 1258 to focus light from the fiber ends 1256, 1262 onto theFaraday rotators 1272, 1276, respectively, and to focus light from theFaraday rotators 1272, 1276 into the fiber ends 1256, 1262. Each of theFaraday rotators 1272, 1276 operates in a well-known manner to causelight input to the Faraday rotator with its polarization at a particularangle to have the polarization rotated so that the polarization is at anew angle rotated by a predetermined amount with respect to the originalangle. For example, in the preferred embodiment, each Faraday rotator1272, 1276 rotates the polarization of the incident light by 45° in thecounterclockwise (ccw) direction. Thus, as illustrated in FIG. 26, lightemitted from the end 1256 of the PM fiber 1252 having its polarizationoriented horizontally, will be rotated by 45° counterclockwise in thefirst Faraday rotator 1272 such that the polarization is oriented at anangle of 45° in the clockwise direction with respect to the originalorientation when it emerges from the first Faraday rotator 1272.

The quarter-wave plate 1274 is positioned between the two Faradayrotators 1272, 1276. The quarter-wave plate 1274 has a firstbirefringent axis 1280 and an orthogonal second birefringent axis 1282.Light propagating in a polarization oriented along one birefringent axis(e.g., the first birefringent axis 1280) has a slower propagationvelocity than the light propagating in a polarization oriented along theother birefringent axis (e.g., the second birefringent axis 1282). Thequarter-wave plate 1274 is oriented so that the first birefringent axis1280, for example, is oriented at 45° in the clockwise direction to thevertical, and is therefore oriented so that the light emerging from thefirst Faraday rotator 1272 is oriented along the first birefringent axis1280 and is orthogonal to the second birefringent axis 1282. Because ofthe difference in propagation velocities along the two axes, thequarter-wave plate 1274 introduces a π/2 or 90° phase shift in the lightpolarized along the first birefringent axis 1280 with respect to thelight polarized along the second birefringent axis 1282. Thus, inaccordance with this example, the light that originally propagated inthe horizontal polarization that was rotated to be in alignment with thefirst birefringent axis 1280 incurs a relative phase shift of 90° withrespect to any light that propagates along the second birefringent axis1282.

After passing through the quarter-wave plate 1274, the light passesthrough the second Faraday rotator 1276 and is again rotated by 45° inthe counterclockwise direction. The light emerging from the secondFaraday rotator 1276 passes through the second collimating lens 1278 andis focused into the second end 1262 of the second PM optical fiber 1258.It should be understood from the foregoing description that any lightoutput from the first PM optical fiber 1252 in the horizontalpolarization enters the second PM optical fiber 1258 in the verticalpolarization. As discussed above, the light entering the second PMoptical fiber 1258 in the vertical polarization will have propagatedalong the slow birefringent axis 1280 of the quarter-wave plate 1274 andwill incur a relative π/2 phase difference with respect to light thatpropagates along the fast birefringent axis 1282.

As indicated by its description, the non-reciprocal phase shifter 1250operates in a non-reciprocal manner because of the operation of theFaraday rotators 1272, 1276. As described above, the light passingthrough the Faraday rotators 1272, 1276 from the first PM fiber 1252 tothe second PM fiber 1258 is rotated 45° counterclockwise by each rotatorwith respect to the direction of propagation of the light shown in FIG.25. If the Faraday rotators were reciprocal, light propagating throughthe Faraday rotators 1272, 1276 in the opposite direction would also berotated in the counterclockwise direction with respect to the directionof propagation of the light; however, because the Faraday rotators arenon-reciprocal, the light is rotated in the opposite direction (i.e.,clockwise with respect to the propagation direction of the light). Thenon-reciprocal effect is illustrated in FIG. 27 for light passing fromthe second end 1262 of the second PM fiber 1258, through thenon-reciprocal phase shifter 1250, to the second end 1256 of the firstPM fiber 1252. Note that when viewed as in FIG. 27, the rotation appearsto again be in the counterclockwise direction; however, the light is nowpropagating toward the viewer. Thus, light emitted from the second end1262 of the second PM optical fiber 1258 in the vertical polarizationpasses through the second collimating lens 1278 and through the secondFaraday rotator 1276 and is rotated to an orientation in alignment withthe second (fast) birefringent axis 1282 of the quarter-wave plate 1274.Thus, the light originally in the vertical polarization does notexperience a relative delay as it propagates through the quarter-waveplate 1274. After passing through the quarter-wave plate, 1274, thelight passes through the first Faraday rotator 1272 such that the lightis rotated an additional 45° to the horizontal polarization. The lightis then focused through the first collimating lens 1270 onto the secondend 1256 of the first PM optical fiber 1252.

From the foregoing, it can be seen that the horizontally polarized lightpassing in the first direction from the first PM fiber 1252 to thesecond PM fiber 1258 via the non-reciprocal phase shifter 1250propagates through the slow birefringent axis 1280 of the quarter-waveplate 1274 and experiences a relative phase delay of 90° or π/2. Thehorizontally polarized light propagating in the first direction isrotated such that that the light is oriented in the verticalpolarization when it enters the second PM fiber 1258. Conversely,vertically polarized light passing from the second PM fiber 1258 to thefirst PM fiber 1252 via the non-reciprocal phase shifter 1250 in thesecond direction propagates through the fast birefringent axis 1282 ofthe quarter-wave plate 1274 and does not experience a relative phasedelay. The vertically polarized light propagating in the seconddirection is rotated such that the light is oriented in the horizontalpolarization when it enters the first PM fiber 1252. As will bediscussed more fully below, the relative phase shift between thehorizontally polarized light propagating in the first direction withrespect to the vertically polarized light propagating in the seconddirection provides a π/2 phase bias.

FIGS. 28 and 29 illustrate an alternative embodiment of thenon-reciprocal phase shifter 1250, in which the first Faraday rotator1272 is positioned between the quarter-wave plate 1274 (now referred toas the first quarter-wave plate) and a second quarter-wave plate 1294.In FIG. 28, light from the second end 1256 of the first PM fiber 1252 iscollimated by the first collimating lens 1270, as before. The light isoriginally in the horizontal polarization. When the light passes throughthe first quarter-wave plate 1274, it is converted to light having acircular polarization. The circular polarized light passes through thefirst Faraday rotator 1272, which causes the circularly polarized lightto incur a phase shift of φ. In the preferred embodiment, the firstFaraday rotator 1272 is selected to cause a phase shift of π/4. Thelight from the Faraday rotator 1272 remains circularly polarized andpasses through the second quarter-wave plate 1294, which converts thecircularly polarized light to linearly polarized light in the verticalpolarization orientation. In addition to being in the verticalpolarization, the light has experienced a phase shift of φ (e.g., π/4).

FIG. 29 illustrates the operation of the alternative embodiment of thenon-reciprocal phase shifter 1250 for light propagating in the oppositedirection. In FIG. 29, vertically polarized light from the second end1262 of the second PM fiber 1260 is collimated by the second collimatinglens 1278 and passes through the second quarter-wave plate 1294. Thesecond quarter-wave plate 1294 converts the vertically polarized lightto light having a circular polarization. The circularly polarized lightpasses through the first Faraday rotator 1272 and experiences a phaseshift as before. Because the light is propagating through the firstFaraday rotator 1272 in the opposite direction, the light experiences anopposite phase shift of −φ (e.g., −π/4). The light from the firstFaraday rotator 1272 then passes through the first quarter-wave plate1274, where the circularly polarized light is converted to linearlypolarized light with a horizontal polarization. Thus, the lightpropagating in the two directions experiences a total relative phaseshift of 2φ (e.g., π/2), which has the same effect as the firstembodiment of the non-reciprocal phase shifter 1250 illustrated in FIGS.26 and 27.

The effect of the non-reciprocal phase shifter 1250 on the orientationof the polarization and the phase delay provides the biasing effectdescribed above and explained again in connection with FIG. 24. As shownin FIG. 24, the light entering the second PM fiber 1258 in the verticalpolarization is combined at the second PM coupler 1266 with the lightthat propagated through the common rung 712 from the first PM coupler1264 to the second PM coupler 1266. For reasons that will becomeapparent in the following discussion, it is desirable that the lightentering the second PM coupler 1266 from the common rung 712 have thesame polarization as the light entering the second PM coupler from thesecond PM fiber 1258. Thus, in the preferred embodiment, either thesecond PM fiber 1258 or the common rung 712 is rotated by 90° so thatthe light in the vertical polarization in the second PM fiber 1258 isoriented in the same direction as the light in the horizontalpolarization of the common rung 712. This is readily accomplished byrotating the second end 1262 of the second PM fiber 1258 proximate tothe second collimating lens 1278 so that the vertically polarized lightenters the second end 1262 with its state of polarization oriented alongthe horizontal polarization axis of the second PM fiber 1258. Thus, thelight that exits the non-reciprocal phase shifter 1250 in the verticalstate of polarization is applied to the coupler 1266 as light in thehorizontal state of polarization with respect to the polarization axesof the coupler 1266. Accordingly, the light from the non-reciprocalphase shifter 1250 has the same state of polarization as the light fromthe common rung 712.

The light that passes through the common rung 712 and the light thatpasses through the non-reciprocal phase shifter 1250 next enter the port1124 of the polarization beam splitter (PBS) 1104. The light in thehorizontal polarization is output from the port 1123 of the PBS 1104 tothe fiber 740. The fiber 740 includes the delay loop 750 and isterminated at the Faraday rotating mirror (FRM) 1106. The delay loop 750and the FRM 1106 operate as discussed above, and the reflected anddelayed pulses are returned to the port 1123 of the PBS 1104 in thevertical polarization. The pulses are output from the port 1121 of thePBS 1104 to the array 716 via the fiber 720 and propagate in theclockwise direction through the sensors 722(i) of the array 716.

The pulses are output from the array 716 via the fiber 714 and thedepolarizer 1110 to the 2×2 coupler 1220 where the clockwise propagatinglight is combined with the counterclockwise propagating light. Thecounterpropagating light also starts out as horizontally polarizedlight. The light is depolarized and passes through the sensor array 716.Light emerging from the sensor array 716 in the vertical polarization isreflected by the PBS 1123 and is discarded via the port 1122 and theterminator 732. Light emerging from the sensor array 716 in thehorizontal polarization passes through the PBS 1123, is delayed by theloop 750, and is rotated to the vertical polarization by the FRM 1106.The return light, which is in the vertical polarization, is reflected bythe PBS 1123 to the port 1124 and is thus directed to the second PMcoupler 1266. A portion of the light passes through the delay loop 1269of the common rung 712 and a portion of the light passes through thenon-reciprocal phase shifter 1250. As discussed above, light enteringthe non-reciprocal phase shifter 1250 in the vertical polarizationpropagates through the fast birefringent axis 1282 of the quarter-waveplate 1274 (FIG. 27) and does not experience a relative phase delay.Thus, the two pulses of counterclockwise light propagate to the coupler1220 where they are combined with the clockwise propagating lightpulses. The light signals that passed through the common rung 712 andthe delay loop 1269 in both directions experience no relative phaseshift and combine as discussed above. The light signals that passedthrough the non-reciprocal phase shifter 1250 in both directionsexperience a relative phase shift of π/2 between the clockwisepropagating signal and the counterclockwise propagating signal and thushave a π/2 phase bias, as discussed above. At both outputs of thecoupler 1220, a portion of the two pulses of light returning from thesensor array 1200 is directed to the polarizer 1224, and the remainingportion is directed to the polarizer 1232. The role of the twopolarizers 1224 and 1232 is to ensure that the light entering the loophas the same polarization as the light leaving the loop, whichguarantees reciprocity. As described earlier, the two pulses reachingthe detector 1230 are in phase quadrature, which allows the use of anumber of signal processing techniques well-known in the art to avoidsignal fading. Similar comments apply to the detector 1226. In theembodiment of FIG. 24, the generation of two pulses in phase quadratureis the main reason for incorporating the rung containing thenon-reciprocal phase shifter 1250.

FIGS. 30-36 illustrate further alternative embodiments of the presentinvention in which a folded Sagnac sensor array utilizespolarization-based biasing for multiple detectors, wherein each detectorhas a bias point which can be set independently of the bias points ofthe other detectors. The embodiments of FIGS. 30-36 include the sensorarray 716, which was described in detail above. It should be understoodthat other configurations of amplified sensor arrays can also be used inplace of the sensor array 716 in the embodiments of FIGS. 30-36.

In a folded Sagnac sensor array 1300 illustrated in FIG. 30, a polarizedfiber superfluorescent source (SFS) 1310 is coupled to a polarizationcontroller 1312 via a fiber 1314. The fiber 1314 further couples thepolarization controller 1312 to a first port of a 2×2 coupler 1316. Asecond port of the coupler 1316 is an output port, which will bediscussed below. A third port of the coupler 1316 is coupled via a fiber1318 to a non-reflective terminator 1320. A fourth port of the coupler1316 is coupled to a first port 1330 of a polarization beam splitter(PBS) 1332 via a common array input/output fiber 1334. A second port1336 of the polarization beam splitter 1332 is coupled to a firsthorizontal polarizer 1338. The first horizontal polarizer 1338 iscoupled to the second array input/output fiber 720 of the array 716. Athird port 1340 of the polarization beam splitter 1332 is connected to acommon delay fiber 1342, which is formed into a delay loop 1344 andwhich is terminated at a Faraday rotating mirror (FRM) 1346. A fourthport 1348 of the polarization beam splitter 1332 is coupled to a secondhorizontal polarizer 1350 and then to a depolarizer 1352. Thedepolarizer 1352 is coupled to the first array input/output fiber 714.

The second port of the coupler 1316 is coupled to a detector subsystem1360 via a fiber 1362. In the embodiment of FIG. 30, the detectorsubsystem 1360 comprises a 1×n coupler 1364 which has a single inputport that receives the light from the second port of the coupler 1316. Afirst output port of the 1×n coupler 1364 is coupled to a polarizationcontroller 1366. The polarization controller 1366 is coupled to apolarizer 1368, which is in turn coupled to a first detector 1370. Asecond output port of the 1×n coupler 1364 is coupled to a polarizationcontroller 1372. The polarization controller 1372 is coupled to apolarizer 1374, which is coupled to a second detector 1376. Additionalpolarization controllers, polarizers and detectors (not shown) can beconnected to additional ports (not shown) of the 1×n coupler 1364.

The folded Sagnac sensor array 1300 of FIG. 30 operates in the followingmanner. The polarized SFS 1310 provides a polarized output signal whichpasses through the polarization controller 1312 via the fiber 1314. Thepolarization controller 1312 is adjustable to vary the polarization to adesired state of polarization. For example, in FIG. 30, the state ofpolarization is adjusted to provide linearly polarized light oriented at45° with respect to the vertical and horizontal axes at the input to thepolarization beam splitter 1332. The light remains in the fiber 1314 andis provided as the input to the coupler 1316. The coupler 1316 couplesapproximately 50 percent of the incoming light to the first output fiber1318 and is thus discarded at the non-reflective terminator 1320. Thecoupler 1316 couples approximately 50 percent of the incoming light tothe common array input/output fiber 1334.

The common array input/output fiber 1334 guides the light to thepolarization beam splitter 1330, which reflects horizontally polarizedlight to the second port 1336 and which passes vertically polarizedlight to the third port 1340. The reflected horizontally polarized lightfrom the second port 1336 passes through the first horizontal polarizer1338 to the second array input/output fiber 720 and propagates in aclockwise direction through the array 716. The clockwise propagatinglight exits the array 716 via the depolarizer 1352 and the arrayinput/output fiber 714. As discussed above, the depolarizer 1352 assuresthat the exiting light is substantially equally distributed in thehorizontal polarization mode and the vertical polarization mode afterpassing through the sensors in the array 716. The clockwise propagatinglight then passes through the second horizontal polarizer 1350, whicheliminates the portion of the light in the vertical polarization. Theclockwise propagating light in the horizontal polarization then entersthe fourth port 1348 of the polarization beam splitter 1330 and isreflected to the third port 1340 to propagate in the common delay fiber1342. The returning clockwise light passes through the delay loop 1344to the Faraday rotating mirror 1346 where it is reflected as verticallypolarized light. The vertically polarized light returns to the thirdport 1340 of the polarization beam splitter 1332 and is passed throughto the first port 1330.

As discussed above, the light which was originally incident at the firstport 1330 of the polarization beam splitter 1332 was oriented atapproximately 45° to the horizontal and vertical polarizations. Thus,approximately 50 percent of the light corresponding to the verticallypolarized component of the light passed through the polarization beamsplitter 1332 to the third port 1340 and thus to the common delay fiber1342. The vertically polarized light propagates through the delay loop1344 and is reflected by the Faraday rotating mirror 1346 ashorizontally polarized light. The reflected horizontally polarized lightpasses through the delay loop 1344 and back to the third port 1340 ofthe polarization beam splitter 1332. Because the light is horizontallypolarized, the light is reflected to the fourth port 1348 of thepolarization beam splitter 1332 and is thus caused to propagate via thefirst array input/output fiber 714, through the second horizontalpolarizer 1350, through the depolarizer 1352 and into the array 716 topropagate therein in a counterclockwise direction. The depolarizer 1352assures that the counterclockwise propagating light has components inall polarizations so that when the counterclockwise propagating lightemerges from the array 716, there will be at least a portion of thelight in the horizontal polarization.

The counterclockwise propagating light emerges from the array 716 viathe second array input/output fiber 720, and the horizontally polarizedcomponent of the light passes through the first horizontal polarizer1338, which eliminates the light at other polarization orientations. Thehorizontally polarized light resulting from the counterclockwisepropagating portion of the light enters the second port 1336 of thepolarization beam splitter 1332 and is reflected to the first port 1330of the polarization beam splitter 1332 where it is combined with thevertically polarized light which resulted from the clockwise propagatingportion of the light.

The combined light propagates to the fourth port of the coupler 1316where approximately 50 percent of the combined light is coupled to thesecond port of the coupler 1316 and thus to the detector subsystem 1360via the fiber 1362. The 1×n coupler 1364 divides the light into Nportions. For example, in FIG. 30, N is equal to 2, and a first portionof the light is coupled to the polarization controller 1366 to propagatethrough the polarizer 1368 to the first detector 1370, and a secondportion of the light is coupled to the polarization controller 1372 topropagate through the polarizer 1374 to the second detector 1376. Theorientations of the polarization controllers 1366, 1372 and thepolarizers 1368, 1374 can be adjusted to bias the optical signalsincident on the first detector 1370 and the second detector 1376 atdifferent phases. For example, the signal applied to the second detector1376 can be biased to be in quadrature with the signal applied to thefirst detector 1370 so that when one signal has minimum sensitivity, theother signal has maximum sensitivity, and vice versa.

As discussed above, each of the two signal portions travels the samedistance through the array 716, through the common delay fiber 1342, andthrough the delay loop 1344. Thus, in the absence of perturbationscaused by acoustic signals or other noise impinging on the sensors inthe array 716, the two portions will be in phase and will constructivelyinterfere to generate a combined optical signal having a linearpolarization of 45°; however, the light has a state of polarizationorthogonal to the original state of polarization. Thus, if the originalstate of polarization was +45°, then the state of polarization of theoutput light (again in the absence of a phase perturbation) is −45°.

In the presence of an acoustic signal, the clockwise propagating lightand the counterclockwise propagating light experience a relative phaseshift. With increasing relative phase shift, the state of polarizationof the two interfering beams changes from −45° linear polarization toleft-hand circular polarization to +45° polarization to right-handcircular polarization and back to −45° polarization. The progressionthrough these four states of polarization define a circle on thePoincaré sphere. The state of polarization at the output of thepolarization beam splitter 1332 corresponds to a point along this circleon the Poincaré sphere whose location on the circle is a function of theacoustically-induced non-reciprocal phase shift.

After traveling from the output of the polarization beam splitter 1332,through the common array input/output fiber 1334, through the coupler1316 and to the detector subsystem 1360, the state of polarization ofthe combined signal is altered arbitrarily by the unknown birefringenceof the fiber 1334. The polarization controller 1366 proximate topolarizer 1368 in front of the first detector 1370 and the polarizationcontroller 1372 proximate to the polarizer 1374 in front of the seconddetector 1376 are used to re-orient the states of polarization to arespective selected state of polarization for each detector 1370, 1376.The polarization controllers 1366, 1372 are set, for example, when noacoustic signals are applied to the array 716, and thus no relativephase shift is introduced to the counterpropagating optical signals.

For example, to provide a bias point of ±90° for the first detector1370, the polarization controller 1376 is set so that when the combinedlight at the output of the polarization beam splitter 1332 has aleft-hand circular state of polarization, the first detector 1370detects either a maximum intensity or a minimum intensity of the light.For other states of polarization of the output light, the first detector1370 detects light having an intensity between the maximum intensity andthe minimum intensity.

As a further example, the second detector 1376 can advantageously be setto a different bias point, such as, for example, 0° and 180°. For thisbias point, the polarization controller 1372 is set so that, when thelight at the output of the polarization beam splitter 1332 has a −45°state of polarization, the second detector 1376 detects either a maximumintensity or a minimum intensity of the light. For other states ofpolarization of the output light, the second detector 1376 detects lighthaving an intensity between the maximum intensity and the minimumintensity.

It should be understood that the light applied to the input of thepolarization beam splitter 1332 can have a state of polarization otherthan +45°. For example, if the input light has an original left-handcircular state of polarization, the polarization controllers 1366, 1372are set accordingly to provide the appropriate bias points to the firstdetector 1370 and the second detector 1376.

FIG. 31 illustrates an alternative configuration of a folded Sagnacsensor array 1300′, which is substantially similar to the folded Sagnacsensor array 1300 of FIG. 30. In the folded Sagnac sensor array 1300′ ofFIG. 31, the depolarizer 1352 is located in the second arrayinput/output fiber 720 rather than in the first array input/output fiber714. Because of the reciprocal structure of the sensor array 716, therelocation of the depolarizer 1352 to the fiber 720 does not change theoverall operation of the folded Sagnac sensor array 1300′ with respectto the operation of the folded Sagnac sensor array 1300. The operationof the folded Sagnac sensor array 1300′ is similar to the operation ofthe folded Sagnac sensor array 1300 and will not be described in detailherein.

FIG. 32 illustrates a further alternative embodiment of a folded Sagnacacoustic sensor array 1400, which is similar to the folded Sagnac sensorarray 1300 of FIG. 30, and like elements have been numbered accordingly.Unlike the folded Sagnac sensor array 1300, the folded Sagnac sensorarray 1400 replaces the 2×2 coupler 1316 with a polarization independentoptical circulator 1410. The optical circulator performs a similarfunction as the 2×2 coupler 1316; however, in the folded Sagnac sensorarray 1300 approximately 50 percent of the input light is lost when theinput light is split at the coupler 1316 and approximately 50 percent ofthe output light is lost when it is split at the coupler 1316. In theembodiment 1400, substantially all the input light is passed from thepolarized SFS 1310 through the circulator 1410 to the polarization beamsplitter 1332 and substantially all the output light is passed from thepolarization beam splitter 1332 through the circulator 1410 to thedetector subsystem 1360.

FIG. 33 illustrates an alternative configuration of a folded Sagnacsensor array 1400′, which is substantially similar to the folded Sagnacsensor array 1400 of FIG. 32. In the folded Sagnac sensor array 1400′ ofFIG. 33, the depolarizer 1352 is located in the second arrayinput/output fiber 720 rather than in the first array input/output fiber714. Because of the reciprocal structure of the sensor array 716, therelocation of the depolarizer 1352 to the fiber 720 does not change theoverall operation of the embodiment 1400′ with respect to the operationof the folded Sagnac sensor array 1400. Thus, the operation of thefolded Sagnac sensor array 1400′ will not be described in detail herein.

FIG. 34 illustrates a further alternative embodiment of a folded Sagnacsensor array 1600 in accordance with the present invention, whichincludes a combined input/output subsystem 1610 which is coupled to thearray 716 in a manner similar to the manner described above inconnection with FIGS. 30-33.

In FIG. 34, a polarized source 1620 provides linearly polarized inputlight along an axis of a polarization maintaining fiber 1622. Thepolarization maintaining fiber 1622 is rotated such that thepolarization axis is oriented at ±45° with respect to the verticalpolarization axis of the input output system 1610. The light from thefiber 1622 is coupled to the input/output subsystem 1610 via a firstcollimating lens 1630. The first collimating lens 1630 directs the lighttoward a first port 1634 of a first polarization beam splitter (PBS)1632, which also has a second port 1636, a third port 1638 and a fourthport 1640. The second port 1636 directs a portion of the input lighttoward a first 45° Faraday rotator (45° FR) 1642. The third port 1638directs a portion of the input light toward a second 45° Faraday rotator1644. As will be described below, the fourth port 1640 directs aselected portion of output light to a detection subsystem 1650.

The light passing through the first Faraday rotator 1642 is collimatedby a second collimating lens 1660 and is coupled into the arrayinput/output fiber 720 and thus propagates to the sensor portion of thearray 716 to propagate in a clockwise direction therein.

The light passing through the second Faraday rotator 1644 passes througha half-wave (λ/2) plate 1662. The half-wave plate 1662 has first andsecond birefringent axes (not shown). One of the birefringent axes isoriented at an angle of 22.5° with respect to the vertical polarizationaxis of the incoming light and at −22.5° with respect to the 45°polarization of the light traveling toward it from the source (i.e., theaxis lies between vertical and the polarization of the light). Thepurpose of this orientation will be described below. The light passingthrough the half-wave plate 1662 enters a first port 1672 of a secondpolarization beam splitter 1670, which also has a second port 1674, athird port 1676 and a fourth port 1678. As discussed below, the secondport 1674 is not coupled to additional elements. Light output from thethird port 1676 is directed toward a third collimating lens 1680. Lightoutput from the from the fourth port 1678 is directed toward a fourthcollimating lens 1682.

The light passing through the fourth collimating lens 1682 is coupledinto the first array input/output fiber 714 and passes through thedepolarizer 1352 into the sensor portion of the array 716 to propagatein a counterclockwise direction therein.

The light passing through the third collimating lens 1680 is focusedonto the end of the common delay fiber 1342, propagates through thedelay loop 1344 to the Faraday rotating mirror 1346, back through thedelay loop 1344 and back to the collimating lens 1680. The reflectedlight is thus directed back into the third port 1676 of the secondpolarization beam splitter 1670.

As discussed above, the light from the fourth port 1640 of the firstpolarization beam splitter 1632 enters the detection subsystem 1650. Thedetection subsystem 1650 comprises a first beam splitter 1690, a secondbeam splitter 1692, a first birefringent element 1694, a secondbirefringent element 1696, a first detector 1698, a second detector1700, a first polarizer 1702, and a second polarizer 1704. A firstpercentage of the light from the fourth port 1640 is reflected by thefirst beam splitter 1690 and passes through the first birefringentelement 1694 and the first polarizer 1702 to the first detector 1698.The remaining portion of the light from the fourth port 1640 passesthrough the first beam splitter 1690 and is incident on the second beamsplitter 1692 where a second percentage of the light is reflected by thesecond beam splitter 1692 to pass through the second birefringentelement 1696 and the second polarizer 1704 to the second detector 1700.The remaining portion of the light passes through the second beamsplitter 1692 to additional elements (not shown). If only two detectorsare provided, the first percentage of coupling is advantageously 50percent and the second percentage is advantageously 100 percent so thatboth detectors 1698, 1700 receive approximately the same amount oflight. If a third detector (not shown) is included, then the firstpercentage is advantageously about 33⅓ percent, and the secondpercentage is advantageously about 50 percent so that the seconddetector 1700 detector also receives approximately 33⅓ percent of theoriginal light. The third detector would then receive the remaining 33⅓percent.

The folded Sagnac sensor array 1600 of FIG. 34 operates in the followingmanner. As discussed above, the light incident on the first lens 1630 isoriented at 45° to the vertical and horizontal polarization axes. Thus,the light passing through the lens 1630 and entering the first port 1634of the first polarization beam splitter 1632 has a component in thehorizontal state of polarization and a component in the vertical stateof polarization. The horizontal component is reflected by thepolarization beam splitter 1632 to the second port 1636, and thevertical component is passed through the polarization beam splitter 1632to the third port 1638.

The horizontal component from the second port 1636 passes through thefirst Faraday rotator 1642, and the state of polarization is rotated by45° in a first direction (e.g., clockwise) so that the light emergingfrom the first Faraday rotator 1642 and incident on the second lens 1660has a linear state of polarization at 45°. The light passes through thesecond lens 1660 and enters the second array input/output fiber 720 topropagate in the clockwise direction through the array 716. The lightmay encounter changes in polarization within the array 716. Thus, asdescribed above, the light exiting the array 716 via the first arrayinput/output fiber 714 passes through the depolarizer 1352, whichassures that at least a portion of the light is in the horizontal andvertical states of polarization.

The clockwise propagating light from the first array input/output fiber714 enters the input/output subsystem 1610 via the fourth lens 1682 andis incident on the second polarization beam splitter 1670. The verticalcomponent of the light passes through the second polarization beamsplitter 1670 is output from the second port 1674 and is discarded. Thehorizontally polarized component of the light is reflected to the thirdport 1676 of the second polarization beam splitter 1670 and passesthrough the third lens 1680 to the common delay fiber 1342 to cause thelight to propagate through the delay loop 1344, be reflected by theFaraday rotating mirror 1346 in the vertical state of polarization, passback through the delay loop 1344 and the common delay fiber 1342 to thethird lens 1680. The reflected light in the vertical state ofpolarization passes from the third port 1676 to the first port 1672 ofthe second polarization beam splitter 1670, passes through the half-waveplate 1662 to the second Faraday rotator 1644 to the third port 1638 ofthe first polarization beam splitter 1632. Because the half-wave plate1662 is oriented with one of its birefringent axes at 22.5° with respectto the vertical polarization axis, the vertical light incident on thehalf-wave plate 1662 is caused to be mirrored about the birefringentaxis so that the state of polarization of the light emerging from thehalf-wave plate 1662 is oriented at 45° with respect to vertical andhorizontal axes. The second Faraday rotator 1644 rotates the state ofpolarization by a further 45° to cause the light emerging from thesecond Faraday rotator 1644 and incident on the third port 1638 of thefirst polarization beam splitter 1632 to have a horizontal state ofpolarization. Thus, the light entering the third port 1638 is reflectedto the fourth port 1640 and enters the detection subsystem 1650 in thehorizontal state of polarization.

As set forth above, the vertical component of the input light incidenton the first port 1634 of the first polarization beam splitter 1632passes through to the third port 1638. The state of polarization of thelight is rotated by 45° by the second Faraday rotator 1644 to a 45°state of polarization with respect to the vertical and horizontalpolarization axes. The state of polarization of the light is thenmirrored about the birefringent axis of the half-wave plate 1662 so thatthe state of polarization of the light emerging from the half-wave plateis again oriented in the vertical direction. It will be understood byone skilled in the art that the non-reciprocal action of the secondFaraday rotator 1644 causes the vertically polarized light that passesfrom left to right through the second Faraday rotator 1644 and thenthrough the half-wave plate 1646 to first be rotated to a 45° state ofpolarization and then to be mirrored back to a vertical state ofpolarization. In contrast, the vertically polarized light that passesfrom right to left is first mirrored by the half-wave plate 1646 to a45° state of polarization and is then rotated by the second Faradayrotator 1644 to a horizontal state of polarization.

The vertically polarized light from the half-wave plate 1662 enters thefirst port 1672 of the second polarization beam splitter 1670 and passesthrough to the third port 1676 to the third lens 1680. The verticallypolarized light passes through the common delay fiber 1342, through thedelay loop 1344, to the Faraday rotating mirror 1346, and is reflectedback through the delay loop 1344 and the common delay fiber 1342 ashorizontally polarized light. The horizontally polarized light passesthrough the third lens 1680 to the third port 1676 of the polarizationbeam splitter 1670. The horizontally polarized light is reflected to thefourth port 1678 and passes through the fourth lens to 1682 to the firstarray input/output fiber 714 and through the depolarizer 1352 topropagate in a counterclockwise direction through the array 716.

The counterclockwise propagating light emerges from the array 716 viathe second array input/output fiber 720 and passes through the secondlens 1660 to the first Faraday rotator 1642. The first Faraday rotator1642 rotates state of polarization of the light by 45°. Since the lightwas effectively depolarized by the depolarizer 1352, the light thatpasses through the first Faraday rotator 1642 to the second port 1634 ofthe first polarization beam splitter 1632 includes light that hashorizontally and vertically polarized components. The horizontallypolarized components of the light are reflected to the first port 1634and are output through the first lens 1630 to the input fiber 1622. Anisolator (not shown) is advantageously included to absorb the light.

The vertically polarized components of the counterclockwise propagatinglight entering the second port 1636 of the first polarization beamsplitter 1632 pass to the fourth port 1640 and are combined with thehorizontally polarized components of the clockwise propagating light. Asdiscussed above in connection with FIG. 30, if the counterpropagatinglight experiences no relative phase shift, the light is combined aslinearly polarized light at a 45° state of polarization. A relativephase shift causes the state of polarization to vary, as furtherdiscussed above.

The birefringent elements 1694, 1696 are included to selectively biasthe light incident on the detectors 1698, 1900 by introducing a relativephase shift for the light in the two different polarizations (e.g., thehorizontal and vertical polarizations, the +45° and −45° polarizations,or the left-hand circular and right-hand circular polarizations). Thebirefringent elements may advantageously comprise linear or circularwaveplates (e.g., quarter-wave plates, half-wave plates, Faradayrotators, or the like).

FIG. 35 illustrates an embodiment of a folded Sagnac acoustic sensorarray 1750 similar to the folded Sagnac acoustic sensor array 1600 ofFIG. 34, and like elements are identified with the same numbers as inFIG. 34. Unlike the embodiment of FIG. 34, the folded Sagnac acousticsensor array 1750 includes an unpolarized light source 1720 instead ofthe polarized light source 1620. In order to utilize the unpolarizedlight source 1720, the folded Sagnac acoustic sensor array 1750 includesa 45° polarizer 1730 between the first collimating lens 1630 and thefirst polarization beam splitter 1632. The 45° polarizer 1730 causes thelight incident on the first port 1634 of the first polarization beamsplitter 1632 to be oriented at 45° and to thus have substantially equalcomponents in the horizontal and vertical polarizations. Thus, thefolded Sagnac acoustic sensor array 1750 of FIG. 35 operates insubstantially the same manner as the folded Sagnac acoustic sensor array1600 of FIG. 34, and the operation of the folded Sagnac acoustic sensorarray 1750 will not be described in further detail.

FIG. 36 illustrates a further embodiment of a folded Sagnac acousticsensor array 1800 similar to the folded Sagnac acoustic sensor arrays1600 and 1750 of FIGS. 34 and 35, respectively, and like elements areidentified with the same numbers as in FIGS. 34 and 35. Unlike theembodiments of FIGS. 34 and 35, the folded Sagnac acoustic sensor array1800, the light signals passing through the polarizers 1702 and 1704 arenot directed to the detectors 1698 and 1700. Rather, the folded Sagnacacoustic sensor array 1800 includes a collimating lens 1810 positionedproximate to the polarizer 1702 and a collimating lens 1812 positionedproximate to the polarizer 1704. The collimating lens 1810 directs thelight from the polarizer 1702 into a first end 1822 of a fiber 1820. Thefiber 1820 has a second end 1824 proximate to the first detector 1698such that the light entering the fiber 1820 from the collimating lens1810 is incident on the first detector 1698. Similarly, the collimatinglens 1812 directs the light from the polarizer 1702 into a first end1832 of a fiber 1830. The fiber 1830 has a second end 1834 proximate tothe second detector 1700 such that the light entering the fiber 1830from the collimating lens 1812 is incident on the second detector 1700.By including the collimating lenses 1810 and 1812 and the fibers 1820and 1830, the fibers are able to transport the light for a distance tothe detectors 1698 and 1700 so that the detectors may be located inremote locations proximate to the detection electronics (not shown).

Note that in FIGS. 34, 35 and 36, the depolarizer 1352 can be relocatedfrom the first array input/output fiber 714 to the second arrayinput/output fiber 720 without significantly affecting the operatingcharacteristics of the folded Sagnac acoustic sensor array 1600, thefolded Sagnac acoustic sensor array 1750 or the folded Sagnac sensorarray 1800.

In the foregoing embodiments of FIGS. 17-36, an amplified sensor array716 receives two counterpropagating signals from and returns twoperturbed counterpropagating signals to a respective front end systemwhich generates the signals and detects the perturbations. In theforegoing embodiments, the sensor array 716 is depicted as a ladderstructure with the sensors 722(i) in respective rungs 718(i). Aplurality of erbium doped fiber amplifiers (EDFAs) 724 distributedwithin the ladder structure to amplify the signals distributed to andreceived from the sensors 722(i).

FIG. 37 illustrates an alternative embodiment of a 16-sensor array 2000for use in combination with the front end systems described in FIGS.30-36 in place of the sensor array 716. In particular, the sensor array2000 is interposed between a first input/output fiber 2002 and a secondinput/output fiber 2004. The first input/output fiber 2002 correspondsto the first input/output fiber 714 in FIGS. 30-36 for example, and thesecond input/output fiber 2004 corresponds to the second input/outputfiber 720 in FIGS. 30-36 Thus, light entering the sensor array 2000 viathe first input/output fiber 2002 propagates from right to left throughthe sensor array 2000 in FIG. 37 and exits via the second input/outputfiber 2004, and the light entering the sensor array 2000 via the secondinput/output fiber 2004 propagates from left to right through the sensorarray 2000 and exits via the first input/output fiber 2002. Thus, theright-to-left propagating light in FIG. 37 corresponds to thecounterclockwise propagating light in FIGS. 30-36, and the left-to-rightpropagating light in FIG. 37 corresponds to the clockwise propagatinglight in FIGS. 30-36.

The sensor array 2000 is implemented as a tree structure that comprisesan outer layer 2010 of amplifiers 2012(1), 2012(2), each having an gainof g₁. Each outer layer amplifier 2012(i) is followed by a respectivesplitter 2014(1), 2012(2).

The sensor array 2000 further comprises an inner layer 2020 ofamplifiers 2022(1), 2022(2), 2022(3), 2022(4), 2022(5), 2022(6),2022(7), 2024(8), each having a respective gain g₂. Each inner layeramplifier 2022(i) is followed by a respective splitter 2024(1), 2024(2),2024(3), 2024(4), 2024(5), 2024(6), 2024(7), 2024(8).

In the illustrated preferred embodiments, each of the splitters 2014(i),2024(i) is advantageously a 4×4 splitter having two sets of fourinput/output ports, with one set of ports at each end of the splitter.Light entering the splitter via one of the input/output ports at one endexits the four input/output ports at the other end in substantiallyequal portions. In FIG. 37, each splitter 2024(i) in the inner layer2020 is illustrated with four ports at a first end and one port at asecond end. It should be understood that the three unused ports (notshown) at the second end are non-reflectively terminated. Thus, lightentering the single used port at the second end is split between thefour ports at the first end, and approximately one-fourth of the lightentering each of the four ports at the first end is coupled to thesingle used port at the second end. The remaining three-fourths of thelight from each of the ports at the first end is lost via the threeunused ports. Thus, each of the 4×4 splitters 2024(i) operates as a1-to-4 splitter for light propagating in one direction and as a 4-to-1combiner for light propagating in the opposite direction.

Because of the layout illustrated in FIG. 37, each layer, 2010, 2020comprises a respective left portion 2010L, 2020L and a respective rightportion 2010R, 2020R.

The right portion 2010R of the outer layer 2010 comprises the amplifier2012(1) followed by the splitter 2014(1). The left portion 2010L of theouter layer 2010 comprises the amplifier 2012(2) followed by thesplitter 2014(2).

The right portion 2020R of the inner layer 2020 comprises the amplifiers2022(1), 2022(2), 2022(3), 2022(4) and the splitters 2024(1), 2024(2),2024(4). The left portion 2020L of the inner layer 2020 comprises theamplifiers 2022(5), 2022(6), 2022(7), 2022(8) and the splitters 2024(5),2024(6), 2024(7), 2024(8).

The outer layer 2010 and the inner layer 2020 are symmetrical about asensor layer 2030 that comprises a plurality of sensors 2032(1) . . .2032(16). The sensors 2032(1) . . . 2032(16) are organized as a fourgroups 2040(1) . . . 2040(4) of sensors.

Each of the four sensors 2032(1) . . . 2032(4) in the first group2040(1) are connected between a respective one of the four input/outputports at the first end of the splitter 2024(1) and a respective one ofthe four input/output ports at the first end of the splitter 2024(5).

Each of the four sensors 2032(5) . . . 2032(8) in the second group2040(2) is connected between a respective one of the four input/outputports at the first end of the splitter 2024(2) and a respective one ofthe four input/output ports at the first end of the splitter 2024(6).

Each of the four sensors 2032(9) . . . 2032(12) in the third group2040(3) is connected between a respective one of the four input/outputports at the first end of the splitter 2024(3) and a respective one ofthe four input/output ports at the first end of the splitter 2024(7).

Each of the four sensors 2032(13) . . . 2032(16) in the fourth group2040(4) is connected between a respective one of the four input/outputports at the first end of the splitter 2024(4) and a respective one ofthe four input/output ports at the first end of the splitter 2024(8).

Within each group 2040(i) of sensors 2032(i), three of the sensorsinclude delay fibers 2042(i) interposed in the path between theinput/output ports of the two splitters 2024(i) coupled to the group.The delay fibers 2042(i) have suitable lengths selected to provideproper timing of the time-division multiplexed (TDM) pulses passingthrough the sensors 2032(i).

The first sensor in each group (i.e., the sensors 2032(1), 2032(5),2032(9), 2032(13)) has no additional delay other than the inherentpropagation delay of the path between the respective two splitters.

The second sensor in each group (i.e., the sensors 2032(2), 2032(6),2032(10), 2032(14)) has an additional delay of τ provided by first andsecond delay fibers 2042(1), each having delay of τ/2.

The third sensor in each group (i.e., the sensors 2032(3), 2032(7),2032(11), 2032(15)) has an additional delay of 2τ provided by third andfourth delay fibers 2042(2), each having delay of τ.

The fourth sensor in each group (i.e., the sensors 2032(4), 2032(8),2032(12), 2032(16)) has an additional delay of 3τ provided by a fifthand sixth delay fibers 2042(3), each having delay of 3τ/2.

The single used input/output port at the second end of each of thesplitters 2024(1) . . . 2024(8) is coupled to a first input/outputterminal of a respective one of the amplifiers 2022(1) . . . 2022(8). Inthe illustrated embodiment, the splitter 2024(1) is coupled to theamplifier 2022(1), the splitter 2024(2) is coupled to the amplifier2022(2), and so on.

A second input/output port of each of the amplifiers 2022(1) . . .2022(4) is coupled to a respective one of the four input/output ports ofthe splitter 2014(1). A second input/output port of each of theamplifiers 2022(5) . . . 2022(8) is coupled to a respective one of thefour input/output ports of the splitter 2014(2).

The amplifiers 2022(1), 2022(5) are coupled to the splitters 2014(1),2014(2) without any additional delay interposed in the path other thaninherent propagation delay.

The amplifiers 2022(2), 2022(6) are coupled to the splitters 2014(1),2014(2) via respective delay fibers 2044(1), which each provide anadditional delay of 2τ. Thus, the total additional delay to all thesensors in the second group 2040(2) is 4τ.

The amplifiers 2022(3), 2022(7) are coupled to the splitters 2014(1),2014(2) via respective delay fibers 2044(2), which each provide anadditional delay of 4τ. Thus, the total additional delay to all thesensors in the third group 2040(3) is 8τ.

The amplifiers 2022(4), 2022(8) are coupled to the splitters 2014(1),2014(2) via respective delay fibers 2044(3), which each provide anadditional delay of 6τ. Thus, the total additional delay to all thesensors in the fourth group 2040(4) is 12τ.

It can be readily determined that the total additional delay between thesplitters 2014(1) and 2014(2) through the first sensor 2032(1) is 0τ.The total additional delays between the splitters 2014(1) and 2014(2)through the remaining sensors 2032(2) . . . 2032(16) are as follows:

2023(2) τ 2032(3) 2τ 2023(4) 3τ 2032(5) 4τ 2023(6) 5τ 2032(7) 6τ 2023(8)7τ 2032(9) 8τ 2023(10) 9τ 2032(11) 10τ 2023(12) 11τ 2032(13) 12τ2023(14) 13τ 2032(15) 14τ 2023(16) 15τ

As discussed above, the value of τ is selected to be sufficient toadequately separate the pulses in time after propagating through eachsensor 2032(i). For example, in one embodiment, the value of τ isadvantageously 60 nanoseconds so that 50-nanosecond pulses are separatedby 10-nanosecond guardbands.

As shown in FIG. 37, the outer layer 4×4 splitters 2014(1) and 2014(2)are configured with the four input/output ports on the first enddirected toward the sensors 2032(i) and connected as described above.Two of the input/output ports on the second end are non-reflectivelyterminated and are not shown in FIG. 37. A first input/output port ofthe amplifier 2012(1) is coupled to one of the remaining input/outputports on the second end of the splitter 2014(1). A first input/outputport of the amplifier 2012(2) is coupled to one of the remaininginput/output ports on the second end of the splitter 2014(2). A secondremaining input/output port on the second end of the splitter 2014(1) iscoupled to a first pump source 2050(1). A second remaining input/outputport on the second end of the splitter 2014(2) is coupled to a secondpump source 2050(2).

A second input/output port of the amplifier 2012(1) is coupled to thefirst input/output fiber 2002 via a first input/output port at a firstend of a first wavelength division multiplexing (WDM) coupler 2060(1). Asecond input/output port at a second end of the first WDM coupler2060(1) is coupled to the first input/output fiber 2002. A thirdinput/output port, also at the second end of the first WDM coupler2060(1), is coupled to receive pump light from a third pump source2062(1).

The first WDM coupler 2060(1) is configured so that all the lightentering the first input/output port at a signal wavelength λ_(s) (e.g.,1,560 nanometers) passes through the coupler to the second input/outputport. Similarly, light at the signal wavelength λ_(s) entering thesecond input/output port passes through the coupler to the firstinput/output port.

Since substantially no coupling occurs at the signal wavelength λ_(s),substantially all the light entering the array 2000 via the firstinput/output fiber 2002 passes through the first WDM coupler 2060(1)from right to left in FIG. 37 and enters the amplifier 2012(1).Similarly, substantially all the light exiting the array 2000 from leftto right via the amplifier 2012(1) passes through the first WDM coupler2060(1) to the first input/output fiber 2002.

In contrast to the light at the signal wavelength, substantially all thepump light at a pump wavelength λ_(p), (e.g., 1,480 nanometers) from thepump source 2062(1) enters the third input/output port of the first WVDMcoupler 2060(1) and is coupled to the second input/output port of thecoupler to thereby propagate to the amplifier 2012(1).

Similarly, a second input/output port of the amplifier 2012(2) iscoupled to the second input/output fiber 2004 via a first input/outputport at a first end of a second wavelength division multiplexing (WDM)coupler 2060(2). A second input/output port at a second end of thesecond WDM coupler 2060(2) is coupled to the second input/output fiber2004. A third input/output port, also at the second end of the secondWDM coupler 2060(2), is coupled to receive pump light from a fourth pumpsource 2062(2).

As described above for the first WDM coupler 2060(1), substantially allthe light at the signal wavelengths λ_(s) passes through the second WDMcoupler 2060(2) without coupling, and substantially all the light fromthe pump source 2062(2) is coupled to the second input/output port ofthe second WDM coupler 2060(2) and propagates to the amplifier 2012(2).

Unlike the two WDM couplers 2060(1) and 2060(2), the 4×4 splitters2014(1), 2014(2) and the splitters 2024(1) . . . 2024(8) are broadbandcouplers. That is, the splitters have substantially the same couplingcoefficients over the range of wavelengths from 1,480 nanometers to1,560 nanometers. Thus, the light at the signal wavelengths λ_(s) andthe light at the pump wavelengths λ_(p) are coupled substantially thesame. In addition, the splitters are preferably characterized so thatcoupling coefficients are substantially similar so that the light at thesignal wavelength λ_(s) and the light at the pump wavelength λ_(p)entering the array 2000 are divided substantially equally among the fouroutputs of each splitter and so that the light at the signal wavelengthfrom each sensor are combined substantially equally.

In operation, the pump source 2062(1) provides pump light to the firstWDM coupler 2060(1). The pump light propagates to the amplifier 2012(1)and is absorbed therein so that the amplifier provides gain to lightentering the amplifier at the signal wavelength. Any residual pump lightnot absorbed by the amplifier 2012(1) passes through the splitter2014(1) and is divided substantially equally among the four amplifiers2022(1) . . . 2022(4). In addition, light from the pump source 2050(1)enters the splitter 2014(1) and is divided substantially equally topropagate to the four amplifiers 2022(1) . . . 2022(4). The pump lightis absorbed in the four amplifiers to provide amplification for thesignal light.

In like manner, light from the pump source 2062(2) is coupled to theamplifier 2012(2) via the second WDM coupler 2060(2), and any residualpump light is split among the four amplifiers 2024(5) . . . 2024(8) viathe splitter 2014(2). In addition, pump light from the pump source2050(2) is split by the splitter 2014(2) and provided to the fouramplifiers 2024(1) . . . 2024(4).

A signal light pulse having a duration of approximately 50 nanosecondsenters the array 2000 via the first input/output fiber 2002 and is firstamplified by the amplifier 2012(1). The amplified light pulse is thensplit into four substantially equal portions by the splitter 2014(1).The first portion is amplified by the amplifier 2022(1) and then entersthe first sensor group 2040(1) via the splitter 2024(1) which dividesthe first portion into four subportions. A first subportion passesthrough the sensor 2032(1). A second subportion passes through thesensor 2032(2). A third subportion passes through the sensor 2032(3). Afourth subportion passes through the sensor 2032(4).

The first subportion encounters no additional delay. The secondsubportion is delayed by the two τ/2 delay fibers 2042(1) for a totaldelay of τ. The third subportion is delayed by the two τ delay fibers2042(2) for a total delay of 2τ. The fourth subportion is delayed by thetwo 3τ/2 delay fibers 2042(3) for a total delay of 3τ.

The four subportions in the first sensor group 2040(1) are recombined bythe splitter 2024(5). However, because of the respective propagationtimes through the sensors and the delays, the first subportion arrivesat the splitter approximately 60 nanoseconds before the secondsubportion. The second subportion arrives approximately 60 nanosecondsbefore the third subportion. The third subportion arrives approximately60 nanoseconds before the fourth subportion. Thus, the output of thesplitter 2024(5) comprises four 50-nanosecond pulses at 60-nanosecondintervals.

The pulses output from the splitter 2024(5) are amplified by theamplifier 2022(5) and then propagate to the splitter 2014(2), where thepulses are combined with pulses that passed through the second sensorgroup 2040(2), the third sensor group 2040(3) and the fourth sensorgroup 2040(4).

Within the second, third and fourth sensor groups, the signal pulse issplit into four subportions, the four subportions are selectivelydelayed and are then recombined as discussed above with respect to thefirst sensor group. In addition, the portion of light entering thesecond sensor group 2040(2) is first delayed by the first 2τ delay2044(1) before being amplified by the amplifier 2022(2) and split by thesplitter 2024(2). The signal pulses output by the second sensor group2040(2) are recombined by the splitter 2046(6) and amplified by theamplifier 2022(6) and are then delayed by a second 2τ delay 2044(1)before arriving at the splitter 2014(2). Because of the extra 4τ-delayprovided by the two 2τ delays 2044(1), the first signal pulse arrivingat the splitter 2014(2) from the second sensor group arrives 60nanoseconds after the fourth signal pulse from the first sensor group.

In like manner, the signal pulses outputted from the third sensor groupare delayed by an additional 8τ by the two delays 2044(2) so that thefirst signal pulse from the third sensor group arrives 60 nanosecondsafter the fourth signal pulse from the second sensor group.

Similarly, the signal pulses outputted from the fourth sensor group aredelayed by an additional 12τ by the two delays 2044(3) so that the firstsignal pulse from the fourth sensor group arrives 60 nanoseconds afterthe fourth signal pulse from the third sensor group.

The signal pulses from the four sensor groups are combined in thesplitter 2014(2). The signal pulses are then amplified by the amplifier2012(2) and are output from the array via the second input/output fiber2004.

In like manner, a 50-nanosecond light pulse entering the sensor array2000 via the second input output fiber 2004 and traveling from left toright through the array is first amplified by the amplifier 2012(2) andsplit by the splitter 2014(2) into four pulses. The four pulses areselectively delayed and further split, as discussed above for the rightto left propagating signal, so that sixteen 50-nanosecond pulses at60-nanosecond intervals arrive at the splitter 2014(1). The sixteenpulses are combined by the splitter 2014(1) and are amplified by theamplifier 2012(1) before exiting the array via the first input/outputfiber 2002.

It should be understood that the pulse widths and the delay times can bevaried. Furthermore, additional splitters, amplifiers and delays can beprovided to increase the number of sensors in the array 2000.

The amplifiers 2012(i), 2022(i) that precede each splitter 2014(i),2024(i) compensate for the splitting loss in the respective coupler. Thetree-structure telemetry of FIG. 37 has the advantage of directingunused pump power from earlier amplifiers into later amplifiers. In anexemplary system wherein the pump light has a pump wavelength λ_(p) of,for example, 1,480 nanometers, and wherein the signal has a signalwavelength λ_(s) of, for example, 1,550 nanometers, the splittingcouplers must be either dual band or broadband to have similar couplingratios across the pump and signal wavelengths. Such couplers arecommercially available from, for example, Gould Fiber Optics, 1121Benfield Boulevard, Millersville, Md. The embodiment of FIG. 37 has afurther advantage that, if additional pump power is required to powerthe later stage amplifiers, the additional pump power can be added atone of the unused ports of one of the 4×4 splitters 2012. Thus, only thefirst-stage amplifiers (i.e., the amplifiers 2012 in the outer layer2010) require additional WDM couplers.

As discussed above, the sensor telemetry must be constructed to provideincreasing differential delays between the sensor paths to preventpulses from overlapping. This time division multiplexing (TDM) scheme isaccomplished by placing delay fibers 2042(i), 2044(i) as shown in FIG.37, where each delay is a fraction of or a multiple of a time intervalτ. The time interval τ is the sensor interrogation window. The durationof each pulse (i.e., the pulse width) is selected to be less than τ toavoid overlap. The prescribed delays include the effect of the delaycoil and the fiber connecting the nearest components. For example, thelongest delay of 8τ is from the first splitting coupler 2014(1) to thesplitting coupler 2024(4). It can be easily seen that the total traveltime through each sensor increments by τ with these delays, inaccordance with TDM requirements.

Unlike previously described telemetries (i.e., array structures), thetree structure EDFA telemetry described above is advantageous because ofits ability to be pumped by a large number of smaller-power pump lasers(instead of a few larger-power pump lasers). In particular, the pumplight can be introduced at various levels in the tree structure ratherthan at one or two locations.

FIG. 38, comprising FIGS. 38A and 38B, illustrates an alternative sensorarray 2100 that includes 16 sensors 2110(i) in an amplifiedtree-structure telemetry using a polarization-based front end 2120similar to the system described above in connection with FIGS. 30-36,for example.

FIG. 38A illustrates a front end 2120 similar to the front end shown inFIG. 34, and like elements are numbered accordingly. Unlike the frontend shown in FIG. 34, the front end 2120 in FIG. 38A includes a bandpassinterference filter 2130 of approximately 10-nanometer width centered at1,550 nanometers positioned between the output of the polarizedsuperfluorescent source 1620 and the collimator 1630. The filter 2130narrows the light from the source 1620 in order to reduce anydeleterious impact of the wavelength-dependence of the Faraday rotators1642, 1644. An additional pair of bandpass interference filters 2132,2134 are positioned in the output path between the fourth port 1640 ofthe first polarization beam splitter 1632 beamsplitter 1634 and thefirst beamsplitter 1690.

The amplified array 2100 in FIG. 38B is similar to the amplifiedtree-structure design shown in FIG. 37. Again, the timing is based on60-nanosecond windows which accommodated 50-nanosecond pulses with10-nanosecond guardbands between pulses. The array timing is measuredduring the construction of the array 2100 to ensure that, not only isthe overall timing of the sensors correct, but also that each sensor2110(i) is placed symmetrically in the sensor telemetry, therebyensuring the same phase sensitivity at each sensor.

The first input/output fiber 2002, the second input/output fiber 2004,the third pump source 2062(1), the fourth pump source 2062(2), the firstWDM coupler 2060(1), the second WDM coupler 2060(2), the first amplifier2012(1) with the gain g₁, and the second amplifier 2012(2) with the gaing₁ operate as described above in connection with FIG. 37, and arelabeled accordingly.

Unlike the array 2000 of FIG. 37, the array 2100 of FIG. 38B utilizesfour levels of 2×2 broadband couplers 2140(i) instead of the two levelsof 4×4 broadband splitters 2014(i), 2024(i). The array 2100 again usestwo levels of amplification.

The counterclockwise propagating signal light incident to the array 2100via the first input/output fiber 2002 and any residual pump light fromthe amplifier 2012(i) passes through a first 2×2 coupler 2140(1), whichsplits the light into two substantially equal portions, a first portion(upper in FIG. 38B) and second portion (lower in FIG. 38B).

The first portion then enters a first input/output port of second 2×2coupler 2140(2), which splits the first portion of the light into twosubstantially equal portions, a third (upper) portion at a thirdinput/output port and a fourth (lower) portion at a fourth input/outputport. Similarly, the second portion then enters a first input/outputport of a third 2×2 coupler 2140(3), which splits the second portion ofthe light into two substantially equal portions, a fifth (upper) portionat a third input/output port and a sixth (lower) portion at a fourthinput/output port.

The third, fourth, fifth and sixth portions enter respective amplifiers2150(1), 2150(2), 2150(3), 2150(4), which are pumped by pump light fromthe first pump source 2050(1) via a 2×2 coupler 2160(1), which dividesthe pump light into two substantially equal portions. A first portion ofthe pump light is provided to a second input/output port of the secondcoupler 2140(2), and a second portion of the pump light is provided to asecond input/output port of the third coupler 2140(3). Thus, the pumplight is further divided by the couplers 2140(2), 2140(3) so thatapproximately one-fourth of the original pump light is provided to eachof the amplifiers 2150(1), 2150(2), 2150(3), 2150(4).

The amplified signal light from the amplifier 2150(1) is provided as aninput to a fourth 2×2 coupler 2140(4). The amplified signal light fromthe amplifier 2150(2) is provided as an input to a fifth 2×2 coupler2140(5). The amplified signal light from the amplifier 2150(3) isprovided as an input to a sixth 2×2 coupler 2140(6). The amplifiedsignal light from the amplifier 2150(4) is provided as an input to aseventh 2×2 coupler 2140(7).

The light incident on each coupler 2140(4), 2140(5), 2140(6), 2140(7) isdivided substantially equally into two portions.

A first portion of the light from the fourth coupler 2140(4) is providedas an input to the first sensor 2110(1) and the second sensor 2110(2)via an eighth 2×2 coupler 2140(8), which divides the first portion ofthe light into two substantially equal portions. A second portion of thelight from the fourth coupler 2140(4) is provided as an input to thethird sensor 2110(3) and the fourth sensor 2110(4) via a ninth 2×2coupler 2140(9), which divides the second portion of the light into twosubstantially equal portions. Thus, each of the sensors 2110(1),2110(2), 2110(3), 2110(4) receives approximately one-sixteenth of theoriginal signal light incident on the first coupler 2140(1).

Similarly, the fifth sensor 2110(5) and the sixth sensor 2110(6) eachreceive approximately one-sixteenth of the original signal light fromthe fifth coupler 2140(5) via a tenth 2×2 coupler 2140(1). The seventhsensor 2110(7) and the eighth sensor 2110(8) each receive approximatelyone-sixteenth of the original signal light from the fifth coupler2140(5) via an eleventh 2×2 coupler 2140(11).

The ninth sensor 2110(9) and the tenth sensor 2110(10) each receiveapproximately one-sixteenth of the original signal light from the sixthcoupler 2140(6) via a twelfth 2×2 coupler 2140(12). The eleventh sensor2110(11) and the twelfth sensor 2110(12) each receive approximatelyone-sixteenth of the original signal light from the sixth coupler2140(6) via a thirteenth 2×2 coupler 2140(13).

The thirteenth sensor 2110(13) and the fourteenth sensor 2110(14) eachreceive approximately one-sixteenth of the original signal light fromthe seventh coupler 2140(7) via a fourteenth 2×2 coupler 2140(14). Thefifteenth sensor 2110(15) and the sixteenth sensor 2110(16) each receiveapproximately one-sixteenth of the original signal light from theseventh coupler 2140(7) via a fifteenth 2×2 coupler 2140(15).

Note that at each level of division and amplification, couplers andamplifiers are numbered from the middle of the figure towards the bottomof the figure so that the lower-numbered couplers and amplifierspropagate the counterclockwise-propagating signal light to thelower-numbered sensors.

The light incident to the array 2100 via the second input/output fiber2004 is similarly divided by a plurality of couplers 2140(16) . . .2140(30) into sixteen substantially equal portions that are provided tothe sensors 2110(1) . . . 2110(16) in the opposite (i.e., clockwisedirection). Pump light from the second pump source 2050(2) is dividedsubstantially equally by a coupler 2160(2), and is further divided bythe couplers 2140(16) and 2140(17) so that four approximately equalportions of the pump light are provided as inputs to four amplifiers2150(5), 2150(6), 2150(7), 2150(8), which operate to amplify the signallight propagating in the clockwise direction. Note that the couplers andthe amplifiers are numbered from the middle toward the upper portion ofFIG. 38B so that the lower-numbered couplers and amplifiers provide thesignal light to the lower-numbered sensors for the clockwise-propagatingsignal light.

In order to provide the time division multiplexing characteristicsdescribed above in connection with FIG. 37, the embodiment of FIG. 38Bincludes a plurality of delay fibers in the array 2100. For example, afirst delay fiber 2180(1) is positioned between the coupler 2140(8) andthe second sensor 2110(2), and a second delay fiber 2180(1) ispositioned between the coupler 2140(23) and the second sensor 2110(2).Similar delay fibers 2180(1) are positioned between the sensors 2110(2),2110(4), 2110(6), 2110(8), 2110(10), 2110(12) 2110(14), 2110(1 couplers.Each delay fiber 2180(1) has a length of approximately 6 meters toprovide 30 nanoseconds of delay so that the light propagating throughthe second sensor in each pair of sensors is separated by approximately60 nanoseconds from the light propagating through the first sensor inthe pair of sensors.

The pairs of signals passing through the first sensor 2110(1) and thesecond sensor 2110(2) are separated from the pairs of signals passingthrough the third sensor 2110(3) and the fourth sensor 2110(4) by anadditional 120 nanoseconds of delay provided by a first delay fiber2180(2) in the path from the coupler 2140(4) to the coupler 2140(9) anda second delay fiber 2180(2) in the path from the coupler 2140(19) tothe coupler 2140(24). Each of the delay fibers 2180(2) has a delay ofapproximately 60 nanoseconds provided by a fiber having a length ofapproximately 12 meters.

In like manner, a first 120-nanosecond delay fiber 2180(3) is interposedin the path between the coupler 2140(5) and the coupler 2140(10), and asecond 120-nanosecond delay fiber 2180(3) is interposed in the pathbetween the coupler 2140(20) the coupler 2140(25). A first180-nanosecond delay fiber 2180(4) is interposed in the path between thecoupler 2140(5) and the coupler 2140(11), and a second 180-nanoseconddelay fiber 2180(4) is interposed in the path between the coupler2140(20) the coupler 2140(26).

For the remaining eight sensors, a first 225-nanosecond delay fiber2180(5) is interposed between the coupler 2140(1) and the coupler2140(3), and a second 225-nanosecond delay fiber 2180(5) is interposedbetween the coupler 2150(16) and the coupler 2140(18).

A first 15-nanosecond delay fiber 2140(6) is interposed between thecoupler 2140(6) and the coupler 2140(12), and a second 15-nanoseconddelay fiber 2140(6) is interposed between the coupler 2140(21) and thecoupler 2140(27) so that the total additional delay to the sensors2110(9) and 2110(10) is 480 nanoseconds. It should be understood thatthe equivalent delay could be provided by increasing the delay of thedelay fibers 2180(5) to 240 nanoseconds and not including the15-nanosecond delay fibers 2140(6).

A first 75-nanosecond delay fiber 2180(7) is interposed in the pathbetween the coupler 2140(6) and the coupler 2140(13), and a second75-nanosecond delay fiber 2180(7) is interposed in the path between thecoupler 2140(21) the coupler 2140(28). A first 135-nanosecond delayfiber 2180(8) is interposed in the path between the coupler 2140(7) andthe coupler 2140(14), and a second 135-nanosecond delay fiber 2180(7) isinterposed in the path between the coupler 2140(22) the coupler2140(29). A first 195-nanosecond delay fiber 2180(9) is interposed inthe path between the coupler 2140(7) and the coupler 2140(15), and asecond 195-nanosecond delay fiber 2180(9) is interposed in the pathbetween the coupler 2140(22) the coupler 2140(30).

Each of the 2×2 couplers is advantageously characterized beforeconstruction so that the couplers can be arranged in the array such thatthe transmission of each sensor path was as equal as possible. The tenamplifiers 2012(1), 2012(2), 2150(1) . . . 2150(8) in FIG. 38B arepreferably erbium-doped fiber amplifiers (EDFAs). In FIG. 38B, the gaing₁ of the amplifiers 2012(1), 2012(2) and the gain g₂ of the amplifiers2150(1) . . . 2150(8) are substantially the same. In the preferredembodiments, the length of the fiber in each amplifier is approximately1.45 meter to provide approximately 8.2 dB of gain.

The total gain seen by a given sensor path from the first input/outputfiber 2002 to the second input/output fiber 2004 in each direction is32.8 dB (8.2 dB×4 amplifiers). This roughly compensates for the totalsplitting loss (3.1 dB per coupler×8 couplers) plus the insertion lossof the WDM coupler and the excess loss due to splices and fiber bending.The transmission through an exemplary sensor array constructed inaccordance with FIG. 38B was measured to be a few dB greater than unity.The amplifiers 2012(1), 2012(2), 2150(1) . . . 2150(8) are pumped by thepump sources 2050(1), 2050(2), 2062(1), 2062(2), which advantageouslycomprise four low-power (e.g., 60 mW) 1,480-nm lasers. As discussedabove, the pump light from the pump sources 2062(1), 2062(2) are coupledinto the array 2100 through the WDM couplers 2060(1), 2060(2) to pumpthe first-stage amplifiers 2012(1), 2012(2). The second-stage amplifiers2150(1) . . . 2150(8) are pumped by the unused power from thefirst-stage amplifiers along with the power added by the pump sources2050(1), 2050(2). The pump power is sufficiently large (and the signalpower is sufficiently low) to pump saturate all of the amplifiers2012(1), 2012(2), 2150(1) . . . 2150(8). Preferably, the pump sources2050(1), 2050(2), 2062(1), 2062(2) are coupled to the array 2100 viaone-way isolators 2190 so that any light returning from the array 2100does not propagate to the pump sources and perturb the source.

In order to test the sensors 2110(i), an additional PZT transducerwrapped with fiber is placed in series with each sensor 2110(i) so thatthe effect of an acoustic signal can be simulated for each sensor byselectively activating one or more of the PZT transducers. (The PZTtransducers are for testing purposes only and are not part of thepreferred embodiment. Therefore, the PZT transducers are not shown inthe drawings.) FIGS. 39A and 39B illustrate the returned pulse trainmeasured at one of the detectors 1698, 1700 (FIG. 38A) for a pulse widthof 50 nanoseconds and repetition frequencies of 0.942 MHz(1.06-microsecond spacing) and 1.042 MHz (0.96-microsecond spacing),respectively. Each signal pulse applied to the array results in a trainof sixteen output pulses which return at 60-nanosecond intervals. Thus,at the lower repetition frequency (FIG. 39A), there are gaps between thetrains of 16 pulses that return from the sensor array. The power fromall of the sensors return with a variation of less than 1.7 dB. Thesepower variations result from the insertion-loss differences of thevarious couplers and splices.

FIG. 39B illustrates the optimal operating condition in which the pulsetrains repeat continuously (i.e., the sixteenth pulse of one pulse trainoccurs approximately 60 nanoseconds before the first pulse of the nextpulse train). The pulse timing was verified to be correct to within 1nanosecond (20 centimeters) for all the sensor paths.

In practice, demultiplexing is achieved post detection by digitalprocessing. In an experimental embodiment, a lithium niobate amplitudemodulator 2192 is connected to one of the input/output fibers (e.g., thesecond input/output fiber 2004 in FIG. 38A) and is used to pass a singlepulse corresponding to the signal to be demultiplexed. The lithiumniobate modulator 2192 provides clean switching, eliminating the largenoise associated with electronic switches used after the photodiodes.

To measure the polarization sensitivity of the array, a large phasemodulation was induced in the fifth sensor 2110(5) by activating therespective series connected PZT transducer for the fifth sensor. Thesignal from the fifth sensor was demultiplexed and its detected powerwas monitored on a digital scope. When a peak-to-peak phase modulationof greater than π is induced in the signal, the visibility of the sensorcan be measured by recording the maximum and minimum voltage obtained onthe digital scope. Such a trace is illustrated in FIG. 40A, yielding ameasured visibility of 0.9. As described above, a significant advantageof the polarization-based front-end design is that polarization-inducedsignal fading is passively eliminated, and the visibility should ideallybe unity on each sensor, independently of the loop birefringence.

To test this prediction, the visibility of the fifth sensor was measuredfor eight random settings of a common polarization controller 2194positioned in one of the input/output fibers (e.g., the firstinput/output fiber 2002 in FIG. 38A). The results of the measurementsare shown in FIG. 40B, which demonstrate the expected polarizationindependence. The visibility is short of unity due mainly to thepresence of distributed amplifier ASE power at the detector.

To illustrate the impact and importance of noise subtraction in aSagnac-based array using a broadband ASE source, a phase modulation at 4kHz was placed on one of the sensors using the respective fiber wrappedPZT to simulate an acoustic signal, and an amplitude modulation at 3.6kHz was placed on the optical signal from the source 1620 using alithium niobate (LiNbO₃) modulator 2196 in the signal path from thesource 1620 to the collimator 1630 in FIG. 38A. The detected responsewas measured in a single-detection configuration using only one detectorand a balanced-detection configuration in which the outputs of the twodetectors 1698, 1700 are connected in such a way that the DC currentsfrom the detectors are subtracted to remove the source amplitude noise.FIG. 41 illustrates the results of the measurements. For thesingle-detection configuration, the amplitude modulation at 3.6 kHz isclearly seen while the phase modulation at 4 kHz is not visible underthe broadband source amplitude noise (ASE—ASE beat noise). In thebalanced-detection configuration, the amplitude modulation is subtractedby over 50 dB. The phase modulation tone becomes easily visible due tothe reduced broadband noise (ASE—ASE beat noise subtraction) and theincreased phase modulation signal. As illustrated, thebalanced-detection configuration achieves a large degree of noisesubtraction.

The optical noise floor of a Sagnac-based array with distributed opticalamplifiers is a function of the optical power at the receiver. The noiseN_(s) ^(bd) for a balanced-detection receiver is given as:$\begin{matrix}{N_{S}^{bd} = {{4{B\left( {s + a} \right)}} + {\frac{8B}{\Delta \quad v_{sa}}({sa})} + {\frac{16B}{\Delta \quad v_{a}}\left( {a_{cw}a_{ccw}} \right)}}} & (19)\end{matrix}$

where B is the detection bandwidth; where s, a_(cw), a_(ccw) are thesource, clockwise distributed amplifier, and counterclockwisedistributed amplifier ASE photon fluxes (#/s); and wherea=a_(cw)+a_(ccw). In Equation 19, Δv_(a) is the bandwidth of thedistributed amplifier photons arriving at the detectors caused byspontaneous emission in the distributed amplifiers, as defined by:$\begin{matrix}{{\Delta \quad v_{a}} = \frac{\left( {\int{{P\left( v_{a} \right)}{v_{a}}}} \right)^{2}}{\int{\left( {P\left( v_{a} \right)} \right)^{2}{v_{a}}}}} & (20)\end{matrix}$

where P(v_(a))dv_(a) is the optical power in the frequency band v_(a) tov_(a)+dv_(a).

Further in Equation 19, Δv_(sa) is defined by: $\begin{matrix}{\frac{1}{\Delta \quad v_{sa}} = {\frac{1}{\Delta \quad v_{s + a}}\left\lbrack {1 + {\frac{s}{2a}\left( {1 - \frac{\Delta \quad v_{s + a}}{\Delta \quad v_{s}}} \right)} + {\frac{a}{2s}\left( {1 - \frac{\Delta \quad v_{s + a}}{\Delta \quad v_{a}}} \right)}} \right\rbrack}} & (21)\end{matrix}$

where Δv_(s) is the bandwidth of source photons arriving at the detectorand Δv_(s+a) is the combined source and distributed amplifier photonsarriving at the receiver defined in accordance with Equation 20 bysubstituting v_(a) and v_(s+a) for v_(a).

The noise expression of Equation 19 is given as a photon countfluctuation per unit bandwidth. Equation 19 can be rewritten as followsfor a detector current noise, (ΔI²): $\begin{matrix}{{\langle{\Delta \quad I^{2}}\rangle} = {{4{{Be}\left( {I_{s} + I_{a}} \right)}} + {\frac{8B}{\Delta \quad v_{sa}}\left( {I_{s}I_{a}} \right)} + {\frac{16B}{\Delta \quad v_{a}}{\left( {I_{acw}I_{accw}} \right).}}}} & (22)\end{matrix}$

where B is the detection bandwidth.

Equation 22 is confirmed experimentally by measuring the optical noiseof the 16-sensor array 2100. An excess noise factor F is first definedas the ratio of the detected noise to the shot-noise limit and is givenby: $\begin{matrix}{{F \equiv \frac{\langle{\Delta \quad I^{2}}\rangle}{4{{eB}\left( {I_{s} + I_{a}} \right)}}} = {1 + {\frac{2}{e\quad \Delta \quad v_{sa}}\left( \frac{I_{s}I_{a}}{I_{s} + I_{a}} \right)} + {\frac{4}{e\quad \Delta \quad v_{a}}{\left( \frac{I_{acw}I_{accw}}{I_{s} + I_{a}} \right).}}}} & (23)\end{matrix}$

By measuring the detector currents I_(s), I_(acw) and I_(accw) alongwith their optical spectra (to get Δv_(sa) and Δv_(a)), the expectedexcess noise factor can be calculated and compared with the measurednoise. To measure distributed amplifier currents and spectra (I_(acw),I_(accw), and Δv_(a)), the source is turned off and the output is simplymeasured. Since the distributed amplifier power cannot be switched off,the source current/spectrum is found by measuring the total outputcurrent/spectrum and subtracting the distributed amplifiercurrent/spectrum. This can be done because the distributed amplifiersare pump saturated and their ASE output is not affected by the presenceof the small source power in the array (experimentally verified).

With these output currents and spectra measured, the predicted excessnoise factor can be calculated through Equation 23. To measure the noisefactor, the output noise is measured directly along with the DC receivercurrents. The shot noise level is then calculated from the DC currents,and the excess noise factor is found by taking the ratio of the measurednoise and the shot noise. The predicted and measured excess noisefactors can then be compared.

In an experimental setup, the measurements are done for a series of DCsource currents at the receiver. These currents are controlled byadjusting the source power with the DC bias of the LiNbO₃ modulator2196. For each DC source current, the measurements are taken for sensorscounts of 4, 8, 12, and 16. The sensor counts were controlled bydarkening one, two, or three 4-sensor subarrays by inducing large bendlosses in the fibers addressing these subarrays. As the sensor count ischanged, the pulse repetition frequency is adjusted to keep a continuoustrain of pulses at the output. Because of this, the average source powerat the detector is the same for each sensor count within a given set ofmeasurements, and the only change is the amount of distributed amplifierASE. Therefore, by taking measurements as a function of the sensorcount, the noise trend as a function of the ASE power can be found andcompared with the predictions of Equation 23.

The results of the foregoing measurements are shown in FIGS. 42A, 42B,42C for three DC source receiver currents. As shown in FIG. 42A for thelowest DC source current (1.4 μA), the agreement between the predictedand measured noise is excellent. For larger source powers and thereforelarger DC source receiver currents (4.05 μA in FIG. 42B and 9.89 μA inFIG. 42C), the measured noise rises above the predicted noise as thesource current is increased. This is due, at least in part, to thepick-up of acoustic room noise in the array, a persistent problem in thetesting of the Sagnac acoustic sensor designs. As the source power isincreased, the sensitivity of the array increases, and room noise pickedup in the fiber making up the sensor array and timing coils becomesnon-negligible. Although this noise makes measurements of the opticalnoise floor difficult in the laboratory, it is not a troublesome sourceof noise in practice since this acoustic pick-up is, in realapplications, the signal which is intended to be measured.

In the Mach-Zehnder based arrays, the sensor response is determined bythe length of fiber wrapped around the hydrophone and by the normalizedhydrophone responsivity, both of which cannot be changed once the sensoris constructed. In the Sagnac sensor arrays described above, theresponse is determined in part by these same parameters and is alsodetermined in part by the length of the delay loop, which is sharedamongst all the sensors. Because this delay loop can be remotely locatedaway from the sensors, where it is accessible, it is feasible to changethe length of the delay loop. It was realized early in the developmentof Sagnac-based arrays that this ability to affect the response of allof the sensors by changing a single parameter, namely the delay looplength, could be of great benefit to the dynamic range performance ofthe array. For example, FIGS. 7, 19 and 21 illustrate embodiments havingtwo delay loops of different lengths to extend the dynamic range of thesensor array.

To illustrate how the response can be tuned by changing the delay looplength, consider a Sagnac array with an effective delay coil time of 50microseconds (e.g., approximately a delay loop length of approximately10 kilometers in an unfolded configuration or a delay loop length ofapproximately 5 kilometers in a folded configuration), a hydrophoneresponsivity, R, of −135 dB rad/μPa, and a flat optical noise floor of 1μrad/{square root over (Hz)}. The minimum detectable acoustic signal,a_(min)(f), can be expressed as: $\begin{matrix}{{a_{\min}(f)} = \frac{\varphi \sqrt{T}}{2\quad R\quad {\sin \left( {\pi \quad f\quad T_{d}} \right)}}} & (24)\end{matrix}$

where φ is the phase noise per unit bandwidth, T is the integrationtime, R is the hydrophone responsivity (rad/μPa), f is the frequency,and T_(d) is the delay coil time. Using the foregoing numbers, theminimum detectable acoustic signal in a 1-second integration time isplotted in FIG. 43A. Also plotted in FIG. 43A is the maximum detectablesignal defined as that which produces a 1 rad/{square root over (Hz)}phase signal. (This is only an approximation of the upper limit, but itserves as a way to compare the relative performance of twoconfigurations).

In FIG. 43B, the same curves are plotted for an array with the sameparameters except that the delay time through the delay loop has beenreduced from 50 microseconds to 0.5 microsecond. The shorter delay timemoves the detection range upward, making the detection range moresuitable for a noisy environment. In contrast, the longer delay time ismore suitable for a quieter environment.

To obtain the combined dynamic range of FIGS. 43A and 43B, an arraywould have to be constructed with both a short delay loop and a longdelay loop. An array that includes two delay coils, designed such thatthe pulses returning from all the sensors and the delay coils areseparable (in time or wavelength), can be easily constructed, asillustrated above in FIGS. 7, 19 and 21. In the illustratedconfigurations, each sensor returns two signals, each of which has adetection ranges corresponding to those of FIGS. 43A and 43B. In a quietenvironment, the signal from the longer delay loop would be used. In thepresence of a large acoustic signal that saturates the response from thelonger delay loop, the output from the shorter delay loop would be used.Delay coil switching such as this would allow the array to function inenvironments in which the magnitude of acoustic signal variesdramatically over time. However, the most useful implementation ofmultiple delay loops in a Sagnac-based array would combine the signalsfrom each delay loop instead of switching between them. By combining thesignals to produce a single output signal which has a true detectionrange given by the union of the detection ranges of the signals fromeach delay loop, the array could simultaneously detect down to the noisefloor with the longer delay loop in one frequency regime, while handlingan extremely large acoustic signal (which saturates the longer delayloop) in another frequency regime. To do this, an algorithm is used thatreceives as inputs a respective output signal from a given sensor fromeach delay loop and that returns a signal which has a detection rangethat comprises the union of the detection ranges of the output signalfrom the short delay loop and the output signal from the long delayloop.

If the signal from the shorter delay loop were the same as that from thelonger delay loop but with a reduced scale factor, then such analgorithm would be straightforward. However, because the two signalshave different frequency responses and sample the acoustic signal atdifferent times, the two signals have little obvious correlation. Toillustrate this, consider an acoustic signal which produces thehydrophone phase variation with time shown in FIG. 44A. FIGS. 44B and44C are plots of the induced phase difference due to the phase signal ofFIG. 44A for the signal which travels through the longer delay loop(FIG. 44B) and the shorter delay loop (FIG. 44C). The correlation is notobvious, and it is not clear how to combine these two signals to producea single signal.

FIGS. 45A and 45B illustrate an acoustic sensor array 2200 which workswith a linear extrapolation algorithm that combines the signals from twodelay loops to reconstruct the phase signal with a dynamic rangedramatically increased over that possible with a single delay loop.FIGS. 45A and 45B correspond to FIGS. 38A and 38B; however, in FIG. 45A,the front end portion of the system includes a first signal source 2210that outputs an first sequence of optical signal pulses at a firstwavelength λ₁ and includes a second signal source 2212 that outputs asecond sequence of optical signal pulses at a second wavelength λ₂. Forexample, the first wavelength λ₁ is advantageously approximately 1,520nanometers, and the second wavelength λ₂ is advantageously approximately1,550 nanometers. The outputs of the two signal sources, 2210, 2212 arecombined via a wavelength division multiplexing (WDM) coupler 2214 sothat a single signal stream comprising the two sequences of signalpulses is applied to the collimator 1630.

The system of FIGS. 45A and 45B further includes a first detectionsubsystem 2220 and a second detection subsystem 2222 to replace thesingle detection subsystem 1650 in FIG. 38A. For example, the twodetection subsystems 2220, 2222 are advantageously coupled to the outputof the beamsplitter 1632 via a collimator 2224, an optical fiber 2226,and a WDM coupler 2230 that couples the optical signals at the firstwavelength λ₁ and the second wavelength λ₂ to a first optical fiber 2232and a second optical fiber 2234, respectively. The first optical fiber2232 propagates the light at the first wavelength λ₁ to the firstdetection subsystem 2220 via a collimator 2236. The second optical fiber2234 propagates the light at the second wavelength λ₂ to the seconddetection subsystem 2222 via a collimator 2238.

The system of FIGS. 45A and 45B further differs in that the systemincludes two delay loops 2240 and 2242 in the delay path from thecollimator 1680 to the Faraday rotating mirror (FRM) 1346 in place ofthe single delay loop 1344 of FIG. 38A. In particular, the first delayloop 2240 is a longer delay loop that provides a delay of approximately100 microseconds, and the second delay loop 2242 is a shorter delay loopthat provides a delay of approximately 100 nanoseconds. As illustrated,the two delay loops 2240, 2242 are in parallel and are coupled to thecommon delay fiber 1342 via a first WDM coupler 2244 and a second WDMcoupler 2246. Thus, light from the third port 1676 of the polarizationbeam splitter 1670 passes through the collimator 1680 and enters thefirst WDM coupler 2244. The first WDM coupler 2244 directs light fromthe first source 2210 at the first wavelength λ₁ to the first delay loop2240 and directs light from the second source 2212 at the secondwavelength λ₂ to the second delay loop 2242. The two portions arerecombined in the second WDM coupler 2246, and the combined portions areincident on the Faraday rotating mirror 1346. The reflected signalportions are again split in the second WDM coupler 2246 such that theportion at the wavelength λ₁ again passes through the first (longer)delay loop 2240 and the portion at the wavelength λ₂ passes through thesecond (shorter) delay loop 2242. Thus, the light incident on the thirdport 1676 of the polarization beam splitter 1670 has a first pulse atthe first wavelength λ₁ that has been delayed by a longer delay than asecond pulse at the second wavelength λ₂. The two EM couplers 2244, 2246have fairly broad passbands which direct light to the two separate delaycoils 2240, 2242. Thus, for example, light from 1,525 nanometers to1,535 nanometers, which includes the first wavelength λ₁, isadvantageously directed to the longer delay loop 2240 while light from1,545 nanometers to 1,555 nanometers, which includes the secondwavelength λ₂, is advantageously directed to the shorter delay loop2242.

From the foregoing description, it can be seen that the signals receivedby the first detection subsystem 2220 provide a first acoustic dynamicrange determined by the delay through the first (longer) delay loop2240, and the signals received by the second detection subsystem 2222provide a second acoustic dynamic range determined by the delay throughthe second (shorter) delay loop 2242. It is assumed that all othercomponents of the system of FIGS. 45A and 45B are broadband to actsufficiently the same across all wavelengths used.

The phase modulation, φ(t), at the hydrophone in a given sensor isdirectly proportional to the acoustic signal amplitude. In aSagnac-based TDM array, each sensor is sampled with a repetition period,τ, and the returned signal is a phase difference, φ(t)−φ(t−T_(d)), whereT_(d) is the coil delay. Thus, for a single-delay coil sensor, aspreviously described, the discrete sampling of a given sensor yields thesamples S_(i) given by: $\begin{matrix}{\quad {\vdots \begin{matrix}{S_{- 1} = {{\varphi \left( {- \tau} \right)} - {{\varphi \left( {{- \tau} - T_{d}} \right)}\quad {mod}\quad 2\quad \pi}}} \\{S_{0} = {{\varphi (0)} - {{\varphi \left( {- T_{d}} \right)}\quad {mod}\quad 2\pi}}} \\{S_{1} = {{\varphi (\tau)} - {{\varphi \left( {\tau - T_{d}} \right)}\quad {mod}\quad 2\pi}}}\end{matrix}\quad \vdots \quad {S_{i} = {{\varphi \left( {i\quad \tau} \right)} - {{\varphi \left( {{i\quad \tau} - T_{d}} \right)}\quad {mod}\quad 2\pi}}}}} & (25)\end{matrix}$

from which the phase difference is recovered through

φ(iτ)−φ(iτ−T _(d))=S _(i)+2πn _(i)  (26)

where S_(i) has been measured by the array and n_(i) is the fringe countwhich is calculated through a fringe count algorithm, F, which takes asinput the current and previous samples,

n _(i) =F(S _(i) , S _(i−1), . . . )  (27)

The standard fringe count algorithm, F, for Mach-Zehnder and Sagnacarrays with a single delay coil is given by: $\begin{matrix}{n_{i} = \left\{ \begin{matrix}{n_{i - 1},} & {{{S_{i} - S_{i - 1}}} \leq \pi} \\{{n_{i - 1} + 1},} & {\left( {S_{i} - S_{i - 1}} \right) < {- \pi}} \\{{n_{i - 1} - 1},} & {\left( {S_{i} - S_{i - 1}} \right) > \pi}\end{matrix} \right.} & (28)\end{matrix}$

and is based on minimizing the phase difference between samples. Itgreatly extends the dynamic range of the sensor beyond that of a sensorwhich does not use fringe counting. The embodiment of FIGS. 45A and 45Benables the use of an improved algorithm which uses information from theadditional, shorter delay loop 2242 in the Sagnac-based TDM array toextend the dynamic range beyond the capabilities of the standard fringecount algorithm of Equation 28.

As discussed above, the Sagnac-based array in FIGS. 45A and 45Bincorporates multiple delay coils 2240, 2242. The broadband WDM couplers2244, 2246 direct fairly broad passbands to the separate delay loops2240, 2242. Thus, in the disclosed embodiment, for example, light in therange from 1,525 nanometers to 1,535 nanometers is advantageouslydirected to the longer delay loop 2240 while light in the range from1,545 nanometers to 1,555 nanometers is advantageously directed to theshorter delay loop 2242.

With multiple delay loops on separate wavelengths, the received dataset, S^(λ1), for a given sensor and the first wavelength λ₁ through thefirst delay loop 2240 is given by:

S ⁻¹ ^(λ1)=φ(−τ)−φ(−τ−T _(d) ^(λ1))mod2π

S ₀ ^(λ1)=φ(0)−φ(0−T _(d) ^(λ1))mod2π

S ₁ ^(λ1)=φ(τ)−φ(τ−T _(d) ^(λ1))mod2π  (29)

S _(i) ^(λ1)=100 (iτ)−φ(iτ−T _(d) ^(λ1))mod2π

The data set, S^(λ2), of the same sensor and the second wavelength λ₁through the second delay loop 2242 is given by:

S ⁻¹ ^(λ2)=φ(−τ+δ)−φ(−τ+δ−T _(d) ^(λ2))mod2π

S ₀ ^(λ2)=φ(δ)−φ(δ−T _(d) ^(λ2))mod2π

S ₁ ^(λ2)=φ(τ+δ)−φ(τ+δ−T _(d) ^(λ2))mod2π  (30)

S _(i) ^(λ2)=φ(iτ+δ)−φ(iτ+δ−T _(d) ^(λ2))mod2π

where δ describes the relative phase between the input (and thusreturned) signals from the wavelengths λ₁ and λ₂. If δ=0, the inputpulses at the two wavelengths are coincident in the input fiber. As theinput pulse at λ₂ is delayed relative to that at the λ₁, δ increases.The delay, δ, can be advantageously set to an arbitrary value based onthe needs of the algorithm.

For the algorithm to work, a constraint is imposed on the longer delayloop time, T_(d) ^(λ1), such that:

T _(d) ^(λ1) =pτ  (31)

where p is an integer such that the longer delay loop time is anintegral multiple of τ. The purpose of this constraint will be discussedbelow.

Based on the foregoing information, a new fringe counting algorithm, F′,is developed to use the information from multiple delay loops toincrease the dynamic range over the capability of the old, single delaycoil algorithm F defined by Equation 28. The algorithm is developed byplotting the signal from a given sensor from the longer delay loop 2240,β(t)=φ(t)−φ(t−T_(d) ^(λ1)) in FIG. 46A. FIG. 46A depicts the standardfringe counting algorithm, F. The array measures the samples S₁ and S₂,and because S₁=S₂>π, 2π is added to the sample S₂ to recover β(2τ), inaccordance with the algorithm of Equation 28. The plot in FIG. 46Aindicates that the fringe count of β(2τ) could be calculated with moreaccuracy if the slope of β(t) were known. If, for example, the slope,β′(t), of β(t) is measured at a frequency 1/τ. If these measurements arephased to give β′(t) at t=τ, t=2τ, etc., the fringe count of β(2τ) canthen be predicted by using a linear extrapolation form β(τ), i.e.,:$\begin{matrix}{{{\beta \left( {2\tau} \right)} \approx {{\beta (\tau)} + {\frac{\tau}{2}\left\lbrack {\beta^{\prime}\left( {2\tau} \right)} \right\rbrack}}},} & (32)\end{matrix}$

as shown in FIG. 46B. The measurement S₂=β(2τ) mod 2τ gives a more exactvalue for the sub-2τ part of β(2τ), but the linear extrapolation ofEquation 32 is more accurate to determine the fringe count as shownbelow in Equation 35 for the new algorithm F′.

It is also reasonable to use the slope at the midpoint (t=1.5τ) insteadof at the endpoints (t=τ, 2τ). In fact, doing this yields substantiallysimilar results with slightly higher low-frequency dynamic range andslightly lower high frequency dynamic range. Since the algorithm hasdramatically increased the low-frequency dynamic range almost to thepoint of not needing any more, we have chosen to use the endpointalgorithm to give a little more high-frequency dynamic range.

In order to complete the algorithm, β′ is measured. This is accomplishedthrough the use of the second, smaller delay loop 2242. The process ofmeasuring β′ starts with an expression for β′(iτ), $\begin{matrix}\begin{matrix}{{\beta^{\prime}\left( {i\quad \tau} \right)} \approx \quad \frac{{\beta \left( {i\quad \tau} \right)} - {\beta \left( {{i\quad \tau} - \Delta} \right)}}{\Delta}} \\{\approx \quad \frac{{\varphi \left( {i\quad \tau} \right)} - {\varphi \left( {{i\quad \tau} - T_{d}^{\lambda \quad 1}} \right)} - {\varphi \left( {{i\quad \tau} - \Delta} \right)} + {\varphi \left( {{i\quad \tau} - \Delta - T_{d}^{\lambda \quad 1}} \right)}}{\Delta}} \\{\approx \quad \frac{\left\lbrack {{\varphi \left( {i\quad \tau} \right)} - {\varphi \left( {{i\quad \tau} - \Delta} \right)}} \right\rbrack - \left\lbrack {{\varphi \left( {{i\quad \tau} - T_{d}^{\lambda 1}} \right)} - {\varphi \left( {{i\quad \tau} - {\Delta \quad T_{d}^{\lambda \quad 1}}} \right)}} \right\rbrack}{\Delta}}\end{matrix} & (33)\end{matrix}$

By identifying Δ as being equal to T_(d) ^(λ2), then the first andsecond bracketed terms in the final line of Equation 33 are recognizedas measurements of the second delay coil with δ=0 in Equation 30.

Writing (φ(iτ)−(φ(iτ−T_(d) ^(λ2)))=S_(i) ^(λ2)+2πn₁ ^(λ2) and(φ(iτ−T_(d) ^(λ1))−(φ(iτ−T_(d) ^(λ1)−T_(d) ^(λ2)))=S_(1−p)^(λ2)+2πS_(i−p) ^(λ2) (see Equation 34) gives: $\begin{matrix}{{\beta^{\prime}\left( {i\quad \tau} \right)} \approx \frac{\left( {S_{u}^{\lambda_{2}} + {2\pi \quad n_{i}^{\lambda_{2}}}} \right) - \left( {S_{i - p}^{\lambda_{2}} - {2\pi \quad n_{i - p}^{\lambda_{2}}}} \right)}{T_{d}^{\lambda \quad 2}} \equiv \beta_{i}^{\prime}} & (34)\end{matrix}$

The standard algorithm of Equation 28 is used to calculate the fringecount, n^(λ2), on the shorter delay loop 2240. β′ is then measured fromthe shorter delay loop 2240, and with this information, Equation 32 isused to measure the fringe count of each sample from the longer delayloop 2242.

In mathematical terms a parameter γ_(1−i) is defined asγ_(i−1)=τ/2(β_(i−1)+β_(i)), and the new linear extrapolation algorithm,F′, is given by: $\begin{matrix}{n_{i}^{\lambda 1}\left\{ \begin{matrix}{{n_{i - 1}^{\lambda \quad 1} + {{Int}\left( {\left\lbrack {S_{i - 1} + \gamma_{i - 1}} \right\rbrack/\left( {2\quad \pi} \right)} \right)}},} & {\left( {{\left\lbrack {S_{i - 1} + \gamma_{i - 1}} \right\rbrack \quad {mod}\quad 2\pi} - S_{i}} \right) \leq \pi} \\{n_{i - 1}^{\lambda \quad 1} + {{Int}\left( {{{\left\lbrack {S_{i - 1} + \gamma_{i - 1}} \right\rbrack/\left( {2\quad \pi} \right)} - 1},} \right.}} & {\left( {{\left\lbrack {S_{i - 1} + \gamma_{i - 1}} \right\rbrack \quad {mod}\quad 2\pi} - S_{i}} \right) < {- \pi}} \\{{n_{i - 1}^{\lambda \quad 1} + {{Int}\left( {\left\lbrack {S_{i - 1} + \gamma_{i - 1}} \right\rbrack/\left( {2\quad \pi} \right)} \right)} + 1},} & {\left( {{\left\lbrack {S_{i - 1} + \gamma_{i - 1}} \right\rbrack \quad {mod}\quad 2\pi} - S_{i}} \right) > \pi}\end{matrix} \right.} & (35)\end{matrix}$

where Int(x) returns the largest integer less than or equal to x.

The limitations of the standard fringe counting algorithm, F, ofEquation 28 and the new two-delay loop linear extrapolation algorithm,F′, of Equation 35 are compared below An accurate analysis of thedynamic range of an interferometric sensor array is quite involved anddifficult to reduce to a single figure of merit. For simplicity, thefollowing comparison compares the maximum tonal phase amplitude at thehydrophone as a function of frequency such that no fringe count errorsoccur. This analysis will illustrate the relative performance of the twoalgorithms as well and the two sensor types (Mach-Zehnder and Sagnac).

Beginning with the standard fringe counting algorithm, F, it can beshown that a fringe count error occurs if the magnitude of thedifference between two samples exceeds π,

|β(t)−β(t−τ)|≧π  (36)

where, as before, β(t)=φ(t)−φ(t−T_(d)) for a Sagnac sensor and β(t)=φ(t)for a Mach-Zehnder sensor. Writing φ(t)=A sin(2πft+θ), then for aMach-Zehnder array, the standard fringe counting algorithm, F,accurately demodulates a tonal phase signal given by A sin(2πft+θ) if:$\begin{matrix}{A < \frac{\pi}{2{{\sin \left( {\pi \quad \tau} \right)}}}} & (37)\end{matrix}$

and for a single delay coil Sagnac array if $\begin{matrix}{A < \frac{\pi}{4{{{\sin \left( {\pi \quad f\quad T_{d}} \right)}{\sin \left( {\pi \quad f\quad \tau} \right)}}}}} & (38)\end{matrix}$

The analysis of the new linear extrapolation algorithm, F′, is moreinvolved. Two conditions must be met for the algorithm to correctlycalculate the fringe count:

1. No fringe count errors on the shorter delay loop 2240. The shorterdelay loop uses the standard fringe counting algorithm, F to calculateits fringe counts, as if there were no longer delay loop in the system.Because of the short length of the delay loop 2240, the standard fringecounting algorithm, F, performs substantially better on it than it doeson the longer delay loop.

2. Linear extrapolation accurately predicts the fringe count. The linearextrapolation algorithm, F′, assumes that the phase difference measuredby the long delay loop 2242 varies linearly. For large enough curvatureon the phase difference, this condition will not be met and thealgorithm will fail.

The requirement for the first condition is given by Equation with T_(d)replaced by T_(d) ^(λ2) as follows: $\begin{matrix}{A < \frac{\pi}{4{{{\sin \left( {\pi \quad f\quad T_{d}^{\lambda \quad 2}} \right)}{\sin \left( {\pi \quad f\quad \tau} \right)}}}}} & (39)\end{matrix}$

For the second condition, a fringe count error occurs when a thelinearly extrapolated value of the next phase difference, {circumflexover (β)}(t+τ), and the actual next phase difference, β(t+τ), differ bymore than π. To see when this occurs, the linearly extrapolatedprediction is written as {circumflex over(β)}(t+τ)=b(t)+τ/2(β′(t)+β′(t+τ)). Thus, the condition for no fringecount error is

π>|{circumflex over(β)}(t+τ)−β(t+τ)|=|β(t)+τ/2(β′(t+τ)+β′(t))−β(t+τ)|.  (40)

Writing β(t)=(A sin(2πft+θ)−A sin(2πf(t−T_(d) ^(λ1))+θ)) and simplifyinggives: $\begin{matrix}{A < \frac{\pi}{4{{\sin \left( {\pi \quad f\quad T_{d}^{\lambda \quad 1}} \right)}\left\lbrack {{\sin \left( {\pi \quad f\quad \tau} \right)} - {\pi \quad f\quad \tau \quad {\cos \left( {\pi \quad f\quad \tau} \right)}}} \right\rbrack}}} & (41)\end{matrix}$

where it is assumed that fT_(d) ^(λ2)<<1. This is always the case sinceT_(d) ^(λ2) is chosen to be very small in order to extend the dynamicrange upward.

For the linear extrapolation algorithm to function, Equation 39 andEquation 41 must be satisfied. For comparison, the maximum value of Abefore fringe count errors occur is plotted as a function of frequencyin FIG. 47 for T_(d) ^(λ1)=100 microseconds, τ=50 microseconds, andT_(d) ^(λ2)=100 nanoseconds. Also shown are the results of a numericsimulation of the dynamic range for a Mach-Zehnder array (plotted withthe points represented by □'s), a Sagnac array with a single delay loop(plotted with the points represented by ⋄'s), and a Sagnac array withtwo delay loops employing the linear extrapolation algorithm presentedabove (plotted with the points represented by Δ's).

Note in FIG. 47 that for the two-delay loop linear extrapolation method,the limit in the Equation 39 and the limit in the Equation 41 must bothbe satisfied. Thus, the simulation follows the minimum of the twocurves. At lower frequencies, the simulation follows the curve ofEquation 39, and at higher frequencies, the simulation follows the curveof Equation 41.

FIG. 47 shows that the linear extrapolation algorithm behaves aspredicted and dramatically increases the dynamic range (e.g., by 3orders of magnitude from DC to 400 Hz, and by 2 orders of magnitude upto 1 kHz) over that possible with a Mach-Zehnder array or a single delayloop Sagnac array.

The foregoing describes an algorithm for using two delay loops 2240,2242 and two wavelengths λ₁, λ₂ in a Sagnac-based array to increase thedynamic range of the sensor substantially above that possible with asingle delay coil, which is already much greater than that possible in aMach-Zehnder array. The improvement in maximum tonal amplitude that canbe detected is quite substantial, calculated to be 2 orders of magnitudeor more up to 1 kHz in one particular example. This algorithm isstraightforward to implement and requires very little additionalcomputation on the receiver side. Thus, there are very few barriers toits implementation if the dynamic range it affords is deemed necessary.

Note further that the foregoing embodiments have been described inconnection with superfluorescent light sources. One skilled in the artwill appreciate that other light sources (e.g., laser sources) may alsobe advantageously used.

Although the foregoing description of the array in accordance with thepresent invention has addressed underwater acoustic sensing, it shouldbe understood that the present invention can be used to sense anymeasurand which can be made to produce non-reciprocal phase modulationsin a fiber. If, for example, the hydrophones were replaced with analternative sensing device which responds to a different measurand, thearray would detect that measurand in the same manner as acoustic wavesare detected. The array of the present invention can be advantageouslyused to sense vibrations, intrusions, impacts, chemicals, temperature,liquid levels and strain. The array of the present invention may also beused to combine a number of different sensors located at either the sameplace or located in different places (e.g., for the detection of variousfaults at various points along the hull of a ship or a building). Otherexemplary applications include the detection and tracking of movingautomobiles on highways or airplanes on airstrips for traffic monitoringand control.

Although described above in connection with particular embodiments ofthe present invention, it should be understood the descriptions of theembodiments are illustrative of the invention and are not intended to belimiting. Various modifications and applications may occur to thoseskilled in the art without departing from the true spirit and scope ofthe invention as defined in the appended claims.

What is claimed is:
 1. A sensor array comprising: a source of light; afirst coupler which receives the light from the source, couples a firstportion of the light to a first coupler port, and couples a secondportion of the light to a second coupler port; and an interferometricloop having a first end couple to the first coupler port to receive thefirst portion of the light, and having a second end coupled to thesecond coupler port to receive the second portion of the light, theinterferometric loop propagating the first portion o the light in afirst direction to the second coupler port, and propagating the secondportion of the light in a second direction opposite the first directionto the first coupler port, the interferometric loop comprising: aplurality of sensors coupled in parallel between t e first end of theinterferometric loop and the second end of the interferometric loop, thesensors perturbing light passing therethrough in response to sensedparameters; a first plurality of couplers that distribute the firstportion of the light approximately equally to each of the sensors, andthat collect the second portion of the light from each of the sensorsand propagate the collected light to the first end of theinterferometric loop; a second plurality of couplers that distribute thesecond portion of the light approximately equally to each of thesensors, and that collect the first portion of the light from each ofthe sensors and propagate the collected light to the second end of theinterferometric loop; at least one first amplifier coupled between thefi St end of the interferometric loop and the first plurality ofcouplers; at least one second amplifier coupled between the second endof the interferometric loop and the second plurality of couplers; and aplurality of delay portions connected between the first and second endsof the interferometric loops and the sensors, the delay portions havingdelays selected so that the light passing through each sensor is delayedby a different amount than the light passing through the other sensors.2. The sensor array as defined in claim 1, wherein: the first pluralityof couplers further comprises a first distribution coupler that receivesthe first portion of the light from the first amplifier and distributesthe first portion of the light to a first plurality of internalamplifiers that are coupled between the first distribution coupler and afirst plurality of internal couplers and that collects the secondportion of the light from the first plurality of internal amplifiers andpropagates the second portion of the light to the first amplifier; thesecond plurality of couplers further comprises a second distributioncoupler that receives the second portion of the light from the secondamplifier and distributes the second portion of the light to a secondplurality of internal amplifiers that are coupled between the seconddistribution coupler and a second plurality of internal couplers, andthat collects the first portion of the light from the second pluralityof internal amplifiers and propagates the first portion of the light tothe second amplifier; the first plurality of internal couplersdistributes the first portion of the light to the plurality of sensorsand collects the second portion of the light from the plurality ofsensors; and the second plurality of internal couplers distributes thesecond portion of the light to the plurality of sensors and collects t efirst portion of the light from the plurality of sensors.
 3. The sensorarray as defined in claim 2, wherein the source of light is a broadbandsource.
 4. The sensor array as defined in claim 3, wherein e broadbandsource is a superfluorescent fiber source.
 5. The sensor array asdefined in claim 2, wherein the first and second amplifiers and thefirst and second plurality of internal amplifiers are erbium-doped fiberamplifiers.
 6. The sensor array as defined in claim 2, wherein the firstand second distribution couplers and the first and second pluralities ofinternal couplers comprise 4×4 couplers.
 7. The sensor array as definedin claim 1, wherein the sensed parameters are acoustic signals.
 8. Amethod of sensing a parameter comprising: propagating light from asource of light through an interferometric loop such that approximatelyequal portions of the light counter-propagate in first and seconddirections in the loop; amplifying and coupling the light propagating inthe first direction of the interferometric loop into a plurality ofsensors such that approximately equal portions of the light propagatingin the first direction are passed through each of the sensors;amplifying and coupling the light propagating in the second direction ofthe interferometric loop into the plurality of sensors such thatapproximately equal portions of the light propagating in the seconddirection are passed through each of the sensors; interfering the lightpropagating in the first direction with the light propagating in thesecond direction to generate a plurality of output signals responsive tolight passing through each sensor in the first and second directions,wherein each of the sensors perturbs light passing therethrough inresponse to a sensed parameter, and has a unique optical path lengthsuch that the light propagating in the first direction interferes withthe light propagating in the second direction at a unique time.
 9. Themethod as defined in claim 8, wherein the sensed parameter is anacoustic signal.